SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1988

An explicit Runge-Kutta method for unsteady rotor/stator interaction

A quasi-three-dimensional rotor/stator analysis has been developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using a four-stage Runge-Kutta scheme with a constant timestep. Results are shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Euler and Navier-Stokes results are compared on the scaled single- and multi-passage machine. The method is relatively fast and the quasi-three-dimensional formulation is applicable to a wide range of turbomachinery geometries.

Jorgenson, Philip C. E.↗

NASA NTRS1988

An explicit Runge-Kutta method for unsteady rotor/stator interaction

Jorgenson, Philip C. E.↗

NASA NTRS1989

The 3-D Euler and Navier-Stokes calculations for aircraft components

An explicit multistage Runge-Kutta type of time-stepping scheme is used for solving transonic flow past a transport type wing/fuselage configuration. Solutions for both Euler and Navier-Stokes equations are obtained for quantitative assessment of boundary layer interaction effects. The viscous solutions are obtained on both a medium resolution grid of approximately 270,000 points and a find grid of 460,000 points to assess the effects of grid density on the solution. Computed pressure distributions are compared with the experimental data.

Vatsa, Veer N.↗

NASA NTRS1987

Improvements to the aircraft Euler method

An unstructured mesh of tetrahedra offers an attractive approach to mesh generation for complex three dimensional shapes. The combination of a method for constructing tetrahedral meshes and a finite element technique for solving the Euler equations has resulted in a powerful new approach to the problem of calculating flows over complex geometries. Recent progress and improvements to the authors' aircraft Euler method are described and results presented to demonstrate its capability.

Jameson, A.↗

NASA NTRS1993

Methodology for CFD Design Analysis of National Launch System Nozzle Manifold

The current design environment dictates that high technology CFD (Computational Fluid Dynamics) analysis produce quality results in a timely manner if it is to be integrated into the design process. The design methodology outlined describes the CFD analysis of an NLS (National Launch System) nozzle film cooling manifold. The objective of the analysis was to obtain a qualitative estimate for the flow distribution within the manifold. A complex, 3D, multiple zone, structured grid was generated from a 3D CAD file of the geometry. A Euler solution was computed with a fully implicit compressible flow solver. Post processing consisted of full 3D color graphics and mass averaged performance. The result was a qualitative CFD solution that provided the design team with relevant information concerning the flow distribution in and performance characteristics of the film cooling manifold within an effective time frame. Also, this design methodology was the foundation for a quick turnaround CFD analysis of the next iteration in the manifold design.

Haire, Scot L.↗

NASA NTRS1993

Calculation methods for rotor wake formation in forward flight

The primary focus of this research centered on improving the vorticity-capturing properties of Eulerian flow solvers in convection-dominated flows. The work was motivated by the fact that excessive numerical diffusion observed in current Euler solvers on grids of feasible size is a major obstacle to accurately computing the flow about helicopter rotors, especially in the forward flight case. The work was originally supported by NASA Grant NAG2-421 for two years. The research then continued for four more years under an ONR contract. The research goal of reduction of numerical diffusion in Euler solvers remained the same, although the target application was now the prediction of internal vortical secondary flows. This was part of a larger ONR-sponsored research program at MIT aimed at reducing the noise of turbopumps. The research work culminated in the PhD thesis titled 'A Coupled Eulerian/Lagrangian Method for the Solution of Three-Dimensional Vortical Flows'.

Drela, Mark↗

NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analyis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions.

Huff, Dennis L.↗

NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the perioidic boundary conditions. Approximate nonreflecting boundary conditons were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distribution for certain flow conditions.

Huff, Dennis L.↗

NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and non-zero inter-blade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions have been coded for the inlet and exit boundary conditions. Sample unsteady solutions have been performed for an oscillating cascade and compared to experimental data. Also, test cases were fun for a flat plate cascade to compare with an unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with non-zero inter-blade phase angles are in good agreement with experimental data and small-perturbation theory. The zero degree inter-blade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus non-reflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions. This code is expected to be used as a tool for reviewing simpler models that do not include the full nonlinear aerodynamics or as a final check for designs against flutter in turbomachinery.

Huff, Dennis L.↗

NASA NTRS1987

Euler calculations for wings using Cartesian grids

A method is presented for the calculation of transonic flows past wings using Cartesian grids. The calculations are based on a finite volume formulation of the Euler equations. Results are presented for a rectangular wing with a flat tip and the ONERA M6 wing. In general, the results are in good agreement with other computations and available experiment. However, Cartesian grids require a greater number of points than body fitted grids in order to resolve the flow properties near the leading edge of a swept wing.

Gaffney, R. L., Jr.↗

DOE OSTI2020

