Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow simulation”

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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Efficient Monte Carlo simulation of rarefied flow in a small nozzle

The flow of carbon dioxide through a small conical nozzle is computed using two different numerical approaches. The first of these is the direct simulation Monte Carlo method in which the physical phenomena are modeled at the molecular level. The second technique solves the Navier-Stokes equations of fluid flow. The purpose of the paper is to compare the solutions generated by these two methods for the conditions prevailing in a small, low density nozzle. The comparison between the two sets of solutions reveals only small differences. However, the flow properties in the thick viscous layer close to the nozzle at the exit plane show important deviations which have implications for the subsequent determination of the plume backflow.

Boyd, Iain D.↗
DOE OSTI2021

Preparing an Incompressible-Flow Fluid Dynamics Code for Exascale-Class Wind Energy Simulations

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices. The result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

algebraic multigrid↗
NASA NTRS1986

Computation of leading-edge vortex flows

The simulation of the leading edge vortex flow about a series of conical delta wings through solution of the Navier-Stokes and Euler equations is studied. The occurrence, the validity, and the usefulness of separated flow solutions to the Euler equations of particular interest. Central and upwind difference solutions to the governing equations are compared for a series of cross sectional shapes, including both rounded and sharp tip geometries. For the rounded leading edge and the flight condition considered, viscous solutions obtained with either central or upwind difference methods predict the classic structure of vortical flow over a highly swept delta wing. Predicted features include the primary vortex due to leading edge separation and the secondary vortex due to crossflow separation. Central difference solutions to the Euler equations show a marked sensitivity to grid refinement. On a coarse grid, the flow separates due to numerical error and a primary vortex which resembles that of the viscous solution is predicted. In contrast, the upwind difference solutions to the Euler equations predict attached flow even for first-order solutions on coarse grids. On a sufficiently fine grid, both methods agree closely and correctly predict a shock-curvature-induced inviscid separation near the leeward plane of symmetry. Upwind difference solutions to the Navier-Stokes and Euler equations are presented for two sharp leading edge geometries. The viscous solutions are quite similar to the rounded leading edge results with vortices of similar shape and size. The upwind Euler solutions predict attached flow with no separation for both geometries. However, with sufficient grid refinement near the tip or through the use of more accurate spatial differencing, leading edge separation results. Once the leading edge separation is established, the upwind solution agrees with recently published central difference solutions to the Euler equations.

Newsome, R. W.↗
NASA NTRS1985

Shock modelling in transonic and supersonic flow

Numerical simulation studies of problems in transonic flow are conducted in connection with the design of aircraft flying near the speed of sound, while recently there has also been a great deal of interest in supersonic and hypersonic flow calculations. Three important models for inviscid compressible gasdynamics are utilized in the study of these flow problems, taking into account the transonic small disturbance equation (TSD), the transonic full potential equation (FP), and the full Euler equations for compressible gasdynamics (EU). The present paper has the objective to provide a description of shock capturing methods recently developed for each model. The design principles behind the algorithms are also examined, giving attention to consistency, conservative form, monotone and sharp discrete shocks, and a bound on the variation, which prevents spurious oscillations.

Osher, S.↗
NASA NTRS1987

Numerical simulation of vortical flows over a strake-delta wing

The vortical flow fields over a double-delta wing with sweep angles 80 and 60 deg are investigated by means of numerical simulations based on the LU-ADI code of Fujii and Obayashi (1986) for the three-dimensional thin-layer Navier-Stokes equations. The results obtained using an 850,000-point grid at angle of attack alpha = 6-40 deg are presented graphically, compared with experimental data, and characterized in detail. Qualitatively accurate simulations are obtained for transition phenomena such as strake/wing-vortex interaction at alpha = 12 deg, bubble-type vortex breakdown near the trailing edge at alpha = 30 deg, and spiral-type breakdown at alpha = 35 deg. It is suggested that quantitative simulations can be achieved at higher grid resolution.

