Search NASA⌕ Search

SEARCH · Search NASA

Results for “flow 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 379 records · Page 21

Parallelization of the Implicit RPLUS Algorithm

The multiblock reacting Navier-Stokes flow-solver RPLUS2D was modified for parallel implementation. Results for non-reacting flow calculations of this code indicate parallelization efficiencies greater than 84% are possible for a typical test problem. Results tend to improve as the size of the problem increases. The convergence rate of the scheme is degraded slightly when additional artificial block boundaries are included for the purpose of parallelization. However, this degradation virtually disappears if the solution is converged near to machine zero. Recommendations are made for further code improvements to increase efficiency, correct bugs in the original version, and study decomposition effectiveness.

Orkwis, Paul D.↗

Numerical Solutions for the CAWAPI Configuration on Unstructured Grids at NASA LaRC, United States

In support of the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) with its goal of improving the Technology Readiness Level of flow solvers by comparing results with measured F-16XL-1 flight data, NASA Langley employed the TetrUSS unstructured grid solver, USM3D, to obtain solutions for all seven flight conditions of interest. A newly available solver version that incorporates a number of turbulence models, including the two-equation linear and non-linear k- , was used in this study. As a first test, a choice was made to utilize only a single grid resolution with the solver for the simulation of the different flight conditions. Comparisons are presented with three turbulence models in USM3D, flight data for surface pressure, boundary-layer profiles, and skin-friction distribution, as well as limited predictions from other solvers. A result of these comparisons is that the USM3D solver can be used in an engineering environment to predict vortex-flow physics on a complex configuration at flight Reynolds numbers with a two-equation linear k- turbulence model.

Lamar, John E.↗

Capsule_Grid User Guide

Software utility CAPSULE_GRID and several related procedures that automate grid generation for computational fluid dynamics calculations about axisymmetric atmospheric entry vehicles are described as a supplement to the documentation within the source codes. This application program and its ancillary software utilities are all in the public domain at the web site https://software.nasa.gov/software/category/all/arc/1/cfdtools with the exception of the 2D hyperbolic volume gridding step for which scripts driving either GRIDGEN1 or HYPGEN2 are available. CAPSULE_GRID was developed at NASA Ames Research Center in support of the DPLR3 real gas flow solver. Both forebodies (only) and full bodies are treated, for axisymmetric 2D flow calculations or for 3D calculations at angle of attack. The expected application is to the smooth outer mold line of a hypersonic atmospheric entry vehicle, for which a single layer of point-matched structured grid blocks normally suffices. Control files for a variety of representative cases are included with illustrations of results, and some insights into the subtleties of curvature-based grid point distribution are provided. Gridding of a capsule on a sting is supported. Additional procedures for treating arbitrarily long 2D wakes, prompted by asteroid studies and relevant to spacecraft such as Stardust, are also documented. Most recently, treating an asymmetric forebody (including an off-center nose, with or without an aft body) has been incorporated as an extra step.

spacecraft↗

Numerical Simulations of A Conceptual MSR-EES Shoulder Recession

The present study demonstrates our in-house material response solver, Icarus's capability to simulate the ablative conditions, including pyrolysis effects due to the interaction between hypersonic boundary layers and the thermal protection system (TPS). A conceptual aeroshell shoulder design, which undergoes mission-relevant flow and material conditions, is selected for the demonstration purpose. LAURA, a structured flow solver, is used to solve flow around the shoulder at several trajectory points of a flight path. The aerothermal dataset obtained from LAURA is used to enforce boundary conditions on the aeroshell wall to simulate ablation processes. The two-layered material system is stacked with HEEET (TPS) and Aluminum (actual material). A set of base parameters, such as the angle of material orientation with respect to the flow direction, convective heat-transfer co-efficient, and aerothermal boundary conditions on the rear end of the HEEET material layer, is selected to conduct a base-case simulation. A small amount of recession was observed, indicating that the design should be fine to go through the selected flow and material conditions. Furthermore, the parameters mentioned above are individually altered and are observed to affect the shoulder design's recession compared to the base case. Verification of our Icarus setup was also carried out using a one-dimensional grid and FIAT, a one-dimensional material-response solver, to build credibility for our results. As the current work is an uncoupled fluid-material-response simulation, a local sharp mesh deformation in the recessed surface near the shoulder corner is observed.

Prakash Shrestha↗

Three-D Flow Analysis of the Alternate SSME HPOT TAD

