SEARCH · Search NASA

Results for “Navier Stokes solver”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Use of Navier-Stokes methods for the calculation of high-speed nozzle flow fields

Flows through three reference nozzles have been calculated to determine the capabilities and limitations of the widely used Navier-Stokes solver, PARC. The nozzles examined have similar dominant flow characteristics as those considered for supersonic transport programs. Flows from an inverted velocity profile (IVP) nozzle, an under expanded nozzle, and an ejector nozzle were examined. PARC calculations were obtained with its standard algebraic turbulence model, Thomas, and the two-equation turbulence model, Chien k-epsilon. The Thomas model was run with the default coefficient of mixing set at both 0.09 and a larger value of 0.13 to improve the mixing prediction. Calculations using the default value substantially underpredicted the mixing for all three flows. The calculations obtained with the higher mixing coefficient better predicted mixing in the IVP and underexpanded nozzle flows but adversely affected PARC's convergence characteristics for the IVP nozzle case. The ejector nozzle case did not converge with the Thomas model and the higher mixing coefficient. The Chien k-epsilon results were in better agreement with the experimental data overall than were those of the Thomas run with the default mixing coefficient, but the default boundary conditions for k and epsilon underestimated the levels of mixing near the nozzle exits.

Georgiadis, Nicholas J.↗

NASA NTRS1994

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗

NASA NTRS1993

Numerical investigation of multi-element airfoils

The flow over multi-element airfoils with flat-plate lift-enhancing tabs was numerically investigated. Tabs ranging in height from 0.25 percent to 1.25 percent of the reference airfoil chord were studied near the trailing edge of the main-element. This two-dimensional numerical simulation employed an incompressible Navier-Stokes solver on a structured, embedded grid topology. New grid refinements were used to improve the accuracy of the solution near the overlapping grid boundaries. The effects of various tabs were studied at a constant Reynolds number on a two-element airfoil with a slotted flap. Both computed and measured results indicated that a tab in the main-element cove improved the maximum lift and lift-to-drag ratio relative to the baseline airfoil without a tab. Computed streamlines revealed that the additional turning caused by the tab may reduce the amount of separated flow on the flap. A three-element airfoil was also studied over a range of Reynolds numbers. For the optimized flap rigging, the computed and measured Reynolds number effects were similar. When the flap was moved from the optimum position, numerical results indicated that a tab may help to reoptimize the airfoil to within 1 percent of the optimum flap case.

Cummings, Russell M.↗

NASA NTRS1994

Hybrid structured/unstructured grid computations for the F/A-18 at high angle of attack

At high angles of attack, vortical flows play a crucial role in the maintenance of lift for fighter aircraft. However, under certain conditions, the vortical flow can have an adverse effect on the aircraft. On the F/A-18 the LEX vortex can impinge on the tail; at high angles of attack the unsteady flow from vortex bursting can cause structural fatigue on the vertical tails. At high angles of attack, the flow field can be quite complex, and a computational analysis challenging. A full configuration analysis with viscous structured grids can be computationally expensive and the task of generating the requisite grids can be quite difficult. To mitigate these difficulties and provide a medium-fidelity analysis tool, a hybrid structured/unstructured approach was adopted for this work. In this analysis, the formation and roll-up of the LEX vortex is computed with a structured Navier-Stokes solver, and the resulting vortex propagated downstream with an unstructured Euler solver.

Robert T Biedron↗

NASA NTRS1994

New Developments in the Method of Space-Time Conservation Element and Solution Element-Applications to Two-Dimensional Time-Marching Problems

