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 649 records · Page 36

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗

Control of flow separation in airfoil/wing design applications

Existing aerodynamic design methods have generally concentrated on the optimization of airfoil or wing shapes to produce a minimum drag while satisfying some basic constraints such as lift, pitching moment, or thickness. Since the minimization of drag almost always precludes the existence of separated flow, the evaluation and validation of these design methods for their robustness and accuracy when separated flow is present has not been aggressively pursued. However, two new applications for these design tools may be expected to include separated flow and the issues of aerodynamic design with this feature must be addressed. The first application of the aerodynamic design tools is the design of airfoils or wings to provide an optimal performance over a wide range of flight conditions (multipoint design). While the definition of 'optimal performance' in the multipoint setting is currently being hashed out, it is recognized that given a wide enough range of flight conditions, it will not be possible to ensure a minimum drag constraint at all conditions, and in fact some amount of separated flow (presumably small) may have to be allowed at the more demanding flight conditions. Thus a multipoint design method must be tolerant of the existence of separated flow and may include some controls upon its extent. The second application is in the design of wings with extended high speed buffet boundaries of their flight envelopes. Buffet occurs on a wing when regions of flow separation have grown to the extent that their time varying pressures induce possible destructive effects upon the wing structure or adversely effect either the aircraft controllability or the passenger comfort. A conservative approach to the expansion of the buffet flight boundary is to simply expand the flight envelope of nonseparated flow under the assumption that buffet will also thus be alleviated. However, having the ability to design a wing with separated flow and thus to control the location, extent, and severity of the separated flow regions may allow aircraft manufacturers to gain an advantage in the early design stages of an aircraft, when configuration changes are relatively inexpensive to make. Continuing the work begun last year, an airfoil design package has been modified to provide some control over the existence and extent of flow separation. This package consists of a 2-D Navier-Stokes flow solver which is coupled to the CDISC (constrained direct/iterative surface curvature) design method. The first modification is a prediction method for determining whether separation is likely based solely upon a given pressure distribution. If separation is predicted but is undesirable, the new routines will modify the pressure distribution to alleviate the problem. This new pressure distribution is then used in the design method to generate a new aerodynamic shape. Since separation may be acceptable in some cases, particularly if the separation does not extend to the trailing edge, another added logic estimates the extent of separation based upon a correlation with calculated separated flow cases. If the flow behind a shock induced separation is not predicted to reattach before the trailing edge, the logic weakens the shock strength and otherwise alters the pressure distribution in order to promote reattachment. This later addition is as yet unreliable due to secondary separation effects, but additional work is being pursued to improve the method.

Gally, Thomas A.↗

Uncertainty Quantification of Turbulence Model Closure Coefficients for Transonic Wall-Bounded Flows

The goal of this work was to quantify the uncertainty and sensitivity of commonly used turbulence models in Reynolds-Averaged Navier-Stokes codes due to uncertainty in the values of closure coefficients for transonic, wall-bounded flows and to rank the contribution of each coefficient to uncertainty in various output flow quantities of interest. Specifically, uncertainty quantification of turbulence model closure coefficients was performed for transonic flow over an axisymmetric bump at zero degrees angle of attack and the RAE 2822 transonic airfoil at a lift coefficient of 0.744. Three turbulence models were considered: the Spalart-Allmaras Model, Wilcox (2006) k-w Model, and the Menter Shear-Stress Trans- port Model. The FUN3D code developed by NASA Langley Research Center was used as the flow solver. The uncertainty quantification analysis employed stochastic expansions based on non-intrusive polynomial chaos as an efficient means of uncertainty propagation. Several integrated and point-quantities are considered as uncertain outputs for both CFD problems. All closure coefficients were treated as epistemic uncertain variables represented with intervals. Sobol indices were used to rank the relative contributions of each closure coefficient to the total uncertainty in the output quantities of interest. This study identified a number of closure coefficients for each turbulence model for which more information will reduce the amount of uncertainty in the output significantly for transonic, wall-bounded flows.

Schaefer, John↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗

RANS and DES Computations for a Wing with Ice Accretion

A computational investigation was performed to assess the effectiveness of Detached Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The wing section considered here was a business jet airfoil (GLC305) with a 22.5-minute glaze ice accretion (944-ice shape). The section was extruded to form a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. More validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy.

Mogili, P.↗

Numerical Study Comparing RANS and LES Approaches on a Circulation Control Airfoil