This paper describes the results of numerical flow analyses performed in support of design development of the Space Shuttle Main Engine Alternate High Pressure Oxidizer Turbine Turn-around duct (TAD). The flow domain has been modeled using a 3D, Navier-Stokes, general purpose flow solver. The goal of this effort is to achieve an alternate TAD exit flow distribution which closely matches that of the baseline configuration. 3D Navier Stokes CFD analyses were employed to evaluate numerous candidate geometry modifications to the TAD flowpath in order to achieve this goal. The design iterations are summarized, as well as a description of the computational model, numerical results and the conclusions based on these calculations.

Kubinski, Cheryl A.↗

Computation of Sound Generated by Flow Over a Circular Cylinder: An Acoustic Analogy Approach

The sound generated by viscous flow past a circular cylinder is predicted via the Lighthill acoustic analogy approach. The two dimensional flow field is predicted using two unsteady Reynolds-averaged Navier-Stokes solvers. Flow field computations are made for laminar flow at three Reynolds numbers (Re = 1000, Re = 10,000, and Re = 90,000) and two different turbulent models at Re = 90,000. The unsteady surface pressures are utilized by an acoustics code that implements Farassat's formulation 1A to predict the acoustic field. The acoustic code is a 3-D code - 2-D results are found by using a long cylinder length. The 2-D predictions overpredict the acoustic amplitude; however, if correlation lengths in the range of 3 to 10 cylinder diameters are used, the predicted acoustic amplitude agrees well with experiment.

Brentner, Kenneth S.↗

High Fidelity Simulations of Unsteady Flow through Turbopumps and Flowliners

High fidelity computations were carried out to analyze the orbiter LH2 feedline flowliner. Computations were performed on the Columbia platform which is a 10,240-processor supercluster consisting of 20 Altix nodes with 512 processor each. Various computational models were used to characterize the unsteady flow features in the turbopump, including the orbiter Low-Pressure-Fuel-Turbopump (LPFTP) inducer, the orbiter manifold and a test article used to represent the manifold. Unsteady flow originating from the orbiter LPFTP inducer is one of the major contributors to the high frequency cyclic loading that results in high cycle fatigue damage to the gimbal flowliners just upstream of the LPFTP. The flow fields for the orbiter manifold and representative test article are computed and analyzed for similarities and differences. The incompressible Navier-Stokes flow solver INS3D, based on the artificial compressibility method, was used to compute the flow of liquid hydrogen in each test article.

Kiris, Cetin C.↗

High Fidelity Simulations for Unsteady Flow Through the Orbiter LH2 Feedline Flowliner

High fidelity computations were carried out to analyze the orbiter M2 feedline flowliner. Various computational models were used to characterize the unsteady flow features in the turbopump, including the orbiter Low-Pressure-Fuel-Turbopump (LPFTP) inducer, the orbiter manifold and a test article used to represent the manifold. Unsteady flow originating from the orbiter LPFTP inducer is one of the major contributors to the high frequency cyclic loading that results in high cycle fatigue damage to the gimbal flowliners just upstream of the LPFTP. The flow fields for the orbiter manifold and representative test article are computed and analyzed for similarities and differences. An incompressible Navier-Stokes flow solver INS3D, based on the artificial compressibility method, was used to compute the flow of liquid hydrogen in each test article.

Kiris, Cetin C.↗

Enabling Metric-Based Mesh Adaptation for Advanced Compressible Flow Simulations Using US3D

This work presents a metric-based mesh adaptation capability added to the flow solver US3D that is used to improve the accuracy of atmospheric entry simulations. Flow past an entry vehicle in free-flight is usually unsteady and characterized by anisotropic flow features like a strong bow shock and small-scale isotropic flow features like turbulent eddies in the wake. The results presented in this paper demonstrate that we are currently able to iteratively apply metric-based mesh adaptation to improve the solution accuracy by simultaneously targeting isotropic and anisotropic flow features. A simplified workflow of running CFD simulations for entry vehicles is presented where we start off with a coarse initial mesh to compute a first approximation of the flow. This coarse mesh is then iteratively refined using the newly developed metric to improve the resolution of the CFD computation. Furthermore, we demonstrate that we are able to adapt large mesh sizes in a timely fashion which allows us to start thinking about applying this adaptive framework to unsteady atmospheric entry CFD simulations.

Mesh adaptation↗

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗

Aeroacoustic Simulation of a Nose Landing Gear in an Open Jet Facility Using FUN3D

