Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulent Flow Simulations”

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 91 records · Page 5

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗

Verification of Anisotropic Mesh Adaptation for Complex Aerospace Applications

The stabilized finite element solver, FUN3D-SFE, along with the grid mechanics package refine are verified for aerospace applications of laminar and turbulent flow simulations. The current verification exercise represents an extension of previous research using FUN3D-SFEwith adjoint-based mesh adaptation to generate highly anisotropic adapted meshes for inviscid problems. Adaptations are performed using a solution-based approach that controls the Lpnorm of Mach number interpolation error and an adjoint-based approach that controls the error in some output functional. Adaptive results are shown for laminar subsonic flow over a delta wing, laminar subsonic flow over ONERA M6 wing, inviscid supersonic flow over a sonic boom test case, and a turbulent flow over a high lift configuration (JAXA Standard Model). Mesh convergence results are also compared with results from FUN3D-FV (Finite Volume) whenever available. For all test cases considered, FUN3D-SFE gives significantly better accurate results than FUN3D-FV on coarse meshes.

Aravind Balan↗

Numerical method for predicting flow characteristics and performance of nonaxisymmetric nozzles. Part 2: Applications

A computer implemented numerical method for predicting the flow in and about an isolated three dimensional jet exhaust nozzle is summarized. The approach is based on an implicit numerical method to solve the unsteady Navier-Stokes equations in a boundary conforming curvilinear coordinate system. Recent improvements to the original numerical algorithm are summarized. Equations are given for evaluating nozzle thrust and discharge coefficient in terms of computed flowfield data. The final formulation of models that are used to simulate flow turbulence effect is presented. Results are presented from numerical experiments to explore the effect of various quantities on the rate of convergence to steady state and on the final flowfield solution. Detailed flowfield predictions for several two and three dimensional nozzle configurations are presented and compared with wind tunnel experimental data.

Thomas, P. D.↗

Progress and supercomputing in computational fluid dynamics; Proceedings of U.S.-Israel Workshop, Jerusalem, Israel, December 1984

Current developments and future trends in the application of supercomputers to computational fluid dynamics are discussed in reviews and reports. Topics examined include algorithm development for personal-size supercomputers, a multiblock three-dimensional Euler code for out-of-core and multiprocessor calculations, simulation of compressible inviscid and viscous flow, high-resolution solutions of the Euler equations for vortex flows, algorithms for the Navier-Stokes equations, and viscous-flow simulation by FEM and related techniques. Consideration is given to marching iterative methods for the parabolized and thin-layer Navier-Stokes equations, multigrid solutions to quasi-elliptic schemes, secondary instability of free shear flows, simulation of turbulent flow, and problems connected with weather prediction.

Murman, E. M.↗

Studying marine stratus with large eddy simulation

Data sets from field experiments over the stratocumulus regime may include complications from larger scale variations, decoupled cloud layers, diurnal cycle, or entrainment instability, etc. On top of the already complicated turbulence-radiation-condensation processes within the cloud-topped boundary layer (CTBL), these complexities may sometimes make interpretation of the data sets difficult. To study these processes, a better understanding is needed of the basic processes involved in the prototype CTBL. For example, is cloud top radiative cooling the primary source of the turbulent kinetic energy (TKE) within the CTBL. Historically, laboratory measurements have played an important role in addressing the turbulence problems. The CTBL is a turbulent field which is probably impossible to generate in laboratories. Large eddy simulation (LES) is an alternative way of 'measuring' the turbulent structure under controlled environments, which allows the systematic examination of the basic physical processes involved. However, there are problems with the LES approach for the CTBL. The LES data need to be consistent with the observed data. The LES approach is discussed, and results are given which provide some insights into the simulated turbulent flow field. Problems with this approach for the CTBL and information from the FIRE experiment needed to justify the LES results are discussed.

Moeng, Chin-Hoh↗

Turbulent transport in the solar nebula

This paper describes the current state of an ongoing project to simulate turbulent flow in a solar nebula, which is the flattened disk of dust and gas out of which a solar system forms. The goal of this project is to determine a model for the transport of mass and angular momentum in the nebula. The nebula flow exhibits compressibility, thermal conduction, viscosity, internal heating through viscous dissipation, a stable shear due to Keplerian rotation, and a gravitational acceleration in the vertical direction which is linear with altitude. These properties combine to give flow patterns not seen in terrestrial applications. Primordial solar systems are known to exist and are presumably undergoing an evolution similar to the early stages of our own solar system; for example, the IRAS infrared telescope has discovered such a protoplanetary system around the star Vega. Solar nebula evolution is the subject of much research in the astrophysical community. In the long run, researchers hope to gain a better understanding of planetary formation and the processes which dissipate the solar nebula with time.

