Search NASA⌕ Search

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.

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

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

X-33 Computational Aeroheating/Aerodynamic Predictions and Comparisons With Experimental Data

This report details a computational fluid dynamics study conducted in support of the phase II development of the X-33 vehicle. Aerodynamic and aeroheating predictions were generated for the X-33 vehicle at both flight and wind-tunnel test conditions using two finite-volume, Navier-Stokes solvers. Aerodynamic computations were performed at Mach 6 and Mach 10 wind-tunnel conditions for angles of attack from 10 to 50 with body-flap deflections of 0 to 20. Additional aerodynamic computations were performed over a parametric range of free-stream conditions at Mach numbers of 4 to 10 and angles of attack from 10 to 50. Laminar and turbulent wind-tunnel aeroheating computations were performed at Mach 6 for angles of attack of 20 to 40 with body-flap deflections of 0 to 20. Aeroheating computations were performed at four flight conditions with Mach numbers of 6.6 to 8.9 and angles of attack of 10 to 40. Surface heating and pressure distributions, surface streamlines, flow field information, and aerodynamic coefficients from these computations are presented, and comparisons are made with wind-tunnel data.

Hollis, Brian R.↗
NASA NTRS2002

Application of the Loci-Based CFD Code Chem at MSFC: Preliminary Results

Contents include the following: 1. Objectives. Concentrate on determining the qualitative accuracy, performance and robustness of the Chen code. 2. What is the Loci-Chem CFD code? Density-based, finite volume, generalized unstructured grid, Navier-Stokes solver. The algorithm was implemented using the Loci framework, which allows implementation issues such as parallel processing to be handled transparently to the coding of the CFD algorithm. 3. Application to Bifurcating Duct problem. Flow splits from single duct to two ducts. 4. Application to single element injector. 5. Application to PSU RBCC rig. 6. 90 degree elbow benchmark problem. 7. Future work.

West, Jeff S.↗
NASA NTRS2003

Nonlinear Aeroacoustics Computations by the Space-Time CE/SE Method

The Space-Time Conservation Element and Solution Element Method, or CE/SE Method for short, is a recently developed numerical method for conservation laws. Despite its second order accuracy in space and time, it possesses low dispersion errors and low dissipation. The method is robust enough to cover a wide range of compressible flows: from weak linear acoustic waves to strong discontinuous waves (shocks). An outstanding feature of the CE/SE scheme is its truly multi-dimensional, simple but effective non-reflecting boundary condition (NRBC), which is particularly valuable for computational aeroacoustics (CAA). In nature, the method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its careful treatment of the surface fluxes and geometry, it is different from the existing schemes. Currently, the CE/SE scheme has been developed to a matured stage that a 3-D unstructured CE/SE Navier-Stokes solver is already available. However, in the present review paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen and sketched in section 2. Then applications of the 2-D and 3-D CE/SE schemes to linear, and in particular, nonlinear aeroacoustics are depicted in sections 3, 4, and 5 to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS2003

A General Simulation Method for Multiple Bodies in Proximate Flight

Methods of unsteady aerodynamic simulation for an arbitrary number of independent bodies flying in close proximity are considered. A novel method to efficiently detect collision contact points is described. A method to compute body trajectories in response to aerodynamic loads, applied loads, and inter-body collisions is also given. The physical correctness of the methods are verified by comparison to a set of analytic solutions. The methods, combined with a Navier-Stokes solver, are used to demonstrate the possibility of predicting the unsteady aerodynamics and flight trajectories of moving bodies that involve rigid-body collisions.

Meakin, Robert L.↗
NASA NTRS2004

Experimental and Computational Evaluation of Flush-Mounted, S-Duct Inlets