Lagrangian and Eulerian drag models that are consistent between Euler-Lagrange and Euler-Euler (two-fluid) approaches for homogeneous systems

The undisturbed flow of a particle is of fundamental importance since it controls both the undisturbed flow force and the perturbation force (which includes quasisteady, added-mass, and history forces). Here we use the pairwise interaction extended point particle framework to evaluate the undisturbed flow of each particle through superposition of the perturbation flow induced by all its neighbors. This approach allows calculation of various statistics related to undisturbed fluid velocity under conditions of both stationary and nonstationary particles. In a random distribution of stationary particles, while the macroscale undisturbed flow is slowly varying, the microscale undisturbed flow that arises due to the perturbation flow of neighbors varies substantially from one particle to another and this in turn leads to large variation in the hydrodynamic force exerted on the particles. The effect of particle motion is generally to increase the particle-to-particle variation in the undisturbed fluid velocity of the particles. We observe that this increase is greater for the transverse component than for the streamwise component. As a result, with increasing random particle motion, the distribution of undisturbed fluid-velocity fluctuation becomes isotropic. Three different normalized forces are defined: $Φ^L$ is the Lagrangian normalized force on an individual particle suitable for application in a microscale-informed Euler-Lagrange simulation, $Φ^E$ is the Eulerian normalized average force suitable for application in an Euler-Euler simulation, and $Φ^{LE}$ is the Lagrangian normalized force on an individual particle suitable for application in the standard Euler-Lagrange simulation. In this work, we establish precise relations between these different definitions. The drag laws developed based on particle-resolved direct numerical simulation results and experiments are appropriate for application only as the Eulerian normalized average force. We introduce the force consistency relation and use it to obtain an expression for $Φ^L$, which when applied to each particle and averaged over all the particles equals $Φ^E$. The results are first obtained in the limit of stationary particles and then extended to the general case of nonstationary particles.

42 ENGINEERING↗

NASA NTRS1986

Design of 9.271-pressure-ratio 5-stage core compressor and overall performance for first 3 stages

Overall aerodynamic design information is given for all five stages of an axial flow core compressor (74A) having a 9.271 pressure ratio and 29.710 kg/sec flow. For the inlet stage group (first three stages), detailed blade element design information and experimental overall performance are given. At rotor 1 inlet tip speed was 430.291 m/sec, and hub to tip radius ratio was 0.488. A low number of blades per row was achieved by the use of low-aspect-ratio blading of moderate solidity. The high reaction stages have about equal energy addition. Radial energy varied to give constant total pressure at the rotor exit. The blade element profile and shock losses and the incidence and deviation angles were based on relevant experimental data. Blade shapes are mostly double circular arc. Analysis by a three-dimensional Euler code verified the experimentally measured high flow at design speed and IGV-stator setting angles. An optimization code gave an optimal IGV-stator reset schedule for higher measured efficiency at all speeds.

Steinke, Ronald J.↗

NASA NTRS1989

Numerical solutions of forward-flight rotor flow using an upwind method

A finite-volume upwind algorithm for solving the 3-D Euler equations with a moving grid has been developed for computing helicopter forward-flight rotor flows. The computed pressure distributions and shock positions of high-speed rotor flow are compared with various experimental data as well as with other numerical results, and the agreement is encouraging. A comparison of quasi-steady solutions with unsteady solutions reveals that when a shock occurs in the flowfield, the assumption of quasi-steady flow may fail due to the time-lag of the shock motion. Similarly, three-dimensional effects cannot be neglected. Sufficient subiterations for each time step are required to avoid numerical lag effects in using the present method. The redistribution of the residual due to the coordinate transformation is discussed. For high-order MUSCL-type schemes, a coordinate-independent solution can be obtained by interpolating primitive variables.

Chen, C. L.↗

NASA NTRS1992

Prediction of vortex breakdown on a delta wing

Recent studies of leading-edge vortex flows with computational fluid dynamics codes using Euler or Navier-Stokes formulations have shown fair agreement with experimental data. These studies have concentrated on simulating the flowfields associated with a sharp-edged flat plate 70 deg delta wing at angles of attack where vortex breakdown or burst is observed over the wing. There are, however, a number of discrepancies between the experimental data and the computed flowfields. The location of vortex breakdown in the computational solutions is seen to differ from the experimental data and to vary with changes in the computational grid and freestream Mach number. There also remain issues as to the validity of steady-state computations for cases which contain regions of unsteady flow, such as in the post-breakdown regions. As a partial response to these questions, a number of laminar Navier-Stokes solutions were examined for the 70 deg delta wing. The computed solutions are compared with an experimental database obtained at low subsonic speeds. The convergence of forces, moments and vortex breakdown locations are also analyzed to determine if the computed flowfields actually reach steady-state conditions.

Agrawal, S.↗

NASA NTRS1993

