Search NASA⌕ Search

SEARCH · Search NASA

Results for “wall modeling”

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

DNS/LES Study of Representative Wall-Bounded Turbulent Flows using SIERRA/Fuego

This report summarizes a series of SIERRA/Fuego validation efforts of turbulent flow models on canonical wall-bounded configurations. In particular, direct numerical simulations (DNS) and large eddy simulations (LES) turbulence models are tested on a periodic channel, a periodic pipe, and an open jet for which results are compared to the velocity profiles obtained theoretically or experimentally. Velocity inlet conditions for channel and pipe flows are developed for application to practical simulations. To show this capability, LES is performed over complex terrain in the form of two natural hills and the results are compared with other flow solvers. The practical purpose of the report is to document the creation of inflow boundary conditions of fully developed turbulent flows for other LES calculations where the role of inflow turbulence is critical.

42 ENGINEERING↗

Comparison of a two-dimensional adaptive-wall technique with analytical wall interference correction techniques

A two dimensional airfoil model was tested in the adaptive wall test section of the NASA Langley 0.3 meter Transonic Cryogenic Tunnel (TCT) and in the ventilated test section of the National Aeronautical Establishment Two Dimensional High Reynold Number Facility (HRNF). The primary goal of the tests was to compare different techniques (adaptive test section walls and classical, analytical corrections) to account for wall interference. Tests were conducted over a Mach number range from 0.3 to 0.8 at chord Reynolds numbers of 10 x 10(exp 6), 15 x 10(exp 6), and 20 x 10(exp 6). The angle of attack was varied from about 12 degrees up to stall. Movement of the top and bottom test section walls was used to account for the wall interference in the HRNF tests. The test results are in good agreement.

Mineck, Raymond E.↗

Euler and Navier-Stokes solutions for hypersonic flows

An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.

Thareja, Rajiv R.↗

Radiative heat transfer in rocket thrust chambers and nozzles

Numerical models based on the discrete exchange factor (DEF) and the zonal methods for radiative analysis of rocket engines containing a radiatively participating medium have been developed. These models implement a new technique for calculating the direct exchange factors to account for possible blockage by the nozzle throat. Given the gas and surface temperature distributions, engine geometry, and radiative properties, the models compute the wall radiative heat fluxes at different axial positions. The results of sample calculations for a typical rocket engine (engine 700 at NASA), which uses RP-1 (a kerosene-type propellant), are presented for a wide range of surface and gas properties. It is found that the heat transfer by radiation can reach up to 50 percent of that due to convection. The maximum radiative heat flux is at the inner side of the engine, where the gas temperature is the highest. While the results of both models are in excellent agreement, the computation time of the DEF method is found to be much smaller.

Hammad, K. J.↗

Numerical computation of aerodynamics and heat transfer in a turbine cascade and a turn-around duct using advanced turbulence models

The objective of this research is to develop turbulence models to predict the flow and heat transfer fields dominated by the curvature effect such as those encountered in turbine cascades and turn-around ducts. A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two layer k-epsilon/ARSM (Algebraic Reynolds Stress Model), Chien's Low Reynolds Number (LRN) k-epsilon model and Coakley's LRN q-omega model. The near wall pressure strain correlation term was included in the ARSM. The formulation is applied to Favre-averaged N-S equations and no thin-layer approximations are made in either the mean flow or turbulence transport equations. Anisotropic scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve convergence. Detailed comparisons were made between computations and data measured in a turbine cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles were predicted well by all the models. The blade heat transfer is predicted well by k-epsilon/ARSM model, as well as the k-epsilon model. It's found that the onset of boundary layer transition on both surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly influenced by the streamline curvature. Detailed computation of the flow in the turn around duct has been carried out and validated against the data by Monson as well as Sandborn. The computed results at various streamwise locations both on the concave and convex sides are compared with flow and turbulence data including the separation zone on the inner well. The k-epsilon/ARSM model yielded relatively better results than the two-equation turbulence models. A detailed assessment of the turbulence models has been made with regard to their applicability to curved flows.

Lakshminarayana, B.↗

Mid-Fidelity Computational Fluid Dynamics Analysis of the Elytron 4S UAV Concept

The Elytron 4S Unmanned Aerial Vehicle (UAV) Concept was developed to combine the advantages of fixed- and rotary-wing technology. The 4S Concept is a box-wing configuration with rotors mounted on a centrally located tiltwing. The UAV is intended to be capable of both conventional takeoff and landing (CTOL) and vertical takeoff and landing (VTOL), and is envisioned to excel in UAV performance because of the combined efficiency of fixed-wing aircraft and the hover and VTOL capabilities of regular drones or quadcopters. A mid-fidelity Unsteady Reynolds-Averaged Navier-Stokes (URANS) approach using Rotorcraft CFD (RotCFD) is performed to analyze and characterize the performance of the aircraft. The flow field is coupled with a rotor model based on blade-element momentum theory to model the 4S UAV rotors. Turbulence is modeled using a realizable k- turbulence model with special wall function. The code is used to generate aerodynamic forces and moments on the body at cruise conditions, and during VTOL. The results and their uncertainties are characterized, and an angle- of-attack and sideslip sweep are computed, both with and without rotors on. Simulations are compared with the wind tunnel tests in the 7- by 10-Foot U.S. Army Wind Tunnel at NASA Ames Research Center, performed in 2017. Results show promising comparison with experimental data, despite a late change in rotor size and rudder size of the physical model that cause the expected deviations from the simulation. A slight change in the net thrust value, when rotors are modeled, is observed because of the rotor diameter increase on the physical model. A noticeable difference in the directional stability was observed because of the increased rudder surface and added strakes. These changes were implemented to improve on the design as simulated, which is observed in the results. The simulation results paved the way to the first successful flight of the UAV Concept.