Thompson, K. W.↗

NPARC Code Upgraded with Two-Equation Turbulence Models

The National PARC (NPARC) Alliance was established by the NASA Lewis Research Center and the Air Force Arnold Engineering Development Center to provide the U.S. aeropropulsion community with a reliable Navier-Stokes code for simulating the nonrotating components of propulsion systems. Recent improvements to the turbulence model capabilities of the NPARC code have significantly improved its capability to simulate turbulent flows. Specifically, the Chien k-epsilon and Wilcox k-omega turbulence models were implemented at Lewis. Lewis researchers installed the Chien k-epsilon model into NPARC to improve the code's ability to calculate turbulent flows with attached wall boundary layers and free shear layers. Calculations with NPARC have demonstrated that the Chien k-epsilon model provides more accurate calculations than those obtained with algebraic models previously available in the code. Grid sensitivity investigations have shown that computational grids must be packed against the solid walls such that the first point off of the wall is placed in the laminar sublayer. In addition, matching the boundary layer and momentum thicknesses entering mixing regions is necessary for an accurate prediction of the free shear-layer growth.

Source record↗

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes - Part 1

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling.

Diskin, Boris↗

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes

Cell-centered and node-centered approaches have been compared for unstructured finite-volume discretization of inviscid fluxes. The grids range from regular grids to irregular grids, including mixed-element grids and grids with random perturbations of nodes. Accuracy, complexity, and convergence rates of defect-correction iterations are studied for eight nominally second-order accurate schemes: two node-centered schemes with weighted and unweighted least-squares (LSQ) methods for gradient reconstruction and six cell-centered schemes two node-averaging with and without clipping and four schemes that employ different stencils for LSQ gradient reconstruction. The cell-centered nearest-neighbor (CC-NN) scheme has the lowest complexity; a version of the scheme that involves smart augmentation of the LSQ stencil (CC-SA) has only marginal complexity increase. All other schemes have larger complexity; complexity of node-centered (NC) schemes are somewhat lower than complexity of cell-centered node-averaging (CC-NA) and full-augmentation (CC-FA) schemes. On highly anisotropic grids typical of those encountered in grid adaptation, discretization errors of five of the six cell-centered schemes converge with second order on all tested grids; the CC-NA scheme with clipping degrades solution accuracy to first order. The NC schemes converge with second order on regular and/or triangular grids and with first order on perturbed quadrilaterals and mixed-element grids. All schemes may produce large relative errors in gradient reconstruction on grids with perturbed nodes. Defect-correction iterations for schemes employing weighted least-square gradient reconstruction diverge on perturbed stretched grids. Overall, the CC-NN and CC-SA schemes offer the best options of the lowest complexity and secondorder discretization errors. On anisotropic grids over a curved body typical of turbulent flow simulations, the discretization errors converge with second order and are small for the CC-NN, CC-SA, and CC-FA schemes on all grids and for NC schemes on triangular grids; the discretization errors of the CC-NA scheme without clipping do not converge on irregular grids. Accurate gradient reconstruction can be achieved by introducing a local approximate mapping; without approximate mapping, only the NC scheme with weighted LSQ method provides accurate gradients. Defect correction iterations for the CC-NA scheme without clipping diverge; for the NC scheme with weighted LSQ method, the iterations either diverge or converge very slowly. The best option in curved geometries is the CC-SA scheme that offers low complexity, second-order discretization errors, and fast convergence.

Diskin, Boris↗

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and complexity are studied for four nominally second-order accurate schemes: a node-centered scheme and three cell-centered schemes - a node-averaging scheme and two schemes with nearest-neighbor and adaptive compact stencils for least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Tests from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The tests of the second class are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes may degenerate on mixed grids, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to that of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping based on a distance function commonly available in practical schemes or modifying the scheme stencil to reflect the direction of strong coupling. The major conclusion is that accuracies of the node centered and the best cell-centered schemes are comparable at equivalent number of degrees of freedom.

Diskin, Boris↗

High-Order Stabilized Finite Elements on Dynamic Meshes

The development of dynamic mesh capability for turbulent flow simulations using the Streamlined Upwind Petrov-Galerkin (SUPG) discretization is described. The current work extends previous research to include high-order spatial accuracy, including the satisfaction of the discrete geometric conservation law (GCL) on curved elements. Two closely-related schemes are described and the ability of these schemes to satisfy the GCL, while also maintaining temporal accuracy and conservation is assessed. Studies indicate that although one scheme discretizes the time derivative in conservative form, both schemes exhibit temporal conservation errors that decrease according to the expected design order of accuracy. The source of the temporal conservation errors is examined, and it is demonstrated that many finite-volume and finite-element schemes can also be expected to have difficulty strictly satisfying conservation in time. The effects on conservation are examined and, while present in the simulations, are seen to be negligible for the problems considered.