A new numerical discretization method for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods, i.e., finite difference, finite volume, finite element, and spectral methods. It is motivated by several important physical/numerical considerations and designed to avoid several key limitations of the above traditional methods. As a result of the above considerations, a set of key principles for the design of numerical schemes was put forth in a previous report. These principles were used to construct several numerical schemes that model a 1-D time-dependent convection-diffusion equation. These schemes were then extended to solve the time-dependent Euler and Navier-Stokes equations of a perfect gas. It was shown that the above schemes compared favorably with the traditional schemes in simplicity, generality, and accuracy. In this report, the 2-D versions of the above schemes, except the Navier-Stokes solver, are constructed using the same set of design principles. Their constructions are simplified greatly by the use of a nontraditional space-time mesh. Its use results in the simplest stencil possible, i.e., a tetrahedron in a 3-D space-time with a vertex at the upper time level and other three at the lower time level. Because of the similarity in their design, each of the present 2-D solvers virtually shares with its 1-D counterpart the same fundamental characteristics. Moreover, it is shown that the present Euler solver is capable of generating highly accurate solutions for a famous 2-D shock reflection problem. Specifically, both the incident and the reflected shocks can be resolved by a single data point without the presence of numerical oscillations near the discontinuity.

Chang, Sin-Chung↗

NASA NTRS1994

Transonic Navier-Stokes calculations about a 65 deg delta wing

A computational study has been conducted in which the CFL3D Navier-Stokes solver coupled with an algebraic and a one-equation nonequilibrium turbulence model has been used to predict the flow over a 65 degree delta wing at transonic conditions for Reynolds numbers ranging from 6 x 10(exp 6) to 120 x 10(exp 6) based on mean aerodynamic chord. Solutions obtained indicated that the computational method when used with the one-equation turbulence model predicts results that compare well with experiment for attached flow conditions. Comparisons with experimental pressure at separated conditions show that the computational method, even though primary flow-field features are predicted well, does not predict secondary flow features.

Londenberg, W. Kelly↗

NASA NTRS1995

Application of Navier-Stokes code PAB3D with kappa-epsilon turbulence model to attached and separated flows

A three-dimensional Navier-Stokes solver was used to determine how accurately computations can predict local and average skin friction coefficients for attached and separated flows for simple experimental geometries. Algebraic and transport equation closures were used to model turbulence. To simulate anisotropic turbulence, the standard two-equation turbulence model was modified by adding nonlinear terms. The effects of both grid density and the turbulence model on the computed flow fields were also investigated and compared with available experimental data for subsonic and supersonic free-stream conditions.

Abdol-Hamid, Khaled S.↗

NASA NTRS1993

Progress in Advanced Spray Combustion Code Integration

A multiyear project to assemble a robust, muitiphase spray combustion code is now underway and gradually building up to full speed. The overall effort involves several university and government research teams as well as Rocketdyne. The first part of this paper will give an overview of the respective roles of the different participants involved, the master strategy, the evolutionary milestones, and an assessment of the state-of-the-art of various key components. The second half of this paper will highlight the progress made to date in extending the baseline Navier-Stokes solver to handle multiphase, multispecies, chemically reactive sub- to supersonic flows. The major hurdles to overcome in order to achieve significant speed ups are delineated and the approaches to overcoming them will be discussed.

Liang, Pak-Yan↗

NASA NTRS1995

Overview of turbulence model development and applications at Rocketdyne

This viewgraph presentation discusses turbulence modeling requirements, development philosophy, and approach; two major areas of concentration (high speed and low speed turbulence modeling); high speed turbulence modeling; compressibility effects; turbulence models adapted to USA code; M = 9.2 flat plate flow; Mach 7.05 flow over axisymmetric flare; Mach 8.6 flow over cold wall edge; low speed turbulence modeling; turbulence models being assessed; turbulence model deck structure and integration with Navier-Stokes solver; nonlinear algebraic-stress model; rotation modified k-epsilon model; and Reynolds stress model.

Hadid, A. H.↗

NASA NTRS1992

Evaluation of a doubly-swept blade tip for rotorcraft noise reduction

