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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Turbulence in unsteady flow at high frequencies

Turbulent flows subjected to oscillations of the mean flow were simulated using a large-eddy simulation computer code for flow in a channel. The objective of the simulations was to provide better understanding of the effects of time-dependent disturbances on the turbulence of a boundary layer and of the underlying physical phenomena regarding the basic interaction between the turbulence and external disturbances. The results confirmed that turbulence is sensitive to certain ranges of frequencies of disturbances. However, no direct connection was found between the frequency of imposed disturbances and the characteristic 'burst' frequency of turbulence. New insight into the nature of turbulence at high frequencies was found. Viscous phenomena near solid walls were found to be the dominant influence for high-frequency perturbations.

Kuhn, Gary D.↗
DOE OSTI2023

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗
DOE OSTI2017

Direct Numerical Simulation of Pebble Bed Flows: Database Development and Investigation of Low-Frequency Temporal Instabilities

Computational analyses of fluid flow through packed pebble bed domains using the Reynolds-averaged NavierStokes framework have had limited success in the past. Because of a lack of high-fidelity experimental or computational data, optimization of Reynolds-averaged closure models for these geometries has not been extensively developed. In the present study, direct numerical simulation was employed to develop a high-fidelity database that can be used for optimizing Reynolds-averaged closure models for pebble bed flows. A face-centered cubic domain with periodic boundaries was used. Flow was simulated at a Reynolds number of 9308 and cross-verified by using available quasi-DNS data. During the simulations, low-frequency instability modes were observed that affected the stationary solution. Furthermore, these instabilities were investigated by using the method of proper orthogonal decomposition, and a correlation was found between the time-dependent asymmetry of the averaged velocity profile data and the behavior of the highest energy eigenmodes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

A Mixed Finite Volume Element Method for Flow Calculations in Porous Media

A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Jim E Jones↗
NASA NTRS2006

PAB3D: Its History in the Use of Turbulence Models in the Simulation of Jet and Nozzle Flows

This is a review paper for PAB3D s history in the implementation of turbulence models for simulating jet and nozzle flows. We describe different turbulence models used in the simulation of subsonic and supersonic jet and nozzle flows. The time-averaged simulations use modified linear or nonlinear two-equation models to account for supersonic flow as well as high temperature mixing. Two multiscale-type turbulence models are used for unsteady flow simulations. These models require modifications to the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k-epsilon) model with a RANS/LES transition function, dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) formulation. All of these models are implemented in the three-dimensional Navier-Stokes code PAB3D. This paper discusses computational methods, code implementation, computed results for a wide range of nozzle configurations at various operating conditions, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions.

Abdol-Hamid, Khaled S.↗
NASA NTRS1996

Study of the velocity gradient tensor in turbulent flow

The behavior of the velocity gradient tensor, A(ij)=delta u(i)/delta x(j), was studied using three turbulent flows obtained from direct numerical simulation The flows studies were: an inviscid calculation of the interaction between two vortex tubes, a homogeneous isotropic flow, and a temporally evolving planar wake. Self-similar behavior for each flow was obtained when A(ij) was normalized with the mean strain rate. The case of the interaction between two vortex tubes revealed a finite sized coherent structure with topological characteristics predictable by a restricted Euler model. This structure was found to evolve with the peak vorticity as the flow approached singularity. Invariants of A(ij) within this structure followed a straight line relationship of the form: gamma(sup 3)+gammaQ+R=0, where Q and R are the second and third invariants of A(ij), and the eigenvalue gamma is nearly constant over the volume of this structure. Data within this structure have local strain topology of unstable-node/saddle/saddle. The characteristics of the velocity gradient tensor and the anisotropic part of a related acceleration gradient tensor H(ij) were also studied for a homogeneous isotropic flow and a temporally evolving planar wake. It was found that the intermediate principal eigenvalue of the rate-of-strain tensor of H(ij) tended to be negative, with local strain topology of the type stable-node/saddle/saddle. There was also a preferential eigenvalue direction. The magnitude of H(ij) in the wake flow was found to be very small when data were conditioned at high local dissipation regions. This result was not observed in the relatively low Reynolds number simulation of homogeneous isotropic flow. A restricted Euler model of the evolution of A(ij) was found to reproduce many of the topological features identified in the simulations.

Cheng, Wei-Ping↗
DOE OSTI2022

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗
DOE OSTI2020

Improved pore network models to simulate single-phase flow in porous media by coupling with lattice Boltzmann method