Koning, Witold J. K.↗

Direct simulation of compressible wall-bounded turbulence

When analyzing many turbulent flows, the effects of compressibility can be neglected. Even some relatively high-speed flows, such as boundary layers generated by a supersonic aircraft, produce turbulent statistics that are similar to those found for the incompressible case. There are other situations, however, in which the non-zero divergence of the turbulence leads to behavior that is fundamentally different from that found at constant density. Examples include flows created by internal combustion engines, hypersonic flight, and supersonic combustion. It is with instances such as these that this project is concerned. In particular, we are interested in the effects of compressibility on turbulence near a smooth solid constant temperature surface; our primary objective is an increased physical understanding that can be used to improve turbulence models of wall-bounded compressible flows. With this in mind, we have begun a direct numerical simulation (DNS) study of turbulence in a plane channel. Because all of the relevant spatial and temporal scales are to be resolved, the simulations require no subgrid scale parameterization. The DNS code developed by Buell to study compressible plane Couette flow was modified to solve the compressible Navier-Stokes equations in the plane channel. The channel was chosen over the Couette flow for two reasons: (1) to avoid using the very large streamwise domains needed to adequately capture the large Couette vortical structures; and (2) to make use of previous experience by considering the compressible version of a well established case and isolate finite Mach number effects by comparing to the incompressible channel.

Coleman, Gary N.↗

On the application of algebraic turbulence models to high Mach number flows

In this paper the implementation of embedded-function algorithms for conventional algebraic turbulence models is described; using this methodology it is shown that the entire wall layer may be removed from the calculation algorithm (as well as near-wall mixing-length models) with no degradation in the accuracy of predicted results for skin friction. Various difficulties that are associated with the application of conventional algebraic models in high-speed flows are discussed. Alternative algebraic models are proposed which are cast in terms of the Howarth-Dorodnitsyn variable and which lead to a simple 'law of the wall' for both velocity and total enthalpy.

Degani, A. T.↗

Numerical analysis of hypersonic turbulent film cooling flows

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

Chen, Y. S.↗

The self streamlining wind tunnel

A two dimensional test section in a low speed wind tunnel capable of producing flow conditions free from wall interference is presented. Flexible top and bottom walls, and rigid sidewalls from which models were mounted spanning the tunnel are shown. All walls were unperforated, and the flexible walls were positioned by screw jacks. To eliminate wall interference, the wind tunnel itself supplied the information required in the streamlining process, when run with the model present. Measurements taken at the flexible walls were used by the tunnels computer check wall contours. Suitable adjustments based on streamlining criteria were then suggested by the computer. The streamlining criterion adopted when generating infinite flowfield conditions was a matching of static pressures in the test section at a wall with pressures computed for an imaginary inviscid flowfield passing over the outside of the same wall. Aerodynamic data taken on a cylindrical model operating under high blockage conditions are presented to illustrate the operation of the tunnel in its various modes.

Goodyer, M. J.↗

Adaptive wall testing sections (AWTS)

The lecture starts with conventional techniques of minimizing wall interference and explains the principle of wall streamlining. The history of AWTS development is highlighted, along with the benefits of wall streamlining, including minimized boundary interference, increased model size, reduced tunnel drive power, noise, and volume, as well as multiple flow field simulations to be performed using one test section. AWTS-associated problems coming from the need to adjust the test-section boundaries for each test condition are assessed, along with the requirements of a boundary-adjustment strategy. Examples of two- and three-dimensional test sections are presented, and attention is focused on residual interference and the effects of compressibility and model lift on flexible-wall contours.

Wolf, Stephen W. D.↗

Navier-Stokes computation of compressible turbulent flows with a second order closure, part 1

A second order closure turbulence model for compressible flows is developed and implemented in a 2D Reynolds-averaged Navier-Stokes solver. From the beginning where a kappa-epsilon turbulence model was implemented in the bidiagonal implicit method of MACCORMACK (referred to as the MAC3 code) to the final stage of implementing a full second order closure in the efficient line Gauss-Seidel algorithm, numerous work was done, individually and collectively. Besides the collaboration itself, the final product of this work is a second order closure derived from the Launder, Reece, and Rodi model to account for near wall effects, which has been called FRAME model, which stands for FRench-AMerican-Effort. During the reporting period, two different problems were worked out. The first was to provide Ames researchers with a reliable compressible boundary layer code including a wide collection of turbulence models for quick testing of new terms, both in two equations and in second order closure (LRR and FRAME). The second topic was to complete the implementation of the FRAME model in the MAC5 code. The work related to these two different contributions is reported. dilatation in presence of stron shocks. This work, which has been conducted during a work at the Center for Turbulence Research with Zeman aimed also to cros-check earlier assumptions by Rubesin and Vandromme.