Fujii, Kozo↗
NASA NTRS1992

Simulator for concurrent processing data flow architectures

A software simulator capability of simulating execution of an algorithm graph on a given system under the Algorithm to Architecture Mapping Model (ATAMM) rules is presented. ATAMM is capable of modeling the execution of large-grained algorithms on distributed data flow architectures. Investigating the behavior and determining the performance of an ATAMM based system requires the aid of software tools. The ATAMM Simulator presented is capable of determining the performance of a system without having to build a hardware prototype. Case studies are performed on four algorithms to demonstrate the capabilities of the ATAMM Simulator. Simulated results are shown to be comparable to the experimental results of the Advanced Development Model System.

Malekpour, Mahyar R.↗
NASA NTRS1979

Numerical simulation of turbulent flows with a three-dimensional vortex-in-cell method

A three-dimensional vortex-in-cell method has been developed for the evaluation of local flow fields due to a family of vortex filaments which employs the principles and architecture of a code developed for magnetic field evaluation in plasma simulations. The computational effort in the new method, as compared to 'vortex pushing' by direct Biot-Savarat interaction, increases directly, rather than quadratically, with the number of vortex elements. The method is well suited for studying large number of vortex filaments or rings and can be used for simulating continuous vorticity.

Buneman, O.↗
DOE OSTI2018

Detailed Chemical Kinetics Based Simulation of Detonation-Containing Flows

Reliable and robust simulations of detonations in inhomogeneous and turbulent environments are of direct importance in the design of rotating detonation engines (RDEs). In particular, computational models will be especially useful in designing and optimizing discrete injectors that introduce fuel and air separately into the detonation chamber, but ensure appropriate level of mixing to sustain detonations but minimize backflow of detonation products and pressure waves into the feed plenums. Since the structure of detonations itself is non-ideal, models have to include a detailed description of this reacting flow in order to be predictive in nature. Here, a highly-scalable open source based solver has been developed for complex detonating flows such that a) the detonation processes are described using detailed chemical kinetics, b) the method is computationally efficient through the introduction of adaptive mesh refinement, and c) the solver can handle complex geometries of relevance to RDE design. Grid convergence of key metrics for detonations is evaluated using canonical flows. Further, the importance of the use of detailed chemical kinetics is illustrated by extracting the composition structure behind a two-dimensional detonation front. Finally, simulations of a practical RDE configuration are used to demonstrate the applicability of this solver to analyzing geometries. The simulation captures the general trends of the experiment well. It is found that the detonation occurs under partially-premixed conditions. Propagation of pressure waves to the injection system is observed which could influence flow behavior in the oxidizer plenum.

03 NATURAL GAS↗
DOE OSTI2014

Large-eddy simulation of turbulent cavitating flow in a micro channel

Large-eddy simulations (LES) of cavitating flow of a Diesel-fuel-like fluid in a generic throttle geometry are presented. Two-phase regions are modeled by a parameter-free thermodynamic equilibrium mixture model, and compressibility of the liquid and the liquid-vapor mixture is taken into account. The Adaptive Local Deconvolution Method (ALDM), adapted for cavitating flows, is employed for discretizing the convective terms of the Navier-Stokes equations for the homogeneous mixture. ALDM is a finite-volume-based implicit LES approach that merges physically motivated turbulence modeling and numerical discretization. Validation of the numerical method is performed for a cavitating turbulent mixing layer. Comparisons with experimental data of the throttle flow at two different operating conditions are presented. The LES with the employed cavitation modeling predicts relevant flow and cavitation features accurately within the uncertainty range of the experiment. The turbulence structure of the flow is further analyzed with an emphasis on the interaction between cavitation and coherent motion, and on the statistically averaged-flow evolution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

2D numerical simulation of unsteady flows for large scale floods prediction in real time

