Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW - CONFORMAL TRANSFORMATION - HODOGRAPH METHOD”

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 91 records · Page 5

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗

Theoretical determination of nozzle admittances using a finite element approach

A finite element method is used to predict the admittances of axisymmetric nozzles. It is assumed that the flow in the nozzle is isentropic and irrotational, and the disturbances are small so that linear analyses apply. An approximate, two dimensional compressible model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear wave equation. The differential form of the acoustic equation is transformed to an integral equation using Galerkin's method, and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. A two-dimensional FEM using linear elements is used to solve the acoustic equation. A one dimensional FEM is also used to solve the reduced equation of Crocco, and the solution verifies the sufficiency of the boundary residual formulation. Comparison between computed admittances and experimental data is shown to be quite good.

Sigman, R. K.↗

A numerical method for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

In the present paper, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BEC) systems has been developed and evaluated. The equations of motion are transformed to a generalized curvilinear coordinate system from which the transformed equations are discretized using finite difference approximations in the transformed domain. The hybrid scheme and a central differencing plus artificial dissipation scheme are used to approximate the convection terms in the governing equations. Effects of these two schemes on the accuracy of numerical predictions are studied. Solutions of the finite difference equations are obtained iteratively by using a pressure-velocity correction algorithm, SIMPLE-C. Numerical examples of two- and three-dimensional, laminar and turbulent flow problems are employed to evaluate the accuracy and efficiency of the present numerical method.

Chen, Y. S.↗

Phase Inversion: Inferring Solar Subphotospheric Flow and Other Asphericity from the Distortion of Acoustic Waves

A method is proposed for analyzing an almost monochromatic train of waves propagating in a single direction in an inhomogeneous medium that is not otherwise changing in time. An effective phase is defined in terms of the Hilbert transform of the wave function, which is related, via the JWKB approximation, to the spatial variation of the background state against which the wave is propagating. The contaminating effect of interference between the truly monochromatic components of the train is eliminated using its propagation properties. Measurement errors, provided they are uncorrelated, are manifest as rapidly varying noise; although that noise can dominate the raw phase-processed signal, it can largely be removed by low-pass filtering. The intended purpose of the analysis is to determine the distortion of solar oscillations induced by horizontal structural variation and material flow. It should be possible to apply the method directly to sectoral modes. The horizontal phase distortion provides a measure of longitudinally averaged properties of the Sun in the vicinity of the equator, averaged also in radius down to the depth to which the modes penetrate. By combining such averages from different modes, the two-dimensional variation can be inferred by standard inversion techniques. After taking due account of horizontal refraction, it should be possible to apply the technique also to locally sectoral modes that propagate obliquely to the equator and thereby build a network of lateral averages at each radius, from which the full three-dimensional structure of the Sun can, in principle, be determined as an inverse Radon transform.

Gough, Douglas↗

Effects of boundary layer refraction and fuselage scattering on fuselage surface noise from advanced turboprop propellers

An acoustic disturbance's propagation through a boundary layer is discussed with a view to the analysis of the acoustic field generated by a propfan rotor incident to the fuselage of an aircraft. Applying the parallel flow assumption, the resulting partial differential equations are reduced to an ordinary acoustic pressure differential equation by means of the Fourier transform. The methods used for the solution of this equation include those of Frobenius and of analytic continuation; both yield exact solutions in series form. Two models of the aircraft fuselage-boundary layer system are considered, in the first of which the fuselage is replaced by a flat plate and the acoustic field is assumed to be two-dimensional, while in the second the fuselage is a cylinder in a fully three-dimensional acoustic field. It is shown that the boundary layer correction improves theory-data comparisons over simple application of a pressure-doubling rule at the fuselage.

Mcaninch, G. L.↗

A finite element approach for predicting nozzle admittances

A finite element method is used to predict the admittances of axisymmetric nozzles. It is assumed that the flow in the nozzle is isentropic and the disturbances are small so that linear analyses apply. An approximate, two dimensional compressible model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear wave equation. The differential form of the acoustic equation is transformed to an integral equation by using Galerkin's method, and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. The boundary conditions are described for both straight and curved sonic lines. A two dimensional FEM with linear elements is used to solve the acoustic equation. A one dimensional FEM is also used to solve the reduced equation of Crocco, and the solution verifies the sufficiency of the boundary residual formulation. Comparison between computed admittances and experimental data is shown to be quite good.