A numerical study over a nominally two-dimensional circulation control airfoil is performed using a large-eddy simulation code and two Reynolds-averaged Navier-Stokes codes. Different Coanda jet blowing conditions are investigated. In addition to investigating the influence of grid density, a comparison is made between incompressible and compressible flow solvers. The incompressible equations are found to yield negligible differences from the compressible equations up to at least a jet exit Mach number of 0.64. The effects of different turbulence models are also studied. Models that do not account for streamline curvature effects tend to predict jet separation from the Coanda surface too late, and can produce non-physical solutions at high blowing rates. Three different turbulence models that account for streamline curvature are compared with each other and with large eddy simulation solutions. All three models are found to predict the Coanda jet separation location reasonably well, but one of the models predicts specific flow field details near the Coanda surface prior to separation much better than the other two. All Reynolds-averaged Navier-Stokes computations produce higher circulation than large eddy simulation computations, with different stagnation point location and greater flow acceleration around the nose onto the upper surface. The precise reasons for the higher circulation are not clear, although it is not solely a function of predicting the jet separation location correctly.

Rumsey, Christopher L.↗

Assessment of LAURA for Laminar Supersonic Shallow Cavities

The ability of the Laura flow solver to predict local heating augmentation factors for shallow cavities is assessed. This assessment is part of a larger e ort within the Space Shuttle return-to-flight program to develop technologies to support on-orbit tile repair decisions. The comparison is made against global phosphor thermography images taken in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel. The cavities are rectangular in shape, with lengths L/H of 14 20 and depths H/ of 1.1 5.2. The fully laminar results, for Re = 300, show good agreement between the data sets. For Re = 503, the wind tunnel data indicates boundary layer transition with turbulent flow both within and downstream of the cavity. The turbulent flow structures are significantly di erent from the laminar predictions, with order of magnitude increases in the heating augmentations. Because of the di erent flow structures, no simple bump factor can be used to correct the laminar calculations to account for the turbulent heating levels. A fine gradation in wind tunnel cases will be required to clearly delineate the laminar-to-turbulent transition point, and hence the limits of applicability of the laminar numerical approach.

Wood, William A.↗

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul↗

USM3D-ME Revised SA-neg Implementation and Assessment for Seventh Drag Prediction Workshop Solutions

Simulations performed by an unstructured-grid, cell-centered, finite-volume, mixed-element flow solver USM3D-MEwere presented at the Seventh Drag Prediction Workshop (DPW-7). Reynolds-averaged Navier-Stokes (RANS) solutions for the NASA high-speed common research model wing-body configuration were computed using a new implementation of the negative variant of the Spalart-Allmaras turbulence model (SA-neg) with hierarchical adaptive nonlinear iteration method (HANIM).ThisUSM3D-ME enhancement improved efficiency and robustness of RANS solutions resulting in fast iterative convergence to machine zero residuals on mixed-element unstructured grids provided by the DPW-7 organizing committee. This paper presents theoretical foundations for the new implementation of the SA-neg model and analyzesUSM3D-ME solutions for flow conditions corresponding to the required DPW-7 cases.

Boris Diskin↗

High-Lift OVERFLOW Analysis of the DLR-F11 Wind Tunnel Model

In response to the 2nd AIAA CFD High Lift Prediction Workshop, the DLR-F11 wind tunnel model is analyzed using the Reynolds-averaged Navier-Stokes flow solver OVERFLOW. A series of overset grids for a bracket-off landing configuration is constructed and analyzed as part of a general grid refinement study. This high Reynolds number (15.1 million) analysis is done at multiple angles-of-attack to evaluate grid resolution effects at operational lift levels as well as near stall. A quadratic constitutive relation recently added to OVERFLOW for improved solution accuracy is utilized for side-of-body separation issues at low angles-of-attack and outboard wing separation at stall angles. The outboard wing separation occurs when the slat brackets are added to the landing configuration and is a source of discrepancy between the predictions and experimental data. A detailed flow field analysis is performed at low Reynolds number (1.35 million) after pressure tube bundles are added to the bracket-on medium grid system with the intent of better understanding bracket/bundle wake interaction with the wing's boundary layer. Localized grid refinement behind each slat bracket and pressure tube bundle coupled with a time accurate analysis are exercised in an attempt to improve stall prediction capability. The results are inconclusive and suggest the simulation is missing a key element such as boundary layer transition. The computed lift curve is under-predicted through the linear range and over-predicted near stall, and the solution from the most complete configuration analyzed shows outboard wing separation occurring behind slat bracket 6 where the experiment shows it behind bracket 5. These results are consistent with most other participants of this workshop.

High-lift↗

A fully 3-D numerical model to predict flood wave propagation and assess efficiency of flood protection measures