A new high Reynolds number test capability for boundary layer ingesting inlets has been developed for the NASA Langley Research Center 0.3-Meter Transonic Cryogenic Tunnel. Using this new capability. an experimental investigation of four S-duct inlet configurations was conducted. A computational study of one of the inlets was also conducted using a Navier-Stokes solver. The objectives of this investigation were to: 1) develop a new high Reynolds number inlet test capability for flush-mounted inlets; 2) provide a database for CFD tool validation; 3) evaluate the performance of S-duct inlets with large amounts of boundary layer ingestion; and 4) provide a baseline inlet for future inlet flow-control studies. Tests were conducted at Mach numbers from 0.25 to 0.83. Reynolds numbers (based on duct exit diameter) from 5.1 million to a full-scale value of 13.9 million, and inlet mass-flow ratios from 0.39 to 1.58 depending on Mach number. Results of the experimental study indicate that inlet pressure recovery generally decreased and inlet distortion generally increased with increasing Mach number. Except at low Mach numbers, increasing inlet mass-flow increased pressure recovery and increased distortion. Increasing the amount of boundary layer ingestion or ingesting a boundary layer with a distorted profile decreased pressure recovery and increased distortion. Finally, increasing Reynolds number had almost no effect on inlet distortion but increased inlet recovery by about one-half percent at a Mach number near cruise. The computational results captured the inlet pressure recovery and distortion trends with Mach number and inlet mass-flow well: the reversal of the pressure recovery trend with increasing inlet mass-flow at low and high Mach numbers was predicted by CFD. However, CFD results were generally more pessimistic (larger losses) than measured experimentally.

Berrier, Bobby L.↗
NASA NTRS2003

Finite-Rate Ablation Boundary Conditions for Carbon-Phenolic Heat-Shield

A formulation of finite-rate ablation surface boundary conditions, including oxidation, nitridation, and sublimation of carbonaceous material with pyrolysis gas injection, has been developed based on surface species mass conservation. These surface boundary conditions are discretized and integrated with a Navier-Stokes solver. This numerical procedure can predict aerothermal heating, chemical species concentration, and carbonaceous material ablation rate over the heatshield surface of re-entry space vehicles. In this study, the gas-gas and gas-surface interactions are established for air flow over a carbon-phenolic heatshield. Two finite-rate gas-surface interaction models are considered in the present study. The first model is based on the work of Park, and the second model includes the kinetics suggested by Zhluktov and Abe. Nineteen gas phase chemical reactions and four gas-surface interactions are considered in the present model. There is a total of fourteen gas phase chemical species, including five species for air and nine species for ablation products. Three test cases are studied in this paper. The first case is a graphite test model in the arc-jet stream; the second is a light weight Phenolic Impregnated Carbon Ablator at the Stardust re-entry peak heating conditions, and the third is a fully dense carbon-phenolic heatshield at the peak heating point of a proposed Mars Sample Return Earth Entry Vehicle. Predictions based on both finite-rate gas- surface interaction models are compared with those obtained using B' tables, which were created based on the chemical equilibrium assumption. Stagnation point convective heat fluxes predicted using Park's finite-rate model are far below those obtained from chemical equilibrium B' tables and Zhluktov's model. Recession predictions from Zhluktov's model are generally lower than those obtained from Park's model and chemical equilibrium B' tables. The effect of species mass diffusion on predicted ablation rate is also examined.

Chen, Y.-K.↗
NASA NTRS2004

On Improving Efficiency of Differential Evolution for Aerodynamic Shape Optimization Applications

Differential Evolution (DE) is a simple and robust evolutionary strategy that has been provEn effective in determining the global optimum for several difficult optimization problems. Although DE offers several advantages over traditional optimization approaches, its use in applications such as aerodynamic shape optimization where the objective function evaluations are computationally expensive is limited by the large number of function evaluations often required. In this paper various approaches for improving the efficiency of DE are reviewed and discussed. Several approaches that have proven effective for other evolutionary algorithms are modified and implemented in a DE-based aerodynamic shape optimization method that uses a Navier-Stokes solver for the objective function evaluations. Parallelization techniques on distributed computers are used to reduce turnaround times. Results are presented for standard test optimization problems and for the inverse design of a turbine airfoil. The efficiency improvements achieved by the different approaches are evaluated and compared.

Madavan, Nateri K.↗
NASA NTRS1999

X-33 Aerodynamic and Aeroheating Computations for Wind Tunnel and Flight Conditions