Sigman, R. K.↗

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗

Application of velocity filtering to optical-flow passive ranging

The performance of the velocity filtering method as applied to optical-flow passive ranging under real-world conditions is evaluated. The theory of the 3-D Fourier transform as applied to constant-speed moving points is reviewed, and the space-domain shift-and-add algorithm is derived from the general 3-D matched filtering formulation. The constant-speed algorithm is then modified to fit the actual speed encountered in the optical flow application, and the passband of that filter is found in terms of depth (sensor/object distance) so as to cover any given range of depths. Two algorithmic solutions for the problems associated with pixel interpolation and object expansion are developed, and experimental results are presented.

Barniv, Yair↗

Design and Analysis Tools for Supersonic Inlets

Computational tools are being developed for the design and analysis of supersonic inlets. The objective is to update existing tools and provide design and low-order aerodynamic analysis capability for advanced inlet concepts. The Inlet Tools effort includes aspects of creating an electronic database of inlet design information, a document describing inlet design and analysis methods, a geometry model for describing the shape of inlets, and computer tools that implement the geometry model and methods. The geometry model has a set of basic inlet shapes that include pitot, two-dimensional, axisymmetric, and stream-traced inlet shapes. The inlet model divides the inlet flow field into parts that facilitate the design and analysis methods. The inlet geometry model constructs the inlet surfaces through the generation and transformation of planar entities based on key inlet design factors. Future efforts will focus on developing the inlet geometry model, the inlet design and analysis methods, a Fortran 95 code to implement the model and methods. Other computational platforms, such as Java, will also be explored.

Slater, John W.↗

Analysis of free turbulent shear flows by numerical methods

Studies are described in which the effort was essentially directed to classes of problems where the phenomenologically interpreted effective transport coefficients could be absorbed by, and subsequently extracted from (by comparison with experimental data), appropriate coordinate transformations. The transformed system of differential equations could then be solved without further specifications or assumptions by numerical integration procedures. An attempt was made to delineate different regimes for which specific eddy viscosity models could be formulated. In particular, this would account for the carryover of turbulence from attached boundary layers, the transitory adjustment, and the asymptotic behavior of initially disturbed mixing regions. Such models were subsequently used in seeking solutions for the prescribed two-dimensional test cases, yielding a better insight into overall aspects of the exchange mechanisms.

Korst, H. H.↗

Grid generation and its application to separated flows

A method for the computation of recirculating flows on a nonorthogonal boundary-fitted coordinate system is presented. The algebraic grid generation technique is used to generate grids in a flow field of arbitrary shape. Computations are performed by using the set of governing equations transformed for use on a generalized curvilinear coordinate system. Turbulent flow computations were made by using the k-epsilon model for most of the flow cases. Three pressure correction algorithms - SIMPLE, SIMPLEC, and PISO - were incorporated for use on a nonorthogonal grid. Comparisons were made on the basis of numerical stability and accuracy. The present method was found to be efficient for computations of recirculating flows.

Maruszewski, J. P.↗

F/A-18 1/9th scale model tail buffet measurements

Wind tunnel tests were carried out on a 1/9th scale model of the F/A-18 at high angles of attack to investigate the characteristics of tail buffet due to bursting of the wing leading edge extension (LEX) vortices. The tests were carried out at the Aeronautical Research Laboratory low-speed wind tunnel facility and form part of a collaborative activity with NASA Ames Research Center, organized by The Technical Cooperative Program (TTCP). Information from the program will be used in the planning of similar collaborative tests, to be carried out at NASA Ames, on a full-scale aircraft. The program covered the measurement of unsteady pressures and fin vibration for cases with and without the wing LEX fences fitted. Fourier transform methods were used to analyze the unsteady data, and information on the spatial and temporal content of the vortex burst pressure field was obtained. Flow visualization of the vortex behavior was carried out using smoke and a laser light sheet technique.

Martin, C. A.↗

Transonic swirling nozzle flow

A numerical model of viscous transonic swirling flow in axisymmetric nozzles is developed. MacCormack's implicit Gauss-Seidel method is applied to the thin-layer Navier-Stokes equations in transformed coordinates. Numerical results are compared with experimental data to validate the method. The effect of swirl and viscosity on nozzle performance are demonstrated by examining wall pressures, Mach contours, and integral parameters.

