SEARCH · Search NASA

Results for “step 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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS2015

The Effect of Backward-Facing Step Height on Instability Growth and Breakdown in Swept Wing Boundary-Layer Transition

A low-speed experiment was performed on a swept at plate model with an imposed pressure gradient to determine the effect of a backward-facing step on transition in a stationary-cross flow dominated flow. Detailed hot-wire boundary-layer measurements were performed for three backward-facing step heights of approximately 36, 45, and 49% of the boundary-layer thickness at the step. These step heights correspond to a subcritical, nearly-critical, and critical case. Three leading-edge roughness configurations were tested to determine the effect of stationary-cross flow amplitude on transition. The step caused a local increase in amplitude of the stationary cross flow for the two larger step height cases, but farther downstream the amplitude decreased and remained below the baseline amplitude. The smallest step caused a slight local decrease in amplitude of the primary stationary cross flow mode, but the amplitude collapsed back to the baseline case far downstream of the step. The effect of the step on the amplitude of the primary cross flow mode increased with step height, however, the stationary cross flow amplitudes remained low and thus, stationary cross flow was not solely responsible for transition. Unsteady disturbances were present downstream of the step for all three step heights, and the amplitudes increased with increasing step height. The only exception is that the lower frequency (traveling crossflow-like) disturbance was not present in the lowest step height case. Positive and negative spikes in instantaneous velocity began to occur for the two larger step height cases and then grew in number and amplitude downstream of reattachment, eventually leading to transition. The number and amplitude of spikes varied depending on the step height and cross flow amplitude. Despite the low amplitude of the disturbances in the intermediate step height case, breakdown began to occur intermittently and the flow underwent a long transition region.

Eppink, Jenna L.↗

NASA NTRS1995

Advanced in Visualization of 3D Time-Dependent CFD Solutions

Numerical simulations of complex 3D time-dependent (unsteady) flows are becoming increasingly feasible because of the progress in computing systems. Unfortunately, many existing flow visualization systems were developed for time-independent (steady) solutions and do not adequately depict solutions from unsteady flow simulations. Furthermore, most systems only handle one time step of the solutions individually and do not consider the time-dependent nature of the solutions. For example, instantaneous streamlines are computed by tracking the particles using one time step of the solution. However, for streaklines and timelines, particles need to be tracked through all time steps. Streaklines can reveal quite different information about the flow than those revealed by instantaneous streamlines. Comparisons of instantaneous streamlines with dynamic streaklines are shown. For a complex 3D flow simulation, it is common to generate a grid system with several millions of grid points and to have tens of thousands of time steps. The disk requirement for storing the flow data can easily be tens of gigabytes. Visualizing solutions of this magnitude is a challenging problem with today's computer hardware technology. Even interactive visualization of one time step of the flow data can be a problem for some existing flow visualization systems because of the size of the grid. Current approaches for visualizing complex 3D time-dependent CFD solutions are described. The flow visualization system developed at NASA Ames Research Center to compute time-dependent particle traces from unsteady CFD solutions is described. The system computes particle traces (streaklines) by integrating through the time steps. This system has been used by several NASA scientists to visualize their CFD time-dependent solutions. The flow visualization capabilities of this system are described, and visualization results are shown.

Lane, David A.↗

NASA NTRS1974

Fracture control methods for space vehicles. Volume 1: Fracture control design methods

A systematic approach for applying methods for fracture control in the structural components of space vehicles consists of four major steps. The first step is to define the primary load-carrying structural elements and the type of load, environment, and design stress levels acting upon them. The second step is to identify the potential fracture-critical parts by means of a selection logic flow diagram. The third step is to evaluate the safe-life and fail-safe capabilities of the specified part. The last step in the sequence is to apply the control procedures that will prevent damage to the fracture-critical parts. The fracture control methods discussed include fatigue design and analysis methods, methods for preventing crack-like defects, fracture mechanics analysis methods, and nondestructive evaluation methods. An example problem is presented for evaluation of the safe-crack-growth capability of the space shuttle crew compartment skin structure.

Liu, A. F.↗

NASA NTRS1976

Viscid-inviscid interaction of incompressible separated flows

A flow model has been devised to deal with the viscid-inviscid interaction of a class of two-dimensional incompressible separated flow problems. It is suggested that the corresponding inviscid flow of these problems is described by the free streamline theory with few unspecified parameters and their values are, in turn, determined by the viscous flow considerations. The problem of a flow past a backward facing step is selected for study in detail. The viscous flow components of turbulent jet mixing, recompression, and reattachment are delineated and studied individually. When they are later combined, it is found that the point of reattachment behaves as a saddle-point-type singularity in the system of differential equations describing the viscous flow process. This feature is employed to the determination of the aforementioned free parameters and thus the establishment of the overall corresponding inviscid flow field. The resulting base pressure coefficient for the specific case agrees reasonably well with the available experimental data.