This report provides an overview of hypersonic Computational Fluid Dynamics research conducted at the NASA Langley Research Center to support the Phase II development of the X-33 vehicle. The X-33, which is being developed by Lockheed-Martin in partnership with NASA, is an experimental Single-Stage-to-Orbit demonstrator that is intended to validate critical technologies for a full-scale Reusable Launch Vehicle. As part of the development of the X-33, CFD codes have been used to predict the aerodynamic and aeroheating characteristics of the vehicle. Laminar and turbulent predictions were generated for the X 33 vehicle using two finite- volume, Navier-Stokes solvers. Inviscid solutions were also generated with an Euler code. Computations were performed for Mach numbers of 4.0 to 10.0 at angles-of-attack from 10 deg to 48 deg with body flap deflections of 0, 10 and 20 deg. Comparisons between predictions and wind tunnel aerodynamic and aeroheating data are presented in this paper. Aeroheating and aerodynamic predictions for flight conditions are also presented.

Hollis, Brian R.↗
NASA NTRS1997

Transition Effects on Heating in the Wake of a Blunt Body

A series of aerodynamic heating tests was conducted on a 70-deg sphere-cone planetary entry vehicle model in a Mach 10 perfect-gas wind tunnel at freestream Reynolds numbers based on diameter of 8.23x104 to 3.15x105. Surface heating distributions were determined from temperature time-histories measured on the model and on its support sting using thin-film resistance gages. The experimental heating data were compared to computations made using an axisymmetric/2D, laminar, perfect-gas Navier-Stokes solver. Agreement between computational and experimental heating distributions to within, or slightly greater than, the experimental uncertainty was obtained on the forebody and afterbody of the entry vehicle as well as on the sting upstream of the free-shear-layer reattachment point. However, the distributions began to diverge near the reattachment point, with the experimental heating becoming increasingly greater than the computed heating with distance downstream from the reattachment point. It was concluded that this divergence was due to transition of the wake free shear layer just upstream of the reattachment point on the sting.

Hollis, Brian R.↗
NASA NTRS1997

Distributed-Memory Computing With the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA)

The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA), a Navier-Stokes solver, has been modified for use in a parallel, distributed-memory environment using the Message-Passing Interface (MPI) standard. A standard domain decomposition strategy is used in which the computational domain is divided into subdomains with each subdomain assigned to a processor. Performance is examined on dedicated parallel machines and a network of desktop workstations. The effect of domain decomposition and frequency of boundary updates on performance and convergence is also examined for several realistic configurations and conditions typical of large-scale computational fluid dynamic analysis.

Riley, Christopher J.↗
NASA NTRS1996

Development of Semi-Span Model Test Techniques

A computational investigation was performed to support the development of a semi-span model test capability in the NASA Langley Research Center's National Transonic Facility. This capability is desirable for the testing of advanced subsonic transport aircraft at full-scale Reynolds numbers. A state-of-the-art three-dimensional Navier-Stokes solver was used to examine methods to improve the flow over a semi-span configuration. First, a parametric study is conducted to examine the influence of the stand-off height on the flow over the semi-span model. It is found that decreasing the stand-off height, below the maximum fuselage radius, improves the aerodynamic characteristics of the semi-span model. Next, active sidewall boundary layer control techniques are examined. Juncture region blowing jets, upstream tangential blowing, and sidewall suction are found to improve the flow over the aft portion of the semi-span model. Both upstream blowing and suction are found to reduce the sidewall boundary layer separation. The resulting near surface streamline patterns are improved, and found to be quite similar to the full-span results. Both techniques however adversely affect the pitching moment coefficient.

Pulnam, L. Elwood↗
NASA NTRS1996

Comparison of Experimental and Computational Aerothermodynamics of a 70-deg Sphere-Cone

Numerical solutions for hypersonic flows of carbon-dioxide and air around a 70-deg sphere-cone have been computed using an axisymmetric non-equilibrium Navier-Stokes solver. Freestream flow conditions for these computations were equivalent to those obtained in an experimental blunt-body heat-transfer study conducted in a high-enthalpy, hypervelocity expansion tube. Comparisons have been made between the computed and measured surface heat-transfer rates on the forebody and afterbody of the sphere-cone and on the sting which supported the test model. Computed forebody heating rates were within the estimated experimental uncertainties of 10% on the forebody and 15% in the wake except for within the recirculating flow region of the wake.

