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 73 records · Page 4

A perspective on unstructured grid flow solvers

This survey paper assesses the status of compressible Euler and Navier-Stokes solvers on unstructured grids. Different spatial and temporal discretization options for steady and unsteady flows are discussed. The integration of these components into an overall framework to solve practical problems is addressed. Issues such as grid adaptation, higher order methods, hybrid discretizations and parallel computing are briefly discussed. Finally, some outstanding issues and future research directions are presented.

Venkatakrishnan, V.↗

Flow Solver for Incompressible 2-D Drive Cavity

This software solves the Navier-Stokes equations for the incompressible driven cavity flow problem. The code uses second-order finite differencing on a staggered grid using the Chorin projection method. The resulting intermediate Poisson equation is efficiently solved using the fast Fourier transform. Time stepping is done using fourth-order Runge-Kutta for stability at high Reynolds numbers. Features include check-pointing, periodic field snapshots, ongoing reporting of kinetic energy and changes between time steps, time histories at selected points, and optional streakline generation.

Kalb, Virginia L.↗

Numerical Predictions of Dust-Induced Heat Flux Augmentation in Hypersonic Blunt-Body Flows Using a Discontinuous Galerkin Multiphase Flow Solver

Recent interest in human-scale missions to Mars has motivated the need for high-fidelity simulations of reentry flows. During a dust storm, there can be high levels of suspended dust in the Martian atmosphere, which cannot only enhance erosion of thermal protection systems but also transfer energy and momentum to the shock layer, thereby significantly augmenting the surface heat flux. Second-order finite-volume schemes are typically employed for hypersonic flow simulations, but such schemes suffer from a number of disadvantages. An attractive alternative is discontinuous Galerkin methods, which benefit from arbitrarily high spatial order of accuracy, geometric flexibility, and other properties. To enable accurate computations of high-speed particle-laden flows, an Euler-Lagrange methodology was developed in which the Eulerian field of the carrier gas is calculated using a discontinuous Galerkin scheme while the disperse phase is treated with Lagrangian particle tracking. We discuss challenges associated with coupling these two formulations and how to handle them. Momentum and energy transfer between the carrier gas and the particle phase is considered, and the importance of accounting for interparticle collisions is assessed. In addition, we describe the physical model of the particle phase and examine effects of its uncertainties on the numerical solution. We demonstrate the performance of the Euler-Lagrange method in representative testcases, with focus on the accurate prediction of particle trajectories and heating augmentation. Quantitative comparisons with experiments are provided.

Ching, Eric J.↗

Development of a flexible and efficient multigrid-based multiblock flow solver

A multigrid-based numerical scheme for solving the Navier-Stokes and Euler equations on block-structured grids is presented for flows over complex aircraft configurations. Self-consistency checks are conducted to ascertain that the multiblock code reproduces previously known solutions from single-block codes by splitting a single-block grid into multiple blocks. The newly-developed multiblock code maintains multigrid performance of an equivalent single-block code.

Vatsa, Veer N.↗

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗

A Flow Solver for Three-Dimensional DRAGON Grids

DRAGONFLOW code has been developed to solve three-dimensional Navier-Stokes equations over a complex geometry whose flow domain is discretized with the DRAGON grid-a combination of Chimera grid and a collection of unstructured grids. In the DRAGONFLOW suite, both OVERFLOW and USM3D are presented in form of module libraries, and a master module controls the invoking of these individual modules. This report includes essential aspects, programming structures, benchmark tests and numerical simulations.

Liou, Meng-Sing↗

MEDUSA - An overset grid flow solver for network-based parallel computer systems

Continuing improvement in processing speed has made it feasible to solve the Reynolds-Averaged Navier-Stokes equations for simple three-dimensional flows on advanced workstations. Combining multiple workstations into a network-based heterogeneous parallel computer allows the application of programming principles learned on MIMD (Multiple Instruction Multiple Data) distributed memory parallel computers to the solution of larger problems. An overset-grid flow solution code has been developed which uses a cluster of workstations as a network-based parallel computer. Inter-process communication is provided by the Parallel Virtual Machine (PVM) software. Solution speed equivalent to one-third of a Cray-YMP processor has been achieved from a cluster of nine commonly used engineering workstation processors. Load imbalance and communication overhead are the principal impediments to parallel efficiency in this application.

Smith, Merritt H.↗

Status Of The UPS Space-Marching Flow Solver

The status of the three-dimensional parabolized Navier-Stokes solver UPS is described. The UPS code, initiated at NASA Ames Research Center in 1986, continues to develop and evolve through application to supersonic and hypersonic flow fields. Hypersonic applications have motivated enhancement of the physical modeling capabilities of the code, specifically real gas modeling, boundary conditions, and turbulence and transition modeling. The UPS code has also been modified to enhance robustness and efficiency in order to be practically used in concert with an optimization code for supersonic transport design. These developments are briefly described along with some relevant results for generic test problems obtained during verification of the enhancements. Included developments and results have previously been published and widely disseminated domestically.

Lawerence, Scott L.↗

Adaptive unstructured triangular mesh generation and flow solvers for the Navier-Stokes equations at high Reynolds number

A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.

Ashford, Gregory A.↗

General Flow-Solver Code for Turbomachinery Applications