Numerical simulations have been performed for a partially-dressed, cavity-closed nose landing gear configuration that was tested in NASA Langley s closed-wall Basic Aerodynamic Research Tunnel (BART) and in the University of Florida s open-jet acoustic facility known as UFAFF. The unstructured-grid flow solver, FUN3D, developed at NASA Langley Research center is used to compute the unsteady flow field for this configuration. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these computations. Time-averaged and instantaneous solutions compare favorably with the measured data. Unsteady flowfield data obtained from the FUN3D code are used as input to a Ffowcs Williams-Hawkings noise propagation code to compute the sound pressure levels at microphones placed in the farfield. Significant improvement in predicted noise levels is obtained when the flowfield data from the open jet UFAFF simulations is used as compared to the case using flowfield data from the closed-wall BART configuration.

Vatsa, Veer N.↗

Effect of Inflow Boundary Conditions on Hovering Tilt-Rotor Flows

A study on the effect of inflow turbulence boundary conditions on the local flow on and around a body in flight has been carried out. The study has been carried out using OVERFLOW2 flow solver using the default Spalart-Allmaras turbulence model in OVERFLOW. Many OVERFLOW turbulent flow simulations have been reported using the SA-fv3 model1 over the years. The present study demonstrates that the turbulence levels imposed as a boundary condition (b.c.) in the far-field as implemented in OVERFLOW are not correct. In fact, very low level of turbulence at the far-field boundary as implemented in OVERFLOW results in the SA-fv3 model predicting transition-like profiles on a given body. By choosing sufficiently high levels of the Reynolds number of turbulence, Rt, as a boundary condition, this anomalous behavior of SA-fv3 model is eliminated. Since numerous papers using OVERFLOW have been presented in the literature including that by the author2 using low level of inflow turbulence, it will be beneficial to the CFD community at large and in particular to the OVERFLOW community to understand this effect of the inflow b.c. in the SA-fv3 model. Various results reported over the years in domains such as high-lift applications, drag prediction applications and rotorcraft flow applications can be revisited using the new turbulence inflow boundary condition prescription as suggested in this study. It has been demonstrated in the literature that the effect of inflow turbulence level on the downstream development of flow is too important to be ignored. Demonstrative results are shown and compared with experiment for the case of hovering XV-15 rotor flow at a tip Reynolds number of 4.9 x 10(exp 6) and Mach number of 0.69 in Figs. 1 and 2 below. Fig. 1(a,b) shows skin friction predictions at radial station of r/R = 0.28 corresponding to 5 different inflow boundary conditions obtained with the SA-fv3 model. As the inflow turbulence is increased monotonically, from a value of Rt = 0.1 to 100.0 as shown in Fig. 1, a fully developed turbulent skin friction profile is realized progressively as shown in Fig. 1(b). Similarly skin friction results shown in Fig. 2(a,b) corresponding to radial stations of r/R = 0.72 reinforce the postulate made here that as the inflow turbulence level is increased gradually, fully turbulent on-the-blade simulation is realized uniformly. In the full paper extensive results on skin friction, pressure and velocity profiles as well as the Figure of Merit will be presented and compared with experiment that will illustrate the substantial effect of the inflow boundary conditions on the development of hovering rotor flow downstream.

Kaul, Upender K.↗

Aeroacoustic Simulation of Nose Landing Gear on Adaptive Unstructured Grids With FUN3D

Numerical simulations have been performed for a partially-dressed, cavity-closed nose landing gear configuration that was tested in NASA Langley s closed-wall Basic Aerodynamic Research Tunnel (BART) and in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D, developed at NASA Langley Research center, is used to compute the unsteady flow field for this configuration. Starting with a coarse grid, a series of successively finer grids were generated using the adaptive gridding methodology available in the FUN3D code. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these computations. Time-averaged and instantaneous solutions obtained on these grids are compared with the measured data. In general, the correlation with the experimental data improves with grid refinement. A similar trend is observed for sound pressure levels obtained by using these CFD solutions as input to a FfowcsWilliams-Hawkings noise propagation code to compute the farfield noise levels. In general, the numerical solutions obtained on adapted grids compare well with the hand-tuned enriched fine grid solutions and experimental data. In addition, the grid adaption strategy discussed here simplifies the grid generation process, and results in improved computational efficiency of CFD simulations.

Vatsa, Veer N.↗

Reconfigurable Model Execution in the OpenMDAO Framework