Anderson, W. Kyle↗

Multiple Scale Simulation for Transitional and Turbulent Flow

A new concept, Multiple Scale Simulation (MSS), is presented in this paper. The basic idea is that the flow is decomposed into several component groups according to spatial and temporal length scales. Each group has its own subdomain, governing system, mesh size, and discretization method. The simulation is then performed groupwise. This approach has been successfully applied in combination with the intergrid dissipation technique for simulation of transitional and turbulent flow in 3-D boundary layers, and it is feasible for 3-D airfoils and other more complex configurations. MSS should prove to ameliorate the scale problems associated with conventional direct numerical simulation.

Liu, Chaoqun↗

Model consistency in large eddy simulation of turbulent channel flows

Combinations of filters and subgrid scale stress models for large eddy simulation of the Navier-Stokes equations are examined by a priori tests and numerical simulations. The structure of the subgrid scales is found to depend strongly on the type of filter used, and consistency between model and filter is essential to ensure accurate results. The implementation of consistent combinations of filter and model gives more accurate turbulence statistics than those obtained in previous investigations in which the models were chosen independently from the filter. Results and limitations of the a priori test are discussed. The effect of grid refinement is also examined.

Piomelli, Ugo↗

Turbulent Vortex-Flow Simulation Over a 65 deg Sharp and Blunt Leading-Edge Delta Wing at Subsonic Speeds

Turbulent thin-layer, Reynolds-Averaged Navier-Stokes solutions, based on a multi-block structured grid, are presented for a 65 deg delta wing having either a sharp leading edge (SLE) or blunt leading edge (BLE) geometry. The primary objective of the study is to assess the prediction capability of the method for simulating the leading-edge flow separation and the ensuing vortex flow characteristics. Computational results are obtained for two angles of attack of approximately 13 and 20 deg, at free-stream Mach number of 0.40 and Reynolds number of 6 million based on the wing mean aerodynamic chord. The effects of two turbulence models of Baldwin-Lomax with Degani-Schiff (BL/DS) and the Spalart-Allmaras (SA) on the numerical results are also discussed. The computations also explore the effects of two numerical flux-splitting schemes, i.e., flux difference splitting (fds) and flux vector splitting (fvs), on the solution development and convergence characteristics. The resulting trends in solution sensitivity to grid resolution for the selected leading-edge geometries, angles of attack, turbulence models and flux splitting schemes are also presented. The validity of the numerical results is evaluated against a unique set of experimental wind-tunnel data that was obtained in the National Transonic Facility at the NASA Langley Research Center.

Ghaffari, Farhad↗

Model consistency in the large eddy simulation of turbulent channel flows

Various combinations of filters and subgrid scale stress models for large eddy simulation of the Navier-Stokes equations are studied by a priori tests and numerical simulations. Consistency between model and filter is found to be essential to ensure accurate results. Results and limitations of the a priori test are discussed. The effect of grid refinement is also examined.

Moin, Parviz↗

The use of direct numerical simulation in the study of turbulent, chemically-reacting flows

Full turbulence simulations are used here to study the effects of chemical heat release on the large-scale structures in turbulent mixing layers. In agreement with laboratory results, it is found that the heat release lowers the rate at which the mixing layer grows and reduces the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer produce changes in the flame vortex structure that act to produce more diffuse vortices than in the constant density case, resulting in lower rotation rates of the large-scale structures. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layer are shown to result from vorticity generation by baroclinic torques. Calculations of the energy in the various wavenumbers shows that the heat release has a stabilizing effect on the growth rate of individual modes. This methodology can be applied for Reynolds numbers less than several hundred and for Damkoehler numbers less than about ten.

Riley, J. J.↗

Compressibility effects on the growth and structure of homogeneous turbulent shear flow

Direct numerical simulations of compressible homogeneous turbulent shear flow are used to provide insight into compressibility effects on turbulence. The simulations show a reduction in the growth rate of the turbulence compared to the incompressible case. Examination of the turbulent kinetic energy budget shows that the reduced growth rate is due to an increase in the dissipation rate due to the divergence of the velocity and to the pressure-dilatation correlation which acts to transfer energy between internal energy and kinetic energy. The structure of the turbulence is also examined. Visualizations of the flowfields reveal the presence of eddy shocklets. These shock structures are important contributors to the increased dissipation of compressible turbulence. A mechanism for the generation of the shocks is suggested.

Blaisdell, G. A.↗

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Raman, Venkatramanan↗