Hollis, Brian R.↗
NASA NTRS1996

Hypervelocity Heat-Transfer Measurements in an Expansion Tube

A series of experiments has been conducted in the NASA HYPULSE Expansion Tube, in both CO2 and air test gases, in order to obtain data for comparison with computational results and to assess the capability for performing hypervelocity heat-transfer studies in this facility. Heat-transfer measurements were made in both test gases on 70 deg sphere-cone models and on hemisphere models of various radii. HYPULSE freestream flow conditions in these test gases were found to be repeatable to within 3-10%, and aerothermodynamic test times of 150 microsec in CO2 and 125 microsec in air were identified. Heat-transfer measurement uncertainty was estimated to be 10-15%. Comparisons were made with computational results from the non-equilibrium Navier-Stokes solver NEQ2D. Measured and computed heat-transfer rates agreed to within 10% on the hemispheres and on the sphere-cone forebodies, and to within 10% in CO2 and 25% in air on the afterbodies and stings of the sphere-cone models.

Hollis, Brian R.↗
NASA NTRS2004

Preliminary Convective-Radiative Heating Environments for a Neptune Aerocapture Mission

Convective and radiative heating environments have been computed for a three-dimensional ellipsled configuration which would perform an aerocapture maneuver at Neptune. This work was performed as part of a one-year Neptune aerocapture spacecraft systems study that also included analyses of trajectories, atmospheric modeling, aerodynamics, structural design, and other disciplines. Complementary heating analyses were conducted by separate teams using independent sets of aerothermodynamic modeling tools (i.e. Navier-Stokes and radiation transport codes). Environments were generated for a large 5.50 m length ellipsled and a small 2.88 m length ellipsled. Radiative heating was found to contribute up to 80% of the total heating rate at the ellipsled nose depending on the trajectory point. Good agreement between convective heating predictions from the two Navier-Stokes solvers was obtained. However, the radiation analysis revealed several uncertainties in the computational models employed in both sets of codes, as well as large differences between the predicted radiative heating rates.

Hollis, Brian R.↗
NASA NTRS1992

Characteristics of the Shuttle Orbiter Leeside Flow During A Reentry Condition

A study of the leeside flow characteristics of the Shuttle Orbiter is presented for a reentry flight condition. The flow is computed using a point-implicit, finite-volume scheme known as the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA). LAURA is a second-order accurate, laminar Navier-Stokes solver, incorporating finite-rate chemistry with a radiative equilibrium wall temperature distribution and finite-rate wall catalysis. The resulting computational solution is analyzed in terms of salient flow features and the surface quantities are compared with flight data.

Kleb, William L.↗
NASA NTRS2004

AIC Computations Using Navier-Stokes Equations on Single Image Supercomputers For Design Optimization

A procedure to accurately generate AIC using the Navier-Stokes solver including grid deformation is presented. Preliminary results show good comparisons between experiment and computed flutter boundaries for a rectangular wing. A full wing body configuration of an orbital space plane is selected for demonstration on a large number of processors. In the final paper the AIC of full wing body configuration will be computed. The scalability of the procedure on supercomputer will be demonstrated.

Guruswamy, Guru↗
NASA NTRS2004

On Improving Efficiency of Differential Evolution for Aerodynamic Shape Optimization Applications

Differential Evolution (DE) is a simple and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems. Although DE offers several advantages over traditional optimization approaches, its use in applications such as aerodynamic shape optimization where the objective function evaluations are computationally expensive is limited by the large number of function evaluations often required. In this paper various approaches for improving the efficiency of DE are reviewed and discussed. These approaches are implemented in a DE-based aerodynamic shape optimization method that uses a Navier-Stokes solver for the objective function evaluations. Parallelization techniques on distributed computers are used to reduce turnaround times. Results are presented for the inverse design of a turbine airfoil. The efficiency improvements achieved by the different approaches are evaluated and compared.

Madavan, Nateri K.↗