NASA's OpenMDAO framework facilitates constructing complex models and computing their derivatives for multidisciplinary design optimization. Decomposing a model into components that follow a prescribed interface enables OpenMDAO to assemble multidisciplinary derivatives from the component derivatives using what amounts to the adjoint method, direct method, chain rule, global sensitivity equations, or any combination thereof, using the MAUD architecture. OpenMDAO also handles the distribution of processors among the disciplines by hierarchically grouping the components, and it automates the data transfer between components that are on different processors. These features have made OpenMDAO useful for applications in aircraft design, satellite design, wind turbine design, and aircraft engine design, among others. This paper presents new algorithms for OpenMDAO that enable reconfigurable model execution. This concept refers to dynamically changing, during execution, one or more of: the variable sizes, solution algorithm, parallel load balancing, or set of variables-i.e., adding and removing components, perhaps to switch to a higher-fidelity sub-model. Any component can reconfigure at any point, even when running in parallel with other components, and the reconfiguration algorithm presented here performs the synchronized updates to all other components that are affected. A reconfigurable software framework for multidisciplinary design optimization enables new adaptive solvers, adaptive parallelization, and new applications such as gradient-based optimization with overset flow solvers and adaptive mesh refinement. Benchmarking results demonstrate the time savings for reconfiguration compared to setting up the model again from scratch, which can be significant in large-scale problems. Additionally, the new reconfigurability feature is applied to a mission profile optimization problem for commercial aircraft where both the parametrization of the mission profile and the time discretization are adaptively refined, resulting in computational savings of roughly 10% and the elimination of oscillations in the optimized altitude profile.

Hwang, John T.↗

High Energy Boundary Conditions for a Cartesian Mesh Euler Solver

Inlets and exhaust nozzles are often omitted or fared over in aerodynamic simulations of aircraft due to the complexities involving in the modeling of engine details such as complex geometry and flow physics. However, the assumption is often improper as inlet or plume flows have a substantial effect on vehicle aerodynamics. A tool for specifying inlet and exhaust plume conditions through the use of high-energy boundary conditions in an established inviscid flow solver is presented. The effects of the plume on the flow fields near the inlet and plume are discussed.

Pandya, Shishir A.↗

Application of the computational aeroacoustics method to an advanced counterrotating propfan configuration

The acoustical near field of a counterrotating propeller configuration of unducted fan (UDF) type is presently predicted by means of a numerical methodology employing the transonic flow field generated by the three-dimensional numerical flow solver designated SSTAGE. This numerical approach simulates the time-averaged, three-dimensional inviscid flow field of the UDF by means of a four-stage Runge-Kutta integration scheme which marches the equations forward in time toward an asymptotic limit. The near field pressure time history is obtained directly from the flow field solution; by converting these histories with a Fourier series representation, specific data on the acoustic signal frequency spectra can be obtained.

Kim, Jinhan↗

Progress in high-lift aerodynamic calculations

The current work presents progress in the effort to numerically simulate the flow over high-lift aerodynamic components, namely, multi-element airfoils and wings in either a take-off or a landing configuration. The computational approach utilizes an incompressible flow solver and an overlaid chimera grid approach. A detailed grid resolution study is presented for flow over a three-element airfoil. Two turbulence models, a one-equation Baldwin-Barth model and a two equation k-omega model are compared. Excellent agreement with experiment is obtained for the lift coefficient at all angles of attack, including the prediction of maximum lift when using the two-equation model. Results for two other flap riggings are shown. Three-dimensional results are presented for a wing with a square wing-tip as a validation case. Grid generation and topology is discussed for computing the flow over a T-39 Sabreliner wing with flap deployed and the initial calculations for this geometry are presented.

Rogers, Stuart E.↗

Investigation of Boundary Layer Development in Multi-Stage Low Pressure

A review of recent literature has been made in the areas of unsteady measurement in turbomachinery and LP turbine design. This has included a detailed review of methods for extracting quantitative shear stress data from surface mounted hot-film arrays. Research goals were identified on the basis of the literature survey and following interaction with the GE and NASA contracts. Progress has also been made in the detailed planning of the research program. An array of hot-film sensors has been mounted at mid-span on a second stage nozzle blade ready for installation in the current turbine buildup. A new technique of applying the hot-film array using double sided adhesive is being tested. The temperature in the vicinity of the hot-films will be monitored using a thermocouple attached to the blade pressure surface. This will increase the accuracy of the data reduction process compared with the methods applied by Halstead et at. (1995) and others. Data obtained from this first set of sensors, along with data from the hot-wire investigation will help to determine the best locations for surface hot-film instrumentation for subsequent builds of the turbine. Other LSRT instrumentation issues which have been addressed include calibration of the torque measurement system and procurement of custom design hot-wire sensors required to compensate for the slope flow path. Preliminary attempts at modeling the boundary layer transition process have been completed. A steady, coupled viscous-inviscid flow solver (the MISES code of Drela and Youngren) was used to calculate the steady-flow boundary layer transition location.

Source record↗