Discussed are the development and validation of a 3-D, non-hydrostatic, Reynolds-Averaged Navier–Stokes (RANS) model using the Volume of Fluid (VOF) approach to simulate flooding events in natural environments. A series of validation simulations in simple and complex geometries are performed to test the different modules of the flow solver. The paper also discusses the differences between simulation results obtained using the 3-D model and those obtained using a standard, 2-D hydrostatic, depth-averaged model. The predictions by the two models show fairly large differences in regions where 3-D effects are expected to be significant (e.g., in regions containing highly-curved river reaches and near hydraulic structures). Then, the validated 3-D model is employed to investigate the efficiency of using flood-protection walls to reduce adverse effects of flooding in the Iowa River near Iowa City, Iowa, USA. Finally, the propagation of a flood wave in a reach of the Iowa River containing a bridge is simulated and the effects of the flow becoming pressurized beneath the bridge deck are discussed. For such applications, in which the regime can change locally between open channel and pressurized flow, the use of fully 3-D, non-hydrostatic models may be the only viable option to accurately predict flow hydrodynamics and its effects on bed shear stresses and sediment entrainment potential, especially around locations where the flow becomes pressurized.

Water Resources↗

CMLM (Co-Optimized Machine-Learned Manifolds) [SWR-23-41]

Co-optimized Machine-Learned Manifolds (CMLM) is a data-driven approach for developing reduced-order manifold models for high-dimensional chemically reacting systems. It involves a specially designed neural network, the training of which simultaneously optimizes linear combinations of species that define the manifold, nonlinear mapping to outputs of interest such as reaction rates, and (optionally) subfilter closure for large eddy simulation. This software package provides an implementation of the CMLM approach in Python using the PyTorch machine learning library. A few example cases are included, showing how the tool can be applied to different types of data from 0D and 1D reacting simulations performed using Cantera. The neural networks can be saved in a format that is readable by the Pele suite of combustion solvers for use in reacting computational fluid dynamics simulations. This software repository contains several python scripts to perform various tasks associated with the Co-optimized Machine Learned Manifolds (CMLM) model, which is described in Perry, Henry de Frahan, and Yellapantula, CNF, 2022 (https://doi.org/10.1016/j.combustflame.2022.112286). This includes not only the code that defines the CMLM model, but also scripts to generate suitable training data, scripts to pre-process the data, scripts to train the CMLM model, and scripts to plot the output, as well as various other helper files. The scripts depend on several commonly used python libraries for data analysis and chemical reaction computations. The trained models that result from this tool are designed to work with the an interface being implemented in the Pele suite of reacting flow solvers (https://github.com/AMReX-Combustion).

Perry, Bruce↗

Numerical Analysis of Wave Characteristics in a Methane-Oxygen Rotating Detonation Engine

Results from the simulation of the U.S. Air Force Research Laboratory methane–oxygen rotating detonation rocket engine from four independent research groups with different flow solvers underpredicted primary detonation wave speeds by a significant margin as compared to experimental values. In a simulation performed by the authors, the average calculated speed of the detonation waves was roughly [Formula: see text] as compared to the experimentally measured value of [Formula: see text]. This paper presents a detailed analysis of the wave characteristics in this simulation to provide a more quantitative understanding of the underlying factors leading to this discrepancy. The results show that weaker counter-rotating shock waves have a significant impact on the behavior of the primary detonation waves. The wave speed of the primary detonation wave is reduced by [Formula: see text] due to collisions with the counter-rotating waves. The presence of these counter-rotating waves also has a strong influence on the flow conditions upstream of the primary detonation waves, as well as the engine heat release rate. In addition, the flow properties upstream of the detonation wave vary significantly in the radial direction. Finally, the subfilter turbulent viscosity is shown to vary radially and with proximity to detonation waves.

Engineering↗

String flash-boiling in gasoline direct injection simulations with transient needle motion

A computational study was performed to investigate the influence of transient needle motion on gasoline direct injection (GDI) internal nozzle flow and near-field sprays. Simulations were conducted with a compressible Eulerian flow solver modeling liquid, vapor, and non-condensable gas phases with a diffuse interface. Variable rate generation and condensation of fuel vapor were captured using the homogeneous relaxation model (HRM). The non-flashing (spray G) and flashing (spray G2) conditions specified by the Engine Combustion Network were modeled using the nominal spray G nozzle geometry and transient needle lift and wobble were based upon ensemble averaged x-ray imaging preformed at Argonne National Lab. The minimum needle lift simulated was 5 μm and dynamic mesh motion was achieved with Laplacian smoothing. The results were qualitatively validated against experimental imaging and the experimental rate of injection profile was captured accurately using pressure boundary conditions and needle motion to actu- ate the injection. Needle wobble was found to have no measurable effect on the flow. Low needle lift is shown to result in vapor generation as fuel rushes past the needle. In conclusion, the internal injector flow is shown to contain many transient and interacting vortices which cause perturbations in the spray angle, fluctuations in the mass flux, and frequently result in string flash-boiling.

33 ADVANCED PROPULSION SYSTEMS↗

Failure of Anisotropic Unstructured Mesh Adaption Based on Multidimensional Residual Minimization

An automated anisotropic unstructured mesh adaptation strategy is proposed, implemented, and assessed for the discretization of viscous flows. The adaption criteria is based upon the minimization of the residual fluctuations of a multidimensional upwind viscous flow solver. For scalar advection, this adaption strategy has been shown to use fewer grid points than gradient based adaption, naturally aligning mesh edges with discontinuities and characteristic lines. The adaption utilizes a compact stencil and is local in scope, with four fundamental operations: point insertion, point deletion, edge swapping, and nodal displacement. Evaluation of the solution-adaptive strategy is performed for a two-dimensional blunt body laminar wind tunnel case at Mach 10. The results demonstrate that the strategy suffers from a lack of robustness, particularly with regard to alignment of the bow shock in the vicinity of the stagnation streamline. In general, constraining the adaption to such a degree as to maintain robustness results in negligible improvement to the solution. Because the present method fails to consistently or significantly improve the flow solution, it is rejected in favor of simple uniform mesh refinement.

Wood, William A.↗

A GPU-based 2D shallow water quality model

In this paper, a 2D shallow water flow solver integrated with a water quality model is presented. The interaction between the main water quality constituents included is based on the Water Quality Analysis Simulation Program. Efficiency is achieved by computing with a combination of a Central Processing Unit (CPU) and a Graphics Processing Unit (GPU) device. This technique is intended to provide robust and accurate simulations with high computation speedups with respect to a single-core CPU in real events. The proposed numerical model is evaluated in cases that include the transport and reaction of water quality components over irregular bed topography and dry–wet fronts, verifying that the numerical solution in these situations conserves the required properties (C-property and positivity). The model can operate in any steady or unsteady form allowing an efficient assessment of the environmental impact of water flows. The field data from an unsteady river reach test case are used to show that the model is capable of predicting the measured temporal distribution of dissolved oxygen and water temperature, proving the robustness and computational efficiency of the model, even in the presence of noisy signals such as wind speed.

54 ENVIRONMENTAL SCIENCES↗

HLPW-5: Overview and Workshop Summary

The Fifth AIAA CFD High-Lift Prediction Workshop was held with the goal of assessing the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in high-lift configurations. A key aspect of this endeavor was the use of Technology Focus Groups (TFG), which included both mesh generation and flow-solver experts working together to accelerate advancements for their particular CFD methodology, by addressing key questions of importance prior to the workshop. The high-lift version of the NASA Common Research Model (CRM-HL) configuration was the focus of this workshop, and was used for three unique test cases. Wind-tunnel data were available for comparison for one of the test cases. Altogether, 365 datasets of CFD results were submitted by 47 teams, with 41 teams contributing to the multiple configurations of Case 1, 40 to Case 2, and 18 to Case 3. This paper provides a high-level summary of the results and conclusions from the workshop. As concluded from past workshops, fixed-grid Reynolds-averaged Navier-Stokes methods continued to be inaccurate and inconsistent for high-lift flows near maximum lift. However, application of mesh-adaptation technology helped to achieve improved consistency. Scale-resolving methods appeared most promising for predicting high-lift flow physics, particularly at maximum lift. Best practices for these methods were refined over the course of the workshop and new challenges were identified.

Adam M Clark↗

Design and Testing of a Blended Wing Body With Boundary Layer Ingestion Nacelles at High Reynolds Numbers

A knowledge-based aerodynamic design method coupled with an unstructured grid Navier-Stokes flow solver was used to improve the propulsion/airframe integration for a Blended Wing Body with boundary-layer ingestion nacelles. A new zonal design capability was used that significantly reduced the time required to achieve a successful design for each nacelle and the elevon between them. A wind tunnel model was built with interchangeable parts reflecting the baseline and redesigned configurations and tested in the National Transonic Facility (NTF). Most of the testing was done at the cruise design conditions (Mach number = 0.85, Reynolds number = 75 million). In general, the predicted improvements in forces and moments as well as the changes in wing pressures between the baseline and redesign were confirmed by the wind tunnel results. The effectiveness of elevons between the nacelles was also predicted surprisingly well considering the crudeness in the modeling of the control surfaces in the flow code. A novel flow visualization technique involving pressure sensitive paint in the cryogenic nitrogen environment used in high-Reynolds number testing in the NTF was also investigated.

Campbell, Richard L.↗