Phantom is a computer code intended primarily for real-fluid turbomachinery problems. It is based on Corsair, an ideal-gas turbomachinery code, developed by the same authors, which evolved from the ROTOR codes from NASA Ames. Phantom is applicable to real and ideal fluids, both compressible and incompressible, flowing at subsonic, transonic, and supersonic speeds. It utilizes structured, overset, O- and H-type zonal grids to discretize flow fields and represent relative motions of components. Values on grid boundaries are updated at each time step by bilinear interpolation from adjacent grids. Inviscid fluxes are calculated to third-order spatial accuracy using Roe s scheme. Viscous fluxes are calculated using second-order-accurate central differences. The code is second-order accurate in time. Turbulence is represented by a modified Baldwin-Lomax algebraic model. The code offers two options for determining properties of fluids: One is based on equations of state, thermodynamic departure functions, and corresponding state principles. The other, which is more efficient, is based on splines generated from tables of properties of real fluids. Phantom currently contains fluid-property routines for water, hydrogen, oxygen, nitrogen, kerosene, methane, and carbon monoxide as well as ideal gases.

Dorney, Daniel↗

Applying and validating the RANS-3D flow-solver for evaluating a subsonic serpentine diffuser geometry

Subsonic inlet ducts for advanced, high-performance aircraft are evolving towards complex three-dimensional shapes for reasons of overall integration and weight. These factors lead to diffuser geometries that may sacrifice inlet performance, unless careful attention to design details and boundary layer management techniques are employed. The ability of viscous computational fluid dynamic (CFD) analysis of such geometries to aid the aircraft configurator in this complex design problem is herein examined. The RANS-3D Reynolds-Averaged Navier-Stokes solver is applied to model the complex flowfield occurring in a representative diffuser geometry and the solutions are compared to experimental results from a static test of the inlet duct. The computational results are shown to compare very favorably with experimental results over a range of mass flow rates, including those involving large amounts of separation in the diffuser. In addition, a novel grid topology is presented, and two turbulence models are evaluated in this study as part of the RANS-3D code.

Fletcher, Michael J.↗

Mesoflow: An Open-Source Reacting Flow Solver for Catalysis at Mesoscale

We present the capabilities and software performance metrics of our open-source continuum solver for catalysis, Mesoflow, developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex catalyst surface morphologies directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous catalyst interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPUs) compared to a single processor for representative problem sizes (2 million cell mesh). We will also present a brief introduction on how to build and use this software for application problems pertaining to catalytic upgrading and gas transport within porous catalyst particles.

adaptive meshing↗

Assessment of Various Flow Solvers Used to Predict the Thermal Environment inside Space Shuttle Solid Rocket Motor Joints

It is very important to accurately predict the gas pressure, gas and solid temperature, as well as the amount of O-ring erosion inside the space shuttle Reusable Solid Rocket Motor (RSRM) joints in the event of a leak path. The scenarios considered are typically hot combustion gas rapid pressurization events of small volumes through narrow and restricted flow paths. The ideal method for this prediction is a transient three-dimensional computational fluid dynamics (CFD) simulation with a computational domain including both combustion gas and surrounding solid regions. However, this has not yet been demonstrated to be economical for this application due to the enormous amount of CPU time and memory resulting from the relatively long fill time as well as the large pressure and temperature rising rate. Consequently, all CFD applications in RSRM joints so far are steady-state simulations with solid regions being excluded from the computational domain by assuming either a constant wall temperature or no heat transfer between the hot combustion gas and cool solid walls.

Wang, Qun-Zhen↗

An incompressible Navier-Stokes flow solver in three-dimensional curvilinear coordinate systems using primitive variables

An implicit, finite-difference computer code has been developed to solve the incompressible Navier-Stokes equations in a three-dimensional, curvilinear coordinate system. The pressure-field solution is based on the pseudo compressibility approach in which the time derivative pressure term is introduced into the mass conservation equation to form a set of hyperbolic equations. The solution procedure employs an implicit, approximate factorization scheme. The Reynolds stresses, that are uncoupled from the implicit scheme, are lagged by one time-step to facilitate implementing various levels of the turbulence model. Test problems for external and internal flows are computed, and the results are compared with existing experimental data. The application of this technique for general three-dimensional problems is then demonstrated.

Kwak, D.↗

A solution adaptive grid procedure for an upwind parabolized flow solver

A solution adaptive grid procedure based on an error equi-distribution scheme is developed and applied to a Parabolized Navier-Stokes solver. An improved method for selecting weighting functions is introduced which involves normalizing a combination of flowfield gradients and curvature of a number of dependent variables and then selecting the largest at each point. The scheme re-distributes grid points line-by-line, with grid point motion controlled by forces analogous to tensional and torsional spring forces with the spring constants set equal to the weighting functions. Torsional terms are functions of the grid point positions along neighboring grid lines and provide grid smoothness and stability. A grid-fitting scheme is introduced for external flows in which the number of grid points in the freestream are reduced to a minimum. Results for several problems are presented to demonstrate the improvements obtainable with the solution adaptive grid procedure.

Harvey, Albert D.↗

Implementation of Implicit Adaptive Mesh Refinement in an Unstructured Finite-Volume Flow Solver

This paper explores the implementation of adaptive mesh refinement in an unstructured, finite-volume solver. Unsteady and steady problems are considered. The effect on the recovery of high-order numerics is explored and the results are favorable. Important to this work is the ability to provide a path for efficient, implicit time advancement. A method using a simple refinement sensor based on undivided differences is discussed and applied to a practical problem: a shock-shock interaction on a hypersonic, inviscid double-wedge. Cases are compared to uniform grids without the use of adapted meshes in order to assess error and computational expense. Discussion of difficulties, advances, and future work prepare this method for additional research. The potential for this method in more complicated flows is described.

Schwing, Alan M.↗