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.

SourceResearch connectionFrom yearThrough year

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

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Two-Phase Turbulence Statistics from High Fidelity Dispersed Droplet Flow Simulations in a Pressurized Water Reactor (PWR) Sub-Channel with Mixing Vanes

In the dispersed flow film boiling regime (DFFB), which exists under post-LOCA (loss-of-coolant accident) conditions in pressurized water reactors (PWRs), there is a complex interplay between droplet dynamics and turbulence in the surrounding steam. Experiments have accredited particular significance to droplet collision with the spacer-grids and mixing vane structures and their consequent positive feedback to the heat transfer recorded in the immediate downstream vicinity. Enabled by high-performance computing (HPC) systems and a massively parallel finite element-based flow solver—PHASTA (Parallel Hierarchic Adaptive Stabilized Transient Analysis)—this work presents high fidelity interface capturing, two-phase, adiabatic simulations in a PWR sub-channel with spacer grids and mixing vanes. Selected flow conditions for the simulations are informed by the experimental data found in the literature, including the steam Reynolds number and collision Weber number (Wec={40,80}), and are characteristic of the DFFB regime. Data were collected from the simulations at an unprecedented resolution, which provides detailed insights into the continuous phase turbulence statistics, highlighting the effects of the presence of droplets and the comparative effect of different Weber numbers on turbulence in the surrounding steam. Further, axial evolution of droplet dynamics was analyzed through cross-sectionally averaged quantities, including droplet volume, surface area and Sauter mean diameter (SMD). The downstream SMD values agree well with the existing empirical correlations for the selected range of Wec. The high-resolution data repository from the simulations herein is expected to be of significance to guide model development for system-level thermal hydraulic codes.

Saini, Nadish↗

NASA NTRS1991

The effects of self-generated and applied magnetic fields on the computation of flow over a Mars return aerobrake

A CFD technique is developed to calculate the electromagnetic phenomena simultaneously with the fluid flow in the shock layer over an axisymmetric blunt body in a thermal-equilibrium chemical-nonequilibrium environment. The flowfield is solved using an explicit time-marching, first-order spatially accurate scheme. The electromagnetic phenomena are coupled to the real-gas flow solver through an iterative procedure. The electromagnetic terms introduce a strong stiffness, which was overcome by using significantly smaller time steps for the electromagnetic conservation equation. The technique is applied in calculating the flow over a Mars return aerobrake vehicle entering the Earth's atmosphere. For the case where no external field is applied, the electromagnetic effects have little impact on the flowfield.

Palmer, Grant↗

NASA NTRS1995

Unstructured viscous flow solution using adaptive hybrid grids

A three dimensional finite volume scheme based on hybrid grids containing both tetrahedral and hexahedral cells is presented. The application to hybrid grids offers the possibility to combine the flexibility of tetrahedral meshes with the accuracy of hexahedral grids. An algorithm to compute a dual mesh for the entire computational domain was developed. The dual mesh technique guarantees conservation in the whole flow field even at interfaces between hexahedral and tetrahedral domains and enables the employment of an accurate upwind flow solver. The hybrid mesh can be adapted to the solution by dividing cells in areas of insufficient resolution. The method is tested on different viscous and inviscid cases for hypersonic, transonic and subsonic flows.

Galle, Martin↗

NASA NTRS1995

Evaluation of Aircraft Platforms for SOFIA by Computational Fluid Dynamics

The selection of an airborne platform for the Stratospheric Observatory for Infrared Astronomy (SOFIA) is based not only on economic cost, but technical criteria, as well. Technical issues include aircraft fatigue, resonant characteristics of the cavity-port shear layer, aircraft stability, the drag penalty of the open telescope bay, and telescope performance. Recently, two versions of the Boeing 747 aircraft, viz., the -SP and -200 configurations, were evaluated by computational fluid dynamics (CFD) for their suitability as SOFIA platforms. In each configuration the telescope was mounted behind the wings in an open bay with nearly circular aperture. The geometry of the cavity, cavity aperture, and telescope was identical in both platforms. The aperture was located on the port side of the aircraft and the elevation angle of the telescope, measured with respect to the vertical axis, was 500. The unsteady, viscous, three-dimensional, aerodynamic and acoustic flow fields in the vicinity of SOFIA were simulated by an implicit, finite-difference Navier-Stokes flow solver (OVERFLOW) on a Chimera, overset grid system. The computational domain was discretized by structured grids. Computations were performed at wind-tunnel and flight Reynolds numbers corresponding to one free-stream flow condition (M = 0.85, angle of attack alpha = 2.50, and sideslip angle beta = 0 degrees). The computational domains consisted of twenty-nine(29) overset grids in the wind-tunnel simulations and forty-five(45) grids in the simulations run at cruise flight conditions. The maximum number of grid points in the simulations was approximately 4 x 10(exp 6). Issues considered in the evaluation study included analysis of the unsteady flow field in the cavity, the influence of the cavity on the flow across empennage surfaces, the drag penalty caused by the open telescope bay, and the noise radiating from cavity surfaces and the cavity-port shear layer. Wind-tunnel data were also available to compare to the CFD results; the data permitted an assessment of CFD as a design tool for the SOFIA program.

Klotz, S. P.↗

NASA NTRS2020

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5 angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for proper validation of the CFD.

Maldonado, Daniel↗

NASA NTRS2020

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5o angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for more thorough validation of the CFD.