The challenge of finding a compromise between computational time and level of accuracy and robustness has traditionally expanded the use simplified models rather than full two-dimensional (2D) models for flood simulation. This work presents a GPU accelerated 2D shallow water model for the simulation of flood events in real time. In particular, an explicit first-order finite volume scheme is detailed to control the numerical instabilities that are likely to appear when used in complex topography. The model is first validated with the benchmark test case of the Toce River (Italy) and numerical fixes are demonstrated to be necessary. The model is next applied to reproduce real events in a reach of the Ebro River (Spain) in order to compare simulation results with field data. The second case deals with a large domain (744 km 2 ) and long flood duration (up to 20 days) allowing an analysis of the performance and speed-up achieved by different GPU devices. We report the high values of fit between observed and simulated results as well as the computational times achieved are encouraging to propose the use of the model as forecasting system.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Dynamic localization and second-order subgrid-scale models in large eddy simulations of channel flow

The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Cabot, William H.↗
NASA NTRS1989

Simulation of turbomachinery flows

Significant advancements have been made in the last five years in the ability to model turbomachinery flows of engineering interest. This advancement can be directly attributed to the second generation of supercomputers like the Cray XMP and Cray 2 and advanced instrumentation techniques. Early on, the National Aeronautics and Space Administration Lewis Research Center recognized the potential gains in turbomachinery performance and life that could be achieved by taking advantage of this technology and instituted a comprehensive research program in turbomachinery flow modeling. This activity combined the areas of fluid flow analysis, computational fluid dynamics, and experimental fluid mechanics. As a result of this activity, Lewis has become an internationally recognized leader in turbomachinery flow modeling. Many of the research activities conducted under this program are utilized by industry. The presentation gives an overview of this program and provides sample illustration of simulation performed to date.

Adamczyk, John J.↗
NASA NTRS2001

Time-Dependent Simulations of Turbopump Flows

The objective of the current effort is to provide a computational framework for design and analysis of the entire fuel supply system of a liquid rocket engine, including high-fidelity unsteady turbopump flow analysis. This capability is needed to support the design of pump sub-systems for advanced space transportation vehicles that are likely to involve liquid propulsion systems. To date, computational tools for design/analysis of turbopump flows are based on relatively lower fidelity methods. An unsteady, three-dimensional viscous flow analysis tool involving stationary and rotational components for the entire turbopump assembly has not been available for real-world engineering applications. The present effort will provide developers with information such as transient flow phenomena at start up, impact of non-uniform inflows, system vibration and impact on the structure. In the proposed paper, the progress toward the capability of complete simulation of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of the hybrid MPI/Open-MP and MLP versions of the INS3D code. The relative motion of the grid systems for the rotor-stator interaction was obtained using overset grid techniques. Time-accuracy of the scheme has been evaluated with simple test cases. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on Origin 2000 systems at NASA Ames Research Center. Results from these time-accurate simulations with moving boundary capability will be presented along with the performance of parallel versions of the code.

Kris, Cetin C.↗
NASA NTRS1991

Towards large eddy and direct simulation of complex turbulent flows

Recent advances in the methodology for direct numerical simulation of turbulent flows and some of the current applications are reviewed. It is argued that high-order finite difference schemes yield solutions with comparable accuracy to the spectral methods with the same number of degrees of freedom. The effects of random inflow conditions on the downstream evolution of turbulence are discussed.

Moin, Parviz↗
NASA NTRS1992

Knowledge representation requirements for model sharing between model-based reasoning and simulation in process flow domains

The paper examines the requirements for the reuse of computational models employed in model-based reasoning (MBR) to support automated inference about mechanisms. Areas in which the theory of MBR is not yet completely adequate for using the information that simulations can yield are identified, and recent work in these areas is reviewed. It is argued that using MBR along with simulations forces the use of specific fault models. Fault models are used so that a particular fault can be instantiated into the model and run. This in turn implies that the component specification language needs to be capable of encoding any fault that might need to be sensed or diagnosed. It also means that the simulation code must anticipate all these faults at the component level.

Throop, David R.↗