Keith, Theo G., Jr.↗

Calculation of Transformer Leakage Inductance by Simplified Flux Path Geometries

Series inductance is a critical energy storage element in isolated power converters. Many have explored methods of estimating the total leakage inductance of transformers as an integrated series inductance without external inductors. However, this poses challenges with medium frequency converters and magnetic ribbon-based cores, e.g., metal amorphous nanocrystalline. The material’s relatively high conductivity allows eddy currents to develop when magnetic flux intersects wide surfaces. These eddy currents result in significant additional losses that increase with increasing power flow. This paper presents a simplified and accurate method of estimating leakage flux while identifying separate leakage flux paths that correlate with the magnetic material surfaces. This enables an understanding the proportion of leakage flux that will contribute to extra losses. The methods presented are applicable to both leakage and fringing (around heterogeneous material interfaces) fluxes. This paper includes detailed FEA studies and 3D flux vector measurements matching the presented analytical models.

transformer↗

Sonic-box method employing local Mach number for oscillating wings with thickness

A computer program was developed to account approximately for the effects of finite wing thickness in the transonic potential flow over an oscillating wing of finite span. The program is based on the original sonic-box program for planar wing which was previously extended to include the effects of the swept trailing edge and the thickness of the wing. Account for the nonuniform flow caused by finite thickness is made by application of the local linearization concept. The thickness effect, expressed in terms of the local Mach number, is included in the basic solution to replace the coordinate transformation method used in the earlier work. Calculations were made for a delta wing and a rectangular wing performing plunge and pitch oscillations, and the results were compared with those obtained from other methods. An input quide and a complete listing of the computer code are presented.

Ruo, S. Y.↗

Analysis of eccentric annular incompressible seals. I - A new solution using fast Fourier transforms for determining hydrodynamic force

A new analysis procedure is presented which solves for the flow variables of an incompressible-flow annular pressure seal in which the rotor has a large static displacement from the centered position. The analysis begins with a set of governing equations based on a turbulent bulk-flow model and Moody's friction equation. No simplification of these bulk-flow equations is required for the solution procedure. Perturbation of the flow variables yields a set of zeroth and first-order equations. The zeroth-order equations (which model the large static displacement) are integrated by means of an efficient new method which employs fast Fourier transforms. Further integration of the zeroth-order pressures yields the hydrodynamic reactive force. Predictions for the hydrodynamic forces from this analysis procedure are in excellent agreement with available experimental results.

Nelson, C. C.↗

Calculation of turbulent boundary layers with heat transfer and pressure gradient utilizing a compressibility transformation. Part 2: Constant property turbulent boundary layer flow with simultaneous mass transfer and pressure gradient

An analysis of the incompressible turbulent boundary layer, developing under the combined effects of mass transfer and pressure gradient, is presented in this paper. A strip-integral method is employed whereby two of the three governing equations are obtained by integrating the combined momentum and continuity equation to 50 percent and 100 percent, respectively, of the boundary-layer height. The latter equation is the usual momentum-integral equation; the former equation requires specification of shear. Accordingly, Clauser's equilibrium eddy-viscosity law is assumed valid at this point. The third and final equation is obtained by specifying that Stevenson's velocity profiles apply throughout the domain of interest, from which a skin-friction law can be derived. Comparisons of the numerical results with the experiments of McQuaid, which include combined effects of variable pressure gradient and mass transfer, show good agreement.

Boccio, J.↗

Theoretical prediction of interference loading on aircraft stores. Part 2: Supersonic speeds

A method is developed for theoretically predicting the loading on pylon-mounted stores in supersonic flow. Linear theory is used, without two dimensional or slender body assumptions, to predict the flow field produced by the aircraft wing, nose, inlet, and pylons. Aircraft shock wave locations are predicted, and their effect on the flow field is included through a transformation of the aircraft geometry. The interference loading is integrated over the store length by considering the local crossflow, its axial and radial derivatives, and buoyancy. Store moment calculations under an F-4 aircraft at Mach 1.2 are compared to wind tunnel data. The method is computerized, and program user information is included.

Danfernandes, F.↗