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 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Incompressible Navier-Stokes Computations in Aerospace and Beyond

Challenges in applying CFD to real world problems will be discussed using incompressible flow examples. Despite the remarkable progress made in CFD technology during the past two decades, new challenges are ahead of us. Some of these will be discussed in conjunction with the incompressible Navier-Stokes solver development and their applications. Computed examples include wing tip vortex formation and propagation, flow simulation of an advanced rocket pump, and an extension of the similar technology to biofluid analysis and design. Numerical issues dealing with these applications will be discussed both from algorithm development and from application point of view. The material will be printed in a bound volume after the symposium.

Kwak, Dochan↗

NASA NTRS2002

Turbomachinery Airfoil Design Optimization Using Differential Evolution

An aerodynamic design optimization procedure that is based on a evolutionary algorithm known at Differential Evolution is described. Differential Evolution is a simple, fast, and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems, including highly nonlinear systems with discontinuities and multiple local optima. The method is combined with a Navier-Stokes solver that evaluates the various intermediate designs and provides inputs to the optimization procedure. An efficient constraint handling mechanism is also incorporated. Results are presented for the inverse design of a turbine airfoil from a modern jet engine and compared to earlier methods. The capability of the method to search large design spaces and obtain the optimal airfoils in an automatic fashion is demonstrated. Substantial reductions in the overall computing time requirements are achieved by using the algorithm in conjunction with neural networks.

Madavan, Nateri K.↗

NASA NTRS1998

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗

NASA NTRS1998

Aerodynamic Influence of a Slat on a Flapped Wing

Three high-lift configurations were computationally studied to assess the aerodynamic influence of slats. A flapped wing was simulated with no slat, a full-span slat, and a three quarter-span slat at a chord based Reynolds number of 3.7 million. The flows were computed using a compressible Navier-Stokes solver on structured grids with the Spalart-Allmaras turbulence model. All cases were compared with experimental data to validate the approach. The slats not only increase the lift generated by the wing but alter the topology of the flowfield considerably. The changes in the flow give insight into the working of a slat and contribute to a better understanding of high-lift flows in general.

Mathias, Donovan L.↗

NASA NTRS1999

Improving the Unsteady Aerodynamic Performance of Transonic Turbines using Neural Networks

A recently developed neural net-based aerodynamic design procedure is used in the redesign of a transonic turbine stage to improve its unsteady aerodynamic performance. The redesign procedure used incorporates the advantages of both traditional response surface methodology and neural networks by employing a strategy called parameter-based partitioning of the design space. Starting from the reference design, a sequence of response surfaces based on both neural networks and polynomial fits are constructed to traverse the design space in search of an optimal solution that exhibits improved unsteady performance. The procedure combines the power of neural networks and the economy of low-order polynomials (in terms of number of simulations required and network training requirements). A time-accurate, two-dimensional, Navier-Stokes solver is used to evaluate the various intermediate designs and provide inputs to the optimization procedure. The procedure yielded a modified design that improves the aerodynamic performance through small changes to the reference design geometry. These results demonstrate the capabilities of the neural net-based design procedure, and also show the advantages of including high-fidelity unsteady simulations that capture the relevant flow physics in the design optimization process.

Rai, Man Mohan↗

NASA NTRS2002

Turbomachinery Airfoil Design Optimization Using Differential Evolution

An aerodynamic design optimization procedure that is based on a evolutionary algorithm known at Differential Evolution is described. Differential Evolution is a simple, fast, and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems, including highly nonlinear systems with discontinuities and multiple local optima. The method is combined with a Navier-Stokes solver that evaluates the various intermediate designs and provides inputs to the optimization procedure. An efficient constraint handling mechanism is also incorporated. Results are presented for the inverse design of a turbine airfoil from a modern jet engine. The capability of the method to search large design spaces and obtain the optimal airfoils in an automatic fashion is demonstrated. Substantial reductions in the overall computing time requirements are achieved by using the algorithm in conjunction with neural networks.

Madavan, Nateri K.↗

NASA NTRS2002

An Evaluation of Parameters Influencing Jet Mixing Using the WIND Navier-stokes Code

The WIND code, a Reynolds-averaged Navier-Stokes solver used for a variety of aerospace flow simulations, was investigated for a Mach 2 nozzle at a series of nozzle stagnation temperatures. Comparisons of WIND calculations are made to experimental measurements of axial velocity, Mach number, and stagnation temperature along the jet centerline. The primary objective was to investigate the capabilities of the two-equation turbulence models available in WIND, version 4.0, for the analysis of heated supersonic nozzle flows. The models examined were the Menter Shear Stress Transport (SST) model and the Chien k-epsilon model, with and without the compressibility correction due to Sarkar. It was observed that all of the turbulence models investigated produced solutions that did not agree well with the experimental measurements. The effects of freestream Mach number and turbulent Prandtl number specifications were also investigated.

Dembowski, Mary Ann↗

NASA NTRS2002

Heat Transfer in a Complex Trailing Edge Passage for a High Pressure Turbine Blade: Simulation Results - Part 2:

A combined experimental and numerical study to investigate the heat transfer distribution in a complex blade trailing edge passage was conducted. The geometry consists of a two pass serpentine passage with taper toward the trailing edge, as well as from hub to tip. The upflow channel has an average aspect ratio of roughly 14:1, while the exit passage aspect ratio is about 5:1. The upflow channel is split in an interrupted way and is smooth on the trailing edge side of the split and turbulated on the other side. A turning vane is placed near the tip of the upflow channel. Reynolds numbers in the range of 31,000 to 61,000, based on inlet conditions, were simulated numerically. The simulation was performed using the Glenn-HT code, a full three-dimensional Navier-Stokes solver using the Wilcox k-omega turbulence model. A structured multi-block grid is used with approximately 4.5 million cells and average y+ values on the order of unity. Pressure and heat transfer distributions are presented with comparison to the experimental data. While there are some regions with discrepancies, in general the agreement is very good for both pressure and heat transfer.