Haminh, Hieu↗

Wall interference correction improvements for the ONERA main wind tunnels

This paper describes improved methods of calculating wall interference corrections for the ONERA large windtunnels. The mathematical description of the model and its sting support have become more sophisticated. An increasing number of singularities is used until an agreement between theoretical and experimental signatures of the model and sting on the walls of the closed test section is obtained. The singularity decentering effects are calculated when the model reaches large angles of attack. The porosity factor cartography on the perforated walls deduced from the measured signatures now replaces the reference tests previously carried out in larger tunnels. The porosity factors obtained from the blockage terms (signatures at zero lift) and from the lift terms are in good agreement. In each case (model + sting + test section), wall corrections are now determined, before the tests, as a function of the fundamental parameters M, CS, CZ. During the windtunnel tests, the corrections are quickly computed from these functions.

Vaucheret, X.↗

A transport equation for eddy viscosity

A transport equation for eddy viscosity is proposed for wall bounded turbulent flows. The proposed model reduces to a quasi-homogeneous form far from surfaces. Near to a surface, the nonhomogeneous effect of the wall is modeled by an elliptic relaxation model. All the model terms are expressed in local variables and are coordinate independent; the model is intended to be used in complex flows. Turbulent channel flow and turbulent boundary layer flows with/without pressure gradient are calculated using the present model. Comparisons between model calculations and direct numerical simulation or experimental data show good agreement.

Durbin, P. A.↗

Comparative Study of Advanced Turbulence Models for Turbomachinery

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been independently tested by Professor C.P. Chen and his group at the University of Alabama at Huntsville (UAH) by interfacing them with own flow solver (MAST).

Hadid, Ali H.↗

Direct assessment of two-dimensional wind-tunnel interference from measurements on two interfaces

A direct assessment of two-dimensional wind-tunnel wall interference using upwash component measured on two interfaces has been formulated by the Prandtl-Glauert equation of the flow field and solved by the Fourier transform technique. The analytic formulae obtained for the interference of upwash and pressure on the model are presented. The formulae have been applied successfully to the analytic models for lifting and blockage interferences induced by the general linear slotted and perforated tunnel wall boundary conditions. The formulae have been derived in terms of Fourier coefficients of the measured upwash and illustrated its merits applying to a wavy-wall model case.

Lo, Ching F.↗

A deep learning upscaling framework: Reactive transport and mineral precipitation in fracture-matrix systems

Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes for which effective medium approximations may become inaccurate, e.g., mineral precipitation and clogging during reactive transport. We have developed a deep learning upscaling framework, in which pore-scale modeling is directly employed in macroscale systems, without relying on effective medium approximations. The upscaling framework is first developed for general multiscale systems and then applied to modeling reactive transport with mineral precipitation in the altered layer in fracture-matrix structures. Solute transport from the fractures to the matrix is modeled as a wall boundary condition for the fractures, which, in turn, is predicted by recurrent neural networks using the concentration histories at the fracture-matrix boundary. Specifically, we consider a meter-scale fracture network embedded in sandstones, where the smallest feature is at the micron scale. Here the proposed framework allows us to span five orders of magnitude in length scales by capturing mineral precipitation in the altered layer of the rock matrix at the pore scale across the entire meter-scale fracture network.

42 ENGINEERING↗

Coarse Grid Modeling of Turbine Film Cooling Flows Using Volumetric Source Terms

The recent trend in numerical modeling of turbine film cooling flows has been toward higher fidelity grids and more complex geometries. This trend has been enabled by the rapid increase in computing power available to researchers. However, the turbine design community requires fast turnaround time in its design computations, rendering these comprehensive simulations ineffective in the design cycle. The present study describes a methodology for implementing a volumetric source term distribution in a coarse grid calculation that can model the small-scale and three-dimensional effects present in turbine film cooling flows. This model could be implemented in turbine design codes or in multistage turbomachinery codes such as APNASA, where the computational grid size may be larger than the film hole size. Detailed computations of a single row of 35 deg round holes on a flat plate have been obtained for blowing ratios of 0.5, 0.8, and 1.0, and density ratios of 1.0 and 2.0 using a multiblock grid system to resolve the flows on both sides of the plate as well as inside the hole itself. These detailed flow fields were spatially averaged to generate a field of volumetric source terms for each conservative flow variable. Solutions were also obtained using three coarse grids having streamwise and spanwise grid spacings of 3d, 1d, and d/3. These coarse grid solutions used the integrated hole exit mass, momentum, energy, and turbulence quantities from the detailed solutions as volumetric source terms. It is shown that a uniform source term addition over a distance from the wall on the order of the hole diameter is able to predict adiabatic film effectiveness better than a near-wall source term model, while strictly enforcing correct values of integrated boundary layer quantities.

Heidmann, James D.↗