Computational techniques for high-speed flows with viscous and chemical effects

Algorithms for solving the Euler and the Navier-Stokes equations in conjunction with chemical kinetic equations are presented. The convective flux is estimated from a quasi one dimensional interpolation procedure. Shock, contact, and expansion waves and thermochemical nonequilibrium phenomena are captured by the Lax-Friedrichs technique. Relaxation techniques were developed to enhance their effectiveness in dealing with spatial and temporal stiffness associated with the physical problems. Both explicit and implicit smoothers were implemented into the standard multigrid time stepping method. Unsteady and steady scalar problems are discussed. A perfect gas and equilibrium air shock tube problem is investigated. Numerical schemes and techniques are compared for the problems of shock and boundary layer interaction and three dimensional viscous, nonequilibrium flow encompassing an aerobrake. The results are comparable in accuracy against other high order non-oscillatory techniques. The multigrid methods are assessed using a Mach 8 flow over a complete planar body, a sphere, and a blunt delta wing at 20 deg incidence. Applying an implicit multigrid method on a nested grid of 128 by 64 nodes the reduction factor is 0.25. The central processing unit reduction factor is 2.2 after both the single and multigrid Runge-Kutta solutions converged to machine zero on a grid of 37 by 41 by 73 nodes.

Li, C. P.↗

NASA NTRS1992

Navier-Stokes analysis of an oxidizer turbine blade with tip clearance

The Gas Generator Oxidizer Turbine (GGOT) Blade is being analyzed by various investigators under the NASA MSFC sponsored Turbine Stage Technology Team design effort. The present work concentrates on the tip clearance region flow and associated losses; however, flow details for the passage region are also obtained in the simulations. The present calculations simulate the rotor blade row in a rotating reference frame with the appropriate coriolis and centrifugal acceleration terms included in the momentum equation. The upstream computational boundary is located about one axial chord from the blade leading edge. The boundary conditions at this location were determined by using a Euler analysis without the vanes to obtain approximately the same flow profiles at the rotor as were obtained with the Euler stage analysis including the vanes. Inflow boundary layer profiles are then constructed assuming the skin friction coefficient at both the hub and the casing. The downstream computational boundary is located about one axial chord from the blade trailing edge, and the circumferentially averaged static pressure at this location was also obtained from the Euler analysis. Results were obtained for the 3-D baseline GGOT geometry at the full scale design Reynolds number. Details of the clearance region flow behavior and blade pressure distributions were computed. The spanwise variation in blade loading distributions are shown, and circumferentially averaged spanwise distributions of total pressure, total temperature, Mach number, and flow angle are shown at several axial stations. The spanwise variation of relative total pressure loss shows a region of high loss in the region near the casing. Particle traces in the near tip region show vortical behavior of the fluid which passes through the clearance region and exits at the downstream edge of the gap.

Gibeling, Howard J.↗

NASA NTRS1994

Canonical-variables multigrid method for steady-state Euler equations

A novel approach for the solution of inviscid flow problems for subsonic compressible flows is described. The approach is based on canonical forms of the equations, in which subsystems governed by hyperbolic operators are separated from those governed by elliptic ones. The discretizations used as well as the iterative techniques for the different subsystems are inherently different. Hyperbolic parts, which describe, in general, propagation phenomena, are discretized using upwind schemes and are solved by marching techniques. Elliptic parts, which are directionally unbiased, are discretized using h-elliptic central discretizations, and are solved by pointwise relaxations together with coarse grid acceleration. The resulting discretization schemes introduce artificial viscosity only for the hyperbolic parts of the system; thus a smaller total artificial viscosity is used, while the multigrid solvers used are much more efficient. Solutions of the subsonic compressible Euler equations are achieved at the same efficiency as the full potential equation.

Taasan, Shlomo↗

NASA NTRS1993

An improved far field drag calculation method for nonlinear CFD codes

An improved method to calculate drag based on far field conditions is presented and demonstrated. The method is illustrated using two CFD codes, PARC and CFL3D, and examining their ability to preserve lift and drag coefficients along grid line defined integration paths. The flow fields were generated by solving the Euler equations for a NACA 0012 airfoil at a free stream Mach number of 0.8 and angle of attack of 1.25 degrees. In comparison to force coefficients obtained by surface pressure integration, neither code acceptably preserved both force coefficients throughout the near and far fields. An investigation into the relationship between numerical prediction error and calculated force coefficients revealed a direct connection between solution mass conservation error and force coefficients error. A method to correct the predicted force coefficients based on integrated mass conservation error is described and demonstrated. The corrected force coefficients for both codes are shown to be more accurate than the uncorrected values. The correction method is applicable to both two and three dimensions and is independent of the algorithm used to generate the flow field.

Sirbaugh, James R.↗