A computational study was performed for a doubly-swept rotor blade tip to determine its benefit for high-speed impulsive (HSI) and blade-vortex interaction (BVI) noise. This design consists of aft and forward sweep. For the HSI-noise computations, unsteady Euler calculations were performed for several variations to a rotor blade geometry. A doubly-swept planform was predicted to increase the delocalizing Mach number to 0.94 (representative of a 200+ kt helicopter). For the BVI-noise problem, it had been hypothesized that the doubly-swept blade tip, by producing a leading-edge vortex, would reduce the tip-vortex effect on BVI noise. A procedure was used in which the tip vortex velocity profile computed by a Navier-Stokes solver was used to compute the inflow associated with BVI. This inflow was used by a Euler solver to compute the unsteady pressures for an acoustic analysis. The results of this study were inconclusive due to the difficulty in accurately predicting the viscous tip vortex downstream of the blade. Also, for the condition studied, no leading-edge vortex formed at the tip.

Wake, Brian E.↗

NASA NTRS1995

Study of potential aerodynamic benefits from spanwise blowing at wingtip

Comprehensive experimental and analytical studies have been conducted to assess the potential aerodynamic benefits from spanwise blowing at the tip of a moderate-aspect-ratio swept wing. Previous studies on low-aspect-ratio wings indicated that blowing from the wingtip can diffuse the tip vortex and displace it outward. The diffused and displaced vortex will induce a smaller downwash at the wing, and consequently the wing will have increased lift and decreased induced drag at a given angle of attack. Results from the present investigation indicated that blowing from jets with a short chord had little effect on lift or drag, but blowing from jets with a longer chord increased lift near the tip and reduced drag at low Mach numbers. A Navier-Stokes solver with modified boundary conditions at the tip was used to extrapolate the results to a Mach number of 0.72. Calculations indicated that lift and drag increase with increasing jet momentum coefficient. Because the momentum of the jet is typically greater than the reduction in the wing drag and the increase in the wing lift due to spanwise blowing is small, spanwise blowing at the wingtip does not appear to be a practical means of improving the aerodynamic efficiency of moderate-aspectratio swept wings at high subsonic Mach numbers.

Mineck, Raymond E.↗

NASA NTRS1995

Use of the PARC code to estimate the off-design transonic performance of an over/under turboramjet nozzle

The transonic performance of a dual-throat, single-expansion-ramp nozzle (SERN) was investigated with a PARC computational fluid dynamics (CFD) code, an external flow Navier-Stokes solver. The nozzle configuration was from a conceptual Mach 5 cruise aircraft powered by four air-breathing turboramjets. Initial test cases used the two-dimensional version of PARC in Euler mode to investigate the effect of geometric variation on transonic performance. Additional cases used the two-dimensional version in viscous mode and the three-dimensional version in both Euler and viscous modes. Results of the analysis indicate low nozzle performance and a highly three-dimensional nozzle flow at transonic conditions. In another comparative study using the PARC code, a single-throat SERN configuration for which experimental data were available at transonic conditions was used to validate the results of the over/under turboramjet nozzle.

Lam, David W.↗

NASA NTRS1994

A comparison of turbulence models in computing multi-element airfoil flows

Four different turbulence models are used to compute the flow over a three-element airfoil configuration. These models are the one-equation Baldwin-Barth model, the one-equation Spalart-Allmaras model, a two-equation k-omega model, and a new one-equation Durbin-Mansour model. The flow is computed using the INS2D two-dimensional incompressible Navier-Stokes solver. An overset Chimera grid approach is utilized. Grid resolution tests are presented, and manual solution-adaptation of the grid was performed. The performance of each of the models is evaluated for test cases involving different angles-of-attack, Reynolds numbers, and flap riggings. The resulting surface pressure coefficients, skin friction, velocity profiles, and lift, drag, and moment coefficients are compared with experimental data. The models produce very similar results in most cases. Excellent agreement between computational and experimental surface pressures was observed, but only moderately good agreement was seen in the velocity profile data. In general, the difference between the predictions of the different models was less than the difference between the computational and experimental data.

Rogers, Stuart E.↗

NASA NTRS1993

Optimization of contoured hypersonic scramjet inlets with a least-squares parabolized Navier-Stokes procedure