Maldonado, Daniel↗

NASA NTRS1996

Implementation of a two-equation k-omega turbulence model in NPARC

The implementation of a two-equation k-omega turbulence model into the NPARC flow solver is described. Motivation for the selection of this model is given, major code modifications are outlined, new imputs to the code are described, and results are presented for several validation cases: an incompressible flow over a smooth flat plate, a subsonic diffuser flow, and a shock-induced separated flow. Comparison of results with the k-epsilon model indicate that the k-omega model predicts simple flows equally well whereas, for adverse pressure gradient flows, the k-omega model outperforms the other turbulence models in NPARC.

Yoder, Dennis A.↗

NASA NTRS2005

High Reynolds Number Investigation of a Flush-Mounted, S-Duct Inlet With Large Amounts of Boundary Layer Ingestion

An experimental investigation of a flush-mounted, S-duct inlet with large amounts of boundary layer ingestion has been conducted at Reynolds numbers up to full scale. The study was conducted in the NASA Langley Research Center 0.3-Meter Transonic Cryogenic Tunnel. In addition, a supplemental computational study on one of the inlet configurations was conducted using the Navier-Stokes flow solver, OVERFLOW. Tests were conducted at Mach numbers from 0.25 to 0.83, Reynolds numbers (based on aerodynamic interface plane diameter) from 5.1 million to 13.9 million (full-scale value), and inlet mass-flow ratios from 0.29 to 1.22, depending on Mach number. Results of the study indicated that increasing Mach number, increasing boundary layer thickness (relative to inlet height) or ingesting a boundary layer with a distorted profile decreased inlet performance. At Mach numbers above 0.4, increasing inlet airflow increased inlet pressure recovery but also increased distortion. Finally, inlet distortion was found to be relatively insensitive to Reynolds number, but pressure recovery increased slightly with increasing Reynolds number.

Berrier, Bobby L.↗

NASA NTRS2005

High Reynolds Number Investigation of a Flush Mounted, S-Duct Inlet With Large Amounts of Boundary Layer Ingestion

Berrier, Bobby L.↗

NASA NTRS2021

Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities such as line generation for general unstructured grids and optimization of the iterative solver. The efficiency of the new developments is assessed in a massively parallel environment. The assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Three benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere cylinder configuration, a transonic flow around the ONERA M6 wing, and a supersonic flow through a square duct are considered.

Mohagna J Pandya↗

NASA NTRSDate not supplied

Roughness Parameter Optimization of the McClain Model in GlennICE

Aircraft icing remains a significant threat to aviation safety. Software that predicts the impingement and ice accretion on full aircraft geometries and aircraft components are in demand and NASA Glenn is committed to produce software that meets this need. One of the key parameters affecting an accurate prediction of iced geometry is the effect of ice roughness on the heat transfer coefficient. While many efforts have been made to implement the roughness in the flow solver, this report takes a correlation for roughness height distribution that is based on experimental measurements and demonstrates how to relate those measurements to an augmentation to the heat transfer coefficient provided by the flow solution. The outcome of this effort was the calibration of defaults for user supplied parameters to this correlation through comparison with 95 large glaze conditions from experiment by adjusting user-supplied parameters in the roughness augmentation equation.

Aircraft Icing↗

DOE OSTI2021

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗

NASA NTRS1988

Three-dimensional incompressible Navier-Stokes computations of internal flows

Several incompressible Navier-Stokes solution methods for obtaining steady and unsteady solutions are discussed. Special attention is given to internal flows which involve distinctly different features from external flows. The characterisitcs of the flow solvers employing the method of pseudocompressibility and a fractional step method are briefly described. This discussion is limited to a primitive variable formulation in generalized curvilinear coordinates. Computed results include simple test cases and internal flow in the Space Shuttle main engine hot-gas manifold.

Kwak, D.↗

NASA NTRS1991

Three-dimensional incompressible Navier-Stokes solver using lower-upper symmetric-Gauss-Seidel algorithm

A numerical method based on the pseudocompressibility concept is developed for solving the three-dimensional incompressible Navier-Stokes equations using the lower-upper symmetric-Gauss-Seidel implicit scheme. Very high efficiency is achieved in a new flow solver, INS3D-LU code, by accomplishing the complete vectorizability of the algorithm on oblique planes of sweep in three dimensions.

Yoon, Seokkwan↗

NASA NTRS1993

Semi-discrete Galerkin solution of the compressible boundary-layer equations with viscous-inviscid interaction

A semi-discrete Galerkin (SDG) method is under development to model attached, turbulent, and compressible boundary layers for transonic airfoil analysis problems. For the boundary-layer formulation the method models the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby providing high resolution near the wall and permitting the use of a uniform finite element grid which automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past RAE 2822 and NACA 0012 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack.

Day, Brad A.↗

NASA NTRS2011

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗

NASA NTRS1999

A CFD Assessment of Several High-Lift Reference H Configuration Using Structured Grids

The objective of this study is to calibrate a Navier-Stokes code for a high-lift Reference H configuration using structured grids. The outline of this presentation will first include a brief description of the grids used and the flow solver. Next the results will be presented in terms of convergence and resources used on the C-90. Predicted force and moment and surface pressure results are compared to experiment and off- and on-surface flow viz. is discussed.

Lessard, Wendy B.↗

NASA NTRS2021

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