Rigby, David L.↗

NASA NTRS2002

Graphite Ablation and Thermal Response Simulation Under Arc-Jet Flow Conditions

The Two-dimensional Implicit Thermal Response and Ablation program, TITAN, was developed and integrated with a Navier-Stokes solver, GIANTS, for multidimensional ablation and shape change simulation of thermal protection systems in hypersonic flow environments. The governing equations in both codes are demoralized using the same finite-volume approximation with a general body-fitted coordinate system. Time-dependent solutions are achieved by an implicit time marching technique using Gauess-Siedel line relaxation with alternating sweeps. As the first part of a code validation study, this paper compares TITAN-GIANTS predictions with thermal response and recession data obtained from arc-jet tests recently conducted in the Interaction Heating Facility (IHF) at NASA Ames Research Center. The test models are graphite sphere-cones. Graphite was selected as a test material to minimize the uncertainties from material properties. Recession and thermal response data were obtained from two separate arc-jet test series. The first series was at a heat flux where graphite ablation is mainly due to sublimation, and the second series was at a relatively low heat flux where recession is the result of diffusion-controlled oxidation. Ablation and thermal response solutions for both sets of conditions, as calculated by TITAN-GIANTS, are presented and discussed in detail. Predicted shape change and temperature histories generally agree well with the data obtained from the arc-jet tests.

Chen, Y.-K.↗

NASA NTRS2003

Aerodynamic Shape Optimization Using Hybridized Differential Evolution

An aerodynamic shape optimization method that uses an evolutionary algorithm known at Differential Evolution (DE) in conjunction with various hybridization strategies is described. DE is a simple and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems. Various hybridization strategies for DE are explored, including the use of neural networks as well as traditional local search methods. A Navier-Stokes solver is used to evaluate the various intermediate designs and provide inputs to the hybrid DE optimizer. The method is implemented on distributed parallel computers so that new designs can be obtained within reasonable turnaround times. Results are presented for the inverse design of a turbine airfoil from a modern jet engine. (The final paper will include at least one other aerodynamic design application). The capability of the method to search large design spaces and obtain the optimal airfoils in an automatic fashion is demonstrated.

Madavan, Nateri K.↗

NASA NTRS2002

Higher Order Time Integration Schemes for the Unsteady Navier-Stokes Equations on Unstructured Meshes

The rapid increase in available computational power over the last decade has enabled higher resolution flow simulations and more widespread use of unstructured grid methods for complex geometries. While much of this effort has been focused on steady-state calculations in the aerodynamics community, the need to accurately predict off-design conditions, which may involve substantial amounts of flow separation, points to the need to efficiently simulate unsteady flow fields. Accurate unsteady flow simulations can easily require several orders of magnitude more computational effort than a corresponding steady-state simulation. For this reason, techniques for improving the efficiency of unsteady flow simulations are required in order to make such calculations feasible in the foreseeable future. The purpose of this work is to investigate possible reductions in computer time due to the choice of an efficient time-integration scheme from a series of schemes differing in the order of time-accuracy, and by the use of more efficient techniques to solve the nonlinear equations which arise while using implicit time-integration schemes. This investigation is carried out in the context of a two-dimensional unstructured mesh laminar Navier-Stokes solver.

Jothiprasad, Giridhar↗

NASA NTRS2002

Progress in Unsteady Turbopump Flow Simulations Using Overset Grid Systems

This viewgraph presentation provides information on unsteady flow simulations for the Second Generation RLV (Reusable Launch Vehicle) baseline turbopump. Three impeller rotations were simulated by using a 34.3 million grid points model. MPI/OpenMP hybrid parallelism and MLP shared memory parallelism has been implemented and benchmarked in INS3D, an incompressible Navier-Stokes solver. For RLV turbopump simulations a speed up of more than 30 times has been obtained. Moving boundary capability is obtained by using the DCF module. Scripting capability from CAD geometry to solution is developed. Unsteady flow simulations for advanced consortium impeller/diffuser by using a 39 million grid points model are currently underway. 1.2 impeller rotations are completed. The fluid/structure coupling is initiated.

Kiris, Cetin C.↗

NASA NTRS2003

Numerical and Experimental Determination of the Geometric Far Field for Round Jets

To reduce ambiguity in the reporting of far field jet noise, three round jets operating at subsonic conditions have recently been studied at the NASA Glenn Research Center. The goal of the investigation was to determine the location of the geometric far field both numerically and experimentally. The combination of the WIND Reynolds-Averaged Navier-Stokes solver and the MGBK jet noise prediction code was used for the computations, and the experimental data was collected in the Aeroacoustic Propulsion Laboratory. While noise sources are distributed throughout the jet plume, at great distances from the nozzle the noise will appear to be emanating from a point source and the assumption of linear propagation is valid. Closer to the jet, nonlinear propagation may be a problem, along with the known geometric issues. By comparing sound spectra at different distances from the jet, both from computational methods that assume linear propagation, and from experiments, the contributions of geometry and nonlinearity can be separately ascertained and the required measurement distance for valid experiments can be established. It is found that while the shortest arc considered here (approx. 8D) was already in the geometric far field for the high frequency sound (St greater than 2.0), the low frequency noise due to its extended source distribution reached the geometric far field at or about 50D. It is also found that sound spectra at far downstream angles does not strictly scale on Strouhal number, an observation that current modeling does not capture.

Koch, L. Danielle↗

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.↗