A new optimization procedure, in which a parabolized Navier-Stokes solver is coupled with a non-linear least-squares optimization algorithm, is applied to the design of a Mach 14, laminar two-dimensional hypersonic subscale flight inlet with an internal contraction ratio of 15:1 and a length-to-throat half-height ratio of 150:1. An automated numerical search of multiple geometric wall contours, which are defined by polynomical splines, results in an optimal geometry that yields the maximum total-pressure recovery for the compression process. Optimal inlet geometry is obtained for both inviscid and viscous flows, with the assumption that the gas is either calorically or thermally perfect. The analysis with a calorically perfect gas results in an optimized inviscid inlet design that is defined by two cubic splines and yields a mass-weighted total-pressure recovery of 0.787, which is a 23% improvement compared with the optimized shock-canceled two-ramp inlet design. Similarly, the design procedure obtains the optimized contour for a viscous calorically perfect gas to yield a mass-weighted total-pressure recovery value of 0.749. Additionally, an optimized contour for a viscous thermally perfect gas is obtained to yield a mass-weighted total-pressure recovery value of 0.768. The design methodology incorporates both complex fluid dynamic physics and optimal search techniques without an excessive compromise of computational speed; hence, this methodology is a practical technique that is applicable to optimal inlet design procedures.

Korte, J. J.↗

NASA NTRS1994

Grid adaptation using chimera composite overlapping meshes

The objective of this paper is to perform grid adaptation using composite overlapping meshes in regions of large gradient to accurately capture the salient features during computation. The chimera grid scheme, a multiple overset mesh technique, is used in combination with a Navier-Stokes solver. The numerical solution is first converged to a steady state based on an initial coarse mesh. Solution-adaptive enhancement is then performed by using a secondary fine grid system which oversets on top of the base grid in the high-gradient region, but without requiring the mesh boundaries to join in any special way. Communications through boundary interfaces between those separated grids are carried out using trilinear interpolation. Application to the Euler equations for shock reflections and to shock wave/boundary layer interaction problem are tested. With the present method, the salient features are well-resolved.

Kao, Kai-Hsiung↗

NASA NTRS1995

Navier-Stokes calculations of rotor-airframe interaction in forward flight

A method for analyzing the mutual aerodynamic interaction between a rotor and an airframe model has been developed. This technique models the rotor implicitly through the source terms of the momentum equations. A three-dimensional, incompressible, laminar, Navier-Stokes solver in cylindrical coordinates was developed for analyzing the rotor-airframe problem. The calculations are performed on a simplified rotor-airframe model at an advance ratio of 0.1. The airframe surface pressure predictions are found to be in good agreement with wind tunnel test data. Results are also presented for velocity and pressure field distributions in the wake of the rotor.

Zori, Laith A. J.↗

NASA NTRS1995

Semi-span model testing in the National Transonic Facility

The objectives of the present research are: (1) to develop a computational approach to support semi-span model test techniques in the NTF; and (2) to integrate this approach with the conduct of an experimental test program. To meet these objectives, the following approach is taken. A state-of-the-art three-dimensional Navier-Stokes solver is employed to compute the flow over both a full-span configuration and a semi-span configuration mounted on the sidewall of the tunnel. The computations are validated by making direct comparisons to experimental data for both configurations. Then, the semi-span computational results are compared to the full-span results to document how the flow over the semi-span configuration differs from that over the full-span configuration. The results of this comparative study will be used to provide a conceptual framework within which a semi-spa model test technique may be implemented in the NTF.

Chokani, Ndaona↗

NASA NTRS1995

Runtime volume visualization for parallel CFD

This paper discusses some aspects of design of a data distributed, massively parallel volume rendering library for runtime visualization of parallel computational fluid dynamics simulations in a message-passing environment. Unlike the traditional scheme in which visualization is a postprocessing step, the rendering is done in place on each node processor. Computational scientists who run large-scale simulations on a massively parallel computer can thus perform interactive monitoring of their simulations. The current library provides an interface to handle volume data on rectilinear grids. The same design principles can be generalized to handle other types of grids. For demonstration, we run a parallel Navier-Stokes solver making use of this rendering library on the Intel Paragon XP/S. The interactive visual response achieved is found to be very useful. Performance studies show that the parallel rendering process is scalable with the size of the simulation as well as with the parallel computer.

Ma, Kwan-Liu↗