Chow, W. L.↗

NASA NTRS1986

Comparison of Euler and full potential marching techniques for flows over complex configurations

Two recently developed aerodynamic prediction techniques based on the steady full potential equation and the unsteady Euler equations have been applied to a variety of three-dimensional supersonic flow problems exhibiting embedded subsonic regions. Both techniques utilize planar Gauss-Seidel relaxation in the marching direction and approximate factorization in the cross-flow plane. A conservative switching scheme and flux bias technique are employed in the full potential method to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. A new unified approach with finite volume, high accuracy (up to third order) Total Variation Diminishing formulation (based on Roe's scheme) is used in the Euler solver. In the supersonic regions of the flow an 'infinitely large' time step is employed, and a finite time step is applied in the subsonic regions of the flow to reach the steady-state as a time-asymptote. Numerical solutions are obtained for a number of complex configurations, including: (1) an elliptic waverider, (2) a realistic fighter configuration, (3) the Space Shuttle, and (4) a Shuttle-like configuration. Both the Full Potential and Euler numerical results are in good agreement with available experimental data.

Szema, K. Y.↗

NASA NTRS1992

An experimental investigation of the separating/reattaching flow over a backstep

Transport characteristics of the turbulent kinetic energy, k, and the shear stream, bar(-uv), were studied in an incompressible boundary layer downstream of the reattachment of the separated flow behind a backward facing step. Hot-wire measurement techniques were used to measure three Reynolds stresses and higher order mean products of velocity fluctuations. These quantities were used to evaluate advection, turbulent diffusion, and production terms of the k and bar(-uv) transport equations. The dissipation rate and pressure-strain terms were obtained as the difference of all the other terms. The Reynolds number based on the step height was 37,000 and the upstream oncoming flow was fully developed turbulent boundary layer with the R(sub theta) = 3600.

Jovic, Srboljub↗

NASA NTRS1995

Fuel Injector With Shear Atomizer

Atomizer for injecting liquid fuel into combustion chamber uses impact and swirl to break incoming stream of fuel into small, more combustible droplets. Slanted holes direct flow of liquid fuel to stepped cylindrical wall. Impact on wall atomizes liquid. Air flowing past vanes entrains droplets of liquid in swirling flow. Fuel injected at pressure lower than customarily needed.

Beal, George W.↗

NASA NTRS1994

Turbulent Heat Transfer Characteristics in the Shear Layer of a Separated Flow

Experiments were performed to study the evolution of the heat transfer structure in a separated free shear layer region of an incompressible separated turbulent boundary layer flow behind a backward-facing step. While there is an abundance of velocity field measurements of separated flows, heat transfer measurements are rather scarce, thus limiting assessment of the heat transfer physics and its accurate modeling. The purpose of the paper is twofold: to improve an understanding of effects of flow separation on heat transfer characteristics, and to provide data for turbulence modeling and computation. The boundary layer upstream of the step was turbulent and fully developed. A constant temperature surface boundary condition was imposed upstream and downstream of the step for the heat transfer study. An internal mixing-layer like flow forms and grows from the step lip within the original boundary layer. The turbulent structure of the flow evolving downstream, however, does not switch immediately to that of a mixing layer over the entire shear layer thickness. Measurements of mean and fluctuating velocity and temperature fields indicate that the internal layer spreads gradually in the transverse direction while the outer part of the original boundary layer is effectively unperturbed. The results in this paper have not been previously reported.

Jovic, S.↗

NASA NTRSDate not supplied

Numerical Investigation of Combustor-Turbine Interactions with a Two-Stage High-Pressure Turbine Including Cooling Airflows

Combustor–turbine interactions (CTIs) are investigated by performing three-dimensional unsteady simulations using a realistic combustor and high-pressure turbine (HPT) configurations from the Energy Efficient Engine (E3) program. To understand CTIs, we compare the numerically predicted flow fields from single-component simulations (Step 1: the combustor + the first-stage stator of turbine; Step 2: the two-stage HPT imposing the time-averaged flow solution from Step 1 as the inflow boundary condition) and a fully coupled combustor–turbine simulation (Step 3) at the sea-level takeoff (SLTO) condition. In addition to three previous simulations where the cooling airflows inside the HPT had been neglected for all Step 1, Step 2, and Step 3, two new simulations of Step 2 and Step 3 take into account the cooling airflows using the source team approach. In this approach, to mimic the cooling airflows, we impose the source term at a specific area of a cooling airflow hole at the solid surfaces without making a mesh of each hole. The objective of this study is twofold. One is to investigate the effect of the presence of the cooling airflows on the aerodynamics of the combustor and HPT as well as the HPT performance. The second is to perform a detailed comparison among the calculated flowfields by two different numerical schemes, the central-difference with the standard Jameson–Schmidt–Turkel (CD-JST) scheme and the AUSM scheme. There is a noticeable difference in the hot-streak distributions at the first-stage stator. In addition, depending on the choice of the numerical scheme and the presence of the cooling airflows, an occurrence of shock waves appearing at the first-stage stator is greatly influenced. Thus, this has a noticeable impact on the HPT performance. It is shown that the estimated turbine efficiencies from Step 3 are about 7 % less than the ones from Step 2.

Combustor-turbine interaction↗

NASA NTRS1988

Body conforming grids for general unsteady airfoil motion

A method for computing a new body conforming grid at each time step of an unsteady flow calculation is presented. The method is based on a spring system analogy. It simply updates the grid from the previous time step, hence it requires that an initial grid be created using another method. The method is designed to be independent of the flow equations used, so it can be easily transferred between aerodynamic codes with the same type of initial grid (i.e., O-type). The initial development was done using an O-type grid. It was then extended to C-type grids which required only minor changes to the routine. Results are presented for both O-type and C-type grids. The quality of the updated grids were examined by observing their movement and the flow results they produced. The effect due to the shape of the prescribed grid lines for the O-type grids was examined for a NACA 0012 airfoil undergoing harmonic pitch. The time history of the pressure coefficients of this case was compared with the interpolation scheme and experimental data.

Levine, M.↗

NASA NTRS1994

Preconditioned conjugate gradient methods for the Navier-Stokes equations

A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulation. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a conventional implicit solver, namely line Gauss-Seidel relaxation (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU time on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow case are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner.

Ajmani, Kumud↗

NASA NTRS1995

Some Developments of the Equilibrium Particle Simulation Method for the Direct Simulation of Compressible Flows

The direct simulation Monte Carlo (DSMC) method is the established technique for the simulation of rarefied gas flows. In some flows of engineering interest, such as occur for aero-braking spacecraft in the upper atmosphere, DSMC can become prohibitively expensive in CPU time because some regions of the flow, particularly on the windward side of blunt bodies, become collision dominated. As an alternative to using a hybrid DSMC and continuum gas solver (Euler or Navier-Stokes solver) this work is aimed at making the particle simulation method efficient in the high density regions of the flow. A high density, infinite collision rate limit of DSMC, the Equilibrium Particle Simulation method (EPSM) was proposed some 15 years ago. EPSM is developed here for the flow of a gas consisting of many different species of molecules and is shown to be computationally efficient (compared to DSMC) for high collision rate flows. It thus offers great potential as part of a hybrid DSMC/EPSM code which could handle flows in the transition regime between rarefied gas flows and fully continuum flows. As a first step towards this goal a pure EPSM code is described. The next step of combining DSMC and EPSM is not attempted here but should be straightforward. EPSM and DSMC are applied to Taylor-Couette flow with Kn = 0.02 and 0.0133 and S(omega) = 3). Toroidal vortices develop for both methods but some differences are found, as might be expected for the given flow conditions. EPSM appears to be less sensitive to the sequence of random numbers used in the simulation than is DSMC and may also be more dissipative. The question of the origin and the magnitude of the dissipation in EPSM is addressed. It is suggested that this analysis is also relevant to DSMC when the usual accuracy requirements on the cell size and decoupling time step are relaxed in the interests of computational efficiency.

Macrossan, M. N.↗

NASA NTRS1989

The use of a Navier-Stokes code in the wing design process

The feasibility was determined of incorporating the Navier-Stokes computational code, CFL3D, into the supersonic wing design process. The approach taken is of two steps. The first step was to calibrate CFL3D against existing experimental data sets obtained on thin sharp edged delta wings. The experimental data identified six flow types which are dependent on the similarity parameters of Mach number and angle of attack normal to the leading edge. The calibration showed CFL3D capable of simulating these various separated and attached flow conditions. The second step was to use CFL3D to study the initial formation of leading edge separation over delta wings at supersonic speeds. This consisted of examining solutions obtained on a 65 deg delta wing at Mach number of 1.6 with varying cross sectional shapes. Reynolds number was held constant at 1000000 and the Baldwin-Lomax turbulence model was used. The study showed that through the use of leading edge radius and/or camber, the onset of leading edge separation can be delayed to a higher angle of attack than observed on a flat sharp edged wing. Based on the geometries studied, three wind tunnel models are being designed to verify these results.

Mcmillin, S. Naomi↗

NASA NTRS2003

Large-Eddy Simulation Code Developed for Propulsion Applications

A large-eddy simulation (LES) code was developed at the NASA Glenn Research Center to provide more accurate and detailed computational analyses of propulsion flow fields. The accuracy of current computational fluid dynamics (CFD) methods is limited primarily by their inability to properly account for the turbulent motion present in virtually all propulsion flows. Because the efficiency and performance of a propulsion system are highly dependent on the details of this turbulent motion, it is critical for CFD to accurately model it. The LES code promises to give new CFD simulations an advantage over older methods by directly computing the large turbulent eddies, to correctly predict their effect on a propulsion system. Turbulent motion is a random, unsteady process whose behavior is difficult to predict through computer simulations. Current methods are based on Reynolds-Averaged Navier- Stokes (RANS) analyses that rely on models to represent the effect of turbulence within a flow field. The quality of the results depends on the quality of the model and its applicability to the type of flow field being studied. LES promises to be more accurate because it drastically reduces the amount of modeling necessary. It is the logical step toward improving turbulent flow predictions. In LES, the large-scale dominant turbulent motion is computed directly, leaving only the less significant small turbulent scales to be modeled. As part of the prediction, the LES method generates detailed information on the turbulence itself, providing important information for other applications, such as aeroacoustics. The LES code developed at Glenn for propulsion flow fields is being used to both analyze propulsion system components and test improved LES algorithms (subgrid-scale models, filters, and numerical schemes). The code solves the compressible Favre-filtered Navier- Stokes equations using an explicit fourth-order accurate numerical scheme, it incorporates a compressible form of Smagorinsky s model for the subgrid-scale turbulence, and it uses generalized curvilinear coordinates to allow analysis of a wide range of geometries. The code runs in parallel on shared memory multiprocessor computers and is written in Fortran 90 with dynamic memory allocation. A sample result for a Mach-1.4 round jet is presented in the figure. Instantaneous Mach number contours in several cross-planes downstream of the nozzle exit are shown, illustrating how an LES captures the large unsteady three-dimensional turbulent structures present in the jet.

DeBonis, James R.↗

NASA NTRS2009

Re-evaluation of an Optimized Second Order Backward Difference (BDF2OPT) Scheme for Unsteady Flow Applications

Recent experience in the application of an optimized, second-order, backward-difference (BDF2OPT) temporal scheme is reported. The primary focus of the work is on obtaining accurate solutions of the unsteady Reynolds-averaged Navier-Stokes equations over long periods of time for aerodynamic problems of interest. The baseline flow solver under consideration uses a particular BDF2OPT temporal scheme with a dual-time-stepping algorithm for advancing the flow solutions in time. Numerical difficulties are encountered with this scheme when the flow code is run for a large number of time steps, a behavior not seen with the standard second-order, backward-difference, temporal scheme. Based on a stability analysis, slight modifications to the BDF2OPT scheme are suggested. The performance and accuracy of this modified scheme is assessed by comparing the computational results with other numerical schemes and experimental data.

Vatsa, Veer N.↗

NASA NTRS1990

A Detailed Investigation of Staged Normal Injection into a Mach 2 Flow

A study of the staged injection of two jets of air behind a rearward facing step into a Mach 2 flow was performed using the SPARK 3-D Navier-Stokes code. Calculated mole fraction distributions were compared with an extensive set of planar mole fraction measurements made with a laser induced iodine fluorescence technique. A statistical measure, the standard deviation, was used to help assess agreement between calculation and experiment. Overall, good agreement was found between calculated and measured values. Generally, agreement was better in the far field of the injectors. The effect of grid resolution was investigated by calculating solutions on grids of 60,000, 200,000, and 450,000 points. Differences in the solutions on the two finer grids were small. However, the mole fraction distributions were distinguishable. The effect of turbulence modeling was investigated by employing three different algebraic models for the jet turbulence: the Baldwin-Lomax model, the Prandtl mixing length model, and the Eggers mixing length model. Overall, the Eggers mixing length model was found to be superior for this case. Finally, the effect of the jet exit conditions was examined. A recently proposed Mach number distribution at the jet exit was found to slightly improve agreement between measurement and calculation.

Eklund, Dean R.↗

NASA NTRS1980

Types of leeside flow over delta wings

It is noted that so far most systematic investigations on the lee side flow over delta wings at supersonic speeds are concerned with flat upper surfaces. On the basis of these results, the paper makes an attempt to characterize the different types of flow over a wing with a delta-shaped upper surface by varying a number of parameters. It is concluded that the work should be considered a first step toward systematizing the flow over delta-shaped lee sides as well.

Szodruch, J.↗

NASA NTRS1989

Solution techniques for incompressible flow problems

A three-step Petrov-Galerkin (PG)/operator spliting scheme for the time-dependent incompressible Navier-Stokes equations is proposed. Each time step is split into two Stokes problems and one nonlinear convection-diffusion problem. Using a PG technique on the two outer Stokes problems ensures a stable scheme despite equal-order interpolation, while using a streamline upwind PG scheme on the inner convection-diffusion problem ensures a numerically stable solution at high Reynolds numbers. Numerical tests of this method have been carried out.

Tezduyar, T. E.↗