In this paper, different pore network models to simulate single-phase flow in porous media are built and their accuracy is evaluated. In addition to the conventional pore network model (CPNM) which consists of regular pore bodies and throat bonds, three improved pore network models (IPNMs) are developed allowing to better describing the real pore and throat geometry. The first improved pore network model (IPNM1) replaces the regular throat bond with a throat bond showing the real throat cross section. The second improvement (IPNM2) uses a series of sub-throat bonds with varying cross sections to better describe the real throat geometry, which is firstly proposed in this paper. The third model (IPNM3) extracts the real pore-throat-pore geometry without simplification. The conductance of fluid flow through these more realistic throat bonds is calculated by the lattice Boltzmann method (LBM). The accuracy and computational efficiency of the different pore network models are evaluated taking the LBM simulation over the whole porous medium as reference solution. The global permeability and detailed pressure distributions in the pores for the different pore network models are validated. The results show that the accuracy of the pore network model increases from CPNM to IPNM3, but at the expense of increasing computational cost. This study suggests that IPNM3 can replace a whole-domain LBM simulation with similar accuracy but much lower computational cost. As a first-order approximation the newly proposed IPNM2 is suggested as good compromise between accuracy and computational cost.

42 ENGINEERING↗
NASA NTRS1997

Eulerian Time-Domain Filtering for Spatial LES

Eulerian time-domain filtering seems to be appropriate for LES (large eddy simulation) of flows whose large coherent structures convect approximately at a common characteristic velocity; e.g., mixing layers, jets, and wakes. For these flows, we develop an approach to LES based on an explicit second-order digital Butterworth filter, which is applied in,the time domain in an Eulerian context. The approach is validated through a priori and a posteriori analyses of the simulated flow of a heated, subsonic, axisymmetric jet.

Pruett, C. David↗
NASA NTRS1985

The investigation of inlet/nozzle flowfield coupling using compact propulsion simulators

An investigation of airframe/propulsion system flowfield coupling has been conducted in the Ames 11-Foot Transonic Wind Tunnel at Mach numbers from 0.40 to 1.40 and at angles-of-attack from -2 to +20 deg. Propulsive flows were simulated by flow-through, jet-effects, and turbo-powered simulators (CMAPS) techniques. The configuration was a 10 percent scale model of a VSTOL canard/wing fighter configuration with twin podded engine nacelles and nonaxisymmetric nozzles. Localized flowfield interactions were present and inlet-nozzle flowfield coupling was identified at supersonic flight conditions. The configuration with simultaneous inlet and nozzle flow simulation (CMAPS) exhibited lower drag supersonically than the flow-through/jet-effects configuration. During this first wind tunnel application of twin simulators in a full aircraft configuration, the CMAPS and control system proved to be a flexible, reliable testing method.

Zilz, D. E.↗
NASA NTRS1986

Numerical simulations of rarefied reentry flows

A review of recent calculations obtained with Bird's direct simulation Monte Carlo (DSMC) method is presented for the transitional flowfield encountered at reentry conditions. Consequently, the emphasis is the real-gas effects resulting from a highly energetic nonequilibrium flow. The DSMC calculations for both wide-angle and slender bodies are compared with continuum calculations. The wide-angle body calculations simulate the windward nose of the Shuttle Orbiter and projected aeroassisted orbital transfer vehicles. For the slender, blunted configurations, both two-dimensional and axisymmetric results are presented. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant.

Moss, James N.↗
NASA NTRS1990

On the application of large-eddy simulation to transitional wall-bounded flows

Large-eddy simulations of transition on a flat-plate boundary layer have been carried out. Standard residual stress models of the Smagorinsky type are excessively dissipative, but rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, leading to results that compare favorably with those of a direct simulation, at a small fraction of the computational effort.

Piomelli, Ugo↗
DOE OSTI2018

Robust system size reduction of discrete fracture networks: a multi-fidelity method that preserves transport characteristics

In this paper, we propose a multi-fidelity system reduction technique that uses weighted graphs paired with three-dimensional discrete fracture network (DFN) modelling for efficient simulation of subsurface flow and transport in fractured media. DFN models are used to simulate flow and transport in subsurface fractured rock with low-permeability. One method to alleviate the heavy computational overhead associated with these simulations is to reduce the size of the DFN using a graph representation of it to identify the primary flow sub-network and only simulate flow and transport thereon. The first of these methods used unweighted graphs constructed solely on DFN topology and could be used for accurate predictions of first-passage times. However, these techniques perform poorly when predicting later stages of the mass breakthrough. We utilize a weighted-graph representation of the DFN where edge weights are based on hydrological parameters in the DFN that allows us to exploit the kinematic quantities derivable a posteriori from the flow solution obtained on the graph representation of the DFN to perform system reduction and predict the later stages of the breakthrough curve with high fidelity. We also propose and demonstrate the use of an adaptive pruning algorithm with error control that produces a pruned DFN sub-network whose predicted mass breakthrough agrees with the original DFN within a user-specified tolerance. Lastly, the method allows for the level of accuracy to be a user-controlled parameter.

58 GEOSCIENCES↗
NASA NTRS2012

Numerical Investigation of PLIF Gas Seeding for Hypersonic Boundary Layer Flows

Numerical simulations of gas-seeding strategies required for planar laser-induced fluorescence (PLIF) in a Mach 10 air flow were performed. The work was performed to understand and quantify adverse effects associated with gas seeding and to compare different flow rates and different types of seed gas. The gas was injected through a slot near the leading edge of a flat plate wedge model used in NASA Langley Research Center's 31- Inch Mach 10 Air Tunnel facility. Nitric oxide, krypton, and iodine gases were simulated at various injection rates. Simulation results showing the deflection of the velocity field for each of the cases are presented. Streamwise distributions of velocity and concentration boundary layer thicknesses as well as vertical distributions of velocity, temperature, and mass distributions are presented for each of the cases. Relative merits of the different seeding strategies are discussed.

Johanson, Craig T.↗
NASA NTRS1993

Time accurate simulations of compressible shear flows

The objectives of this research are to employ direct numerical simulation (DNS) to study the phenomenon of mixing (or lack thereof) in compressible free shear flows and to suggest new means of enhancing mixing in such flows. The shear flow configurations under investigation are those of parallel mixing layers and planar jets under both non-reacting and reacting nonpremixed conditions. During the three-years of this research program, several important issues regarding mixing and chemical reactions in compressible shear flows were investigated.

Givi, Peyman↗
DOE OSTI2019

Comprehensive study of the interactions between the critical dimensionless numbers associated with multiphase flow in 3D porous media

Multiphase flow in porous media is of signifcant interest in many engineering applications, such as geologic carbon sequestration, enhanced oil recovery, and groundwater contamination and remediation. In order to advance the fundamental understanding of multiphase flow in complex three-dimensional (3D) porous media, the interactions between the critical dimensionless numbers, including the contact angle, viscosity ratio, and capillary (Ca) number, were explored using X-ray micro-computed tomography (micro-CT) scanning and lattice Boltzmann (LB) modeling. In this study, the 3D pore structure information was extracted from micro-CT images and then used as interior boundary conditions of flow modeling in a pore-scale LB simulator to simulate multiphase flow within the pore space. A Berea sandstone sample was scanned and then two-phase flow LB simulations were performed based on the micro-CT images. The LB-simulated water/CO 2 distributions agreed well with the micro-CT scanned images. Simulation findings showed that a decreasing contact angle causes a decrease in wetting-fluid relative permeability and an increase in non-wetting fluid relative permeability. A rising Ca number increases both wetting and non-wetting fluid relative permeabilities. An increasing viscosity ratio (the ratio of non-wetting fluid viscosity to wetting fluid viscosity) facilitates the increase of non-wetting fluid relative permeability and mitigates the reduction of wetting fluid relative permeability, when the contact angle decreases continuously. The primary novel finding of this study is that the viscosity ratio affects the rate of change of the relative permeability curves for both phases when the contact angle changes continuously. To date, we believe it is the first time that comprehensive interactions between these dimensionless numbers are demonstrated in a sandstone sample based on real 3D structures. We also investigated the role of the changes of flow direction and sample location on relative permeability curves. Simulation results showed that the change in non-wetting fluid relative permeability was larger when the flow direction was switched from vertical to horizontal, which indicated that there was stronger anisotropy in larger pore networks that were primarily occupied by the non-wetting fluid.

42 ENGINEERING↗
NASA NTRS1988

Numerical simulation of hypersonic inlet flows with equilibrium or finite rate chemistry

An efficient numerical program incorporated with comprehensive high temperature gas property models has been developed to simulate hypersonic inlet flows. The computer program employs an implicit lower-upper time marching scheme to solve the two-dimensional Navier-Stokes equations with variable thermodynamic and transport properties. Both finite-rate and local-equilibrium approaches are adopted in the chemical reaction model for dissociation and ionization of the inlet air. In the finite rate approach, eleven species equations coupled with fluid dynamic equations are solved simultaneously. In the local-equilibrium approach, instead of solving species equations, an efficient chemical equilibrium package has been developed and incorporated into the flow code to obtain chemical compositions directly. Gas properties for the reaction products species are calculated by methods of statistical mechanics and fit to a polynomial form for C(p). In the present study, since the chemical reaction time is comparable to the flow residence time, the local-equilibrium model underpredicts the temperature in the shock layer. Significant differences of predicted chemical compositions in shock layer between finite rate and local-equilibrium approaches have been observed.

Yu, Sheng-Tao↗
NASA NTRS2006

Blunt Body Aerodynamics for Hypersonic Low Density Flows

Numerical simulations are performed for the Apollo capsule from the hypersonic rarefied to the continuum regimes. The focus is on flow conditions similar to those experienced by the Apollo 6 Command Module during the high altitude portion of its reentry. The present focus is to highlight some of the current activities that serve as a precursor for computational tool assessments that will be used to support the development of aerodynamic data bases for future capsule flight environments, particularly those for the Crew Exploration Vehicle (CEV). Results for aerodynamic forces and moments are presented that demonstrate their sensitivity to rarefaction; that is, free molecular to continuum conditions. Also, aerodynamic data are presented that shows their sensitivity to a range of reentry velocities, encompassing conditions that include reentry from low Earth orbit, lunar return, and Mars return velocities (7.7 to 15 km/s). The rarefied results obtained with direct simulation Monte Carlo (DSMC) codes are anchored in the continuum regime with data from Navier-Stokes simulations.

Moss, James N.↗