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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Effect of driven-wall motion on a turbulent boundary layer

An acoustically driven wall model with amplitudes greater than 0.01 cm was used to study the effect of unsteady wall boundary conditions on the structure of a zero-pressure-gradient turbulent boundary layer. Hot-wire measurements showed large variations in fluctuation quantities (including Reynolds stress) along each model wave. These longitudinal variations decreased rapidly with increasing distance from the driven wall, but the average level of these fluctuations decreased only slightly even outside the boundary layer. When the acoustical fluctuations were subtracted, the remaining fluctuations were essentially the same as for the undriven model. It is concluded that there is no significant interaction between the acoustic and turbulent fluctuations.

Weinstein, L. M.↗

NASA NTRS1995

Analysis of a High-Lift Multi-Element Airfoil using a Navier-Stokes Code

A thin-layer Navier-Stokes code, CFL3D, was utilized to compute the flow over a high-lift multi-element airfoil. This study was conducted to improve the prediction of high-lift flowfields using various turbulence models and improved glidding techniques. An overset Chimera grid system is used to model the three element airfoil geometry. The effects of wind tunnel wall modeling, changes to the grid density and distribution, and embedded grids are discussed. Computed pressure and lift coefficients using Spalart-Allmaras, Baldwin-Barth, and Menter's kappa-omega - Shear Stress Transport (SST) turbulence models are compared with experimental data. The ability of CFL3D to predict the effects on lift coefficient due to changes in Reynolds number changes is also discussed.

Whitlock, Mark E.↗

NASA NTRS2015

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Simulation↗

DOE OSTI2023

Modeling of the spray-induced wall stress acting on the ignition assistance device

This research introduces a novel wall-stress model called the Spray-Induced Wall Stress (SIWS) model, which considers the effects of spray-wall impingement and the resulting formation of wall stress within the Lagrangian spray modeling framework. The primary objective of this paper is to provide a mathematical description of the fundamental physics underlying the model. Subsequently, the proposed model is validated using existing experimental data. Here, the remainder of the study focuses on the practical application of the model to an ignition assistance device. Specifically, this device is installed in a compression ignition engine and designed to enhance ignition in aviation-fueled high-altitude aircraft propulsion systems. The research sheds light on the mechanical impulse caused by the high-speed impact of the spray jet, leading to the accumulation of mechanical stress on the rigid body of the ignition assistance device. Previous studies on fluid-structure interaction have only considered the interaction between the gas phase and the solid wall. However, the SIWS model incorporates the additional impact of the impinging liquid spray jet. Consequently, the simulated stress distribution on the ignition assistance device can be estimated by considering both the gas-phase-induced term and the spray-induced term simultaneously.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1971

Numerical calculation of turbulent boundary layers including suction or injection with binary diffusion.

The modelling of turbulent transport processes with an eddy viscosity representation of the inner, wall region for cases with suction or injection is considered in detail. It is shown how the van Driest solid wall model can be logically extended for such cases. Further, a new inductive model is developed. Finite difference calculations using the new inductive model are shown to be in good agreement with experimental data over a wide range of conditions. Lastly, a new method for calculating the wall shear within a finite difference solution is developed. This procedure greatly relaxes the requirement for fine grid spacing in the vertical direction.

Schetz, J. A.↗

NASA NTRS1991

Application of a new K-tau model to near wall turbulent flows

A recently developed K-tau model for near wall turbulent flows is applied to two severe test cases. The turbulent flows considered include the incompressible flat plate boundary layer with the adverse pressure gradients and incompressible flow past a backward facing step. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.

Thangam, S.↗

NASA NTRS1991

Application of a new K-tau model to near wall turbulent flows

A recently developed K-tau model for near wall turbulent flows is applied to a variety of test cases. The turbulent flows considered include the incompressible flat plate boundary layer with adverse pressure gradients, incompressible flow past a backward facing step, and the supersonic flat plate boundary layer at zero pressure gradient. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.

Abid, R.↗

NASA NTRS2001

Turbulence Modeling Workshop

Advances in turbulence modeling are needed in order to calculate high Reynolds number flows near the onset of separation and beyond. To this end, the participants in this workshop made the following recommendations. (1) A national/international database and standards for turbulence modeling assessment should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be taken of other efforts already under-way, such as that of the European Research Community on Flow, Turbulence, and Combustion (ERCOFTAC) consortium. Carefully selected "unit" experiments will be needed, as well as advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least understood and is a key component of two-equation and higher models. Second priority should be given to the development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes (RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in turbulence modeling development, validation, and implementation should include these approaches as well.

Rubinstein, R.↗

DOE OSTI2024

Quantum Barkhausen noise induced by domain wall cotunneling

Most macroscopic magnetic phenomena (including magnetic hysteresis) are typically understood classically. Here, we examine the dynamics of a uniaxial rare-earth ferromagnet deep within the quantum regime, so that domain wall motion, and the associated hysteresis, is initiated by quantum nucleation, which then grows into large-scale domain wall motion, which is observable as an unusual form of Barkhausen noise. We observe noncritical behavior in the resulting avalanche dynamics that only can be explained by going beyond traditional renormalization group methods or classical domain wall models. We find that this “quantum Barkhausen noise” exhibits two distinct mechanisms for domain wall movement, each of which is quantum-mechanical, but with very different dependences on an external magnetic field applied transverse to the spin (Ising) axis. These observations can be understood in terms of the correlated motion of pairs of domain walls, nucleated by cotunneling of plaquettes (sections of domain wall), with plaquette pairs correlated by dipolar interactions; this correlation is suppressed by the transverse field. Similar macroscopic correlations may be expected to appear in the hysteresis of other systems with long-range interactions.

Science & Technology - Other Topics↗

NASA NTRS1990

A critical evaluation of two-equation models for near wall turbulence

A variety of two-equation turbulence models,including several versions of the K-epsilon model as well as the K-omega model, are analyzed critically for near wall turbulent flows from a theoretical and computational standpoint. It is shown that the K-epsilon model has two major problems associated with it: the lack of natural boundary conditions for the dissipation rate and the appearance of higher-order correlations in the balance of terms for the dissipation rate at the wall. In so far as the former problem is concerned, either physically inconsistent boundary conditions have been used or the boundary conditions for the dissipation rate have been tied to higher-order derivatives of the turbulent kinetic energy which leads to numerical stiffness. The K-omega model can alleviate these problems since the asymptotic behavior of omega is known in more detail and since its near wall balance involves only exact viscous terms. However, the modeled form of the omega equation that is used in the literature is incomplete-an exact viscous term is missing which causes the model to behave in an asymptotically inconsistent manner. By including this viscous term and by introducing new wall damping functions with improved asymptotic behavior, a new K-tau model (where tau is identical with 1/omega is turbulent time scale) is developed. It is demonstrated that this new model is computationally robust and yields improved predictions for turbulent boundary layers.

Speziale, Charles G.↗

NASA NTRS1993

Flowfield dynamics in blunt fin-induced shock wave/turbulent boundary layer interactions

Fluctuating wall pressure measurements were made on centerline upstream of a blunt fin in a Much 5 turbulent boundary layer. By examining the ensemble average wall pressure distributions for different fixed shock foot positions, it was shown that local fluctuating wall pressure measurements are due to a distinct pressure distribution, P(i), which undergoes a stretching and flattening effect as its upstream boundary translates aperiodically between the upstream influence and separation lines. The locations of the maxima and minima in the centerline wall standard deviation distribution can be accurately predicted using this distribution, providing quantitative confirmation of the model. This model also explains the observed cross-correlations and ensemble average measurements within the interaction. Using the P(i) model, wall pressure signals from under the separated flow region were able to reproduce the position-time history of the separation shock foot. The negative time delay peak in the cross-correlation between the predicted and actual shock foot histories shows that the separated region fluctuations precede shock foot motion.

Brusniak, L.↗

NASA NTRS1994

Flowfield dynamics in blunt fin-induced shock wave/turbulent boundary layer interactions

Fluctuating wall pressure measurements have been made on centerline upstream of a blunt fin in a Mach 5 turbulent boundary layer. By examining the ensemble averaged wall pressure distributions for different separation shock foot positions, it has been shown that local fluctuating wall pressure measurements are due to a distinct pressure distribution, Rho(sub i), which undergoes a stretching and flattening effect as its upstream boundary translates aperiodically between the upstream influence and separation lines. The locations of the maxima and minima in the wall pressure standard deviation can be accurately predicted using this distribution, providing quantitative confirmation of the model. This model also explains the observed cross-correlations and ensemble average measurements within the interaction. Using the Rho(sub i) model, wall pressure signals from under the separated flow region were used to reproduce the position-time history of the separation shock foot. Further, the negative time delay peak in the cross-correlation between the predicted and actual shock foot histories suggests that the separated region fluctuations precede shock foot motion. The unsteady behavior of the primary horseshoe vortex and its relation to the unsteady separation shock are described.

Dolling, David S.↗

NASA NTRS1990

A critical evaluation of two-equation models for near wall turbulence

A basic theoretical and computational study of two-equation models for near-wall turbulent flows was conducted. Two major problems established for the K-epsilon model are discussed, the lack of natural boundary conditions for the dissipation rate and the appearance of higher-order correlations in the balance of terms for the dissipation rate at the wall. The K-omega equation is shown to have two problems also: an exact viscous term is missing, and the destruction of the dissipation term is not properly damped near the wall. A new K-tau model (where tau = 1/omega is the turbulent time scale) was developed by inclusion of the exact viscous term, and by introduction of new wall damping functions with improved asymptotic behavior. A preliminary test of the new model yields improved predictions for the flat-plate turbulent boundary layer.

Speziale, Charles G.↗

DOE OSTI2023

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗

NASA NTRSDate not supplied

High-Lift Common Research Model: RANS, HRLES and WMLES Perspectives for CLmax Prediction Using LAVA

A unified assessment of three turbulence treatments: Reynolds Averaged Navier-Stokes(RANS), Hybrid RANS/LES (HRLES) and Equilibrium Wall-Modelled Large Eddy Simulation(WMLES) is presented for the High-Lift Common Research Model (CRM-HL). For the free-air configuration, steady-state RANS simulations show very accurate drag polar predictions in the low-𝛼linear regime. However, strong grid sensitivity is reported near the maximum lift-state(𝐶𝐿max), with finer-grids showing larger errors and predicting erroneous flow topologies on the wing. Our RANS simulations show that several corrections for the Spalart-Allmaras (SA)turbulence model widely used in the community lead to more erroneous results compared to the baseline closure, without exception. Both scale-resolving methods (HRLES and WMLES)address these drawbacks and predict an outboard separation pattern on the main element that is in good agreement with the oil flow photographs taken from the QinetiQ wind tunnel experiments, when LES-appropriate grids and numerical discretizations are used. While RANS simulations with the baseline SA closure do not show any wing-root separation post𝐶𝐿max, both HRLES and WMLES show onset of corner flow separation with varying degrees of progression, along with a weak pitch break in the wing-contribution of the overall pitching moment. This post-𝐶𝐿max pitch break seen in the free-air simulations is weaker than the break observed in experiments, with a weaker break reported in WMLES for each iteration of grid-refinement. In-tunnel simulations using both SA-baseline RANS and WMLES show a much stronger post-𝐶𝐿max break with the WMLES predictions showing excellent agreement with the experiment in terms of both the flow-topology observed and the pressure-coefficients at various spanwise stations. Sensitivity to the tunnel wall boundary layer is characterized via comparisons between viscous and inviscid treatments for the tunnel walls. WMLES predictions show moderate sensitivity at the predicted inboard flow-state at 𝐶𝐿max along with the progression towards a post-𝐶𝐿max stall; however, this stalled state at 𝛼≈20◦(inside the tunnel) obtained with both tunnel wall treatments appears to be largely identical.

TTT↗

DOE OSTI2022

Aero-Optics of Hypersonic Turbulent Boundary Layers

Aero-optics refers to optical distortions due to index-of-refraction gradients that are induced by aerodynamic density gradients. At hypersonic flow conditions, the bulk velocity is many times the speed of sound and density gradients may originate from shock waves, compressible turbulent structures, acoustic waves, thermal variations, etc. Due to the combination of these factors, aero-optic distortions are expected to differ from those common to sub-sonic and lower super-sonic speeds. This report summarizes the results from a 2019-2022 Laboratory Directed Research and Development (LDRD) project led by Sandia National Laboratories in collaboration with the University of Notre Dame, New Mexico State University, and the Georgia Institute of Technology. Efforts extended experimental and simulation methodologies for the study of turbulent hypersonic boundary layers. Notable experimental advancements include development of spectral de-aliasing techniques for highspeed wavefront measurements, a Spatially Selective Wavefront Sensor (SSWFS) technique, new experimental data at Mach 8 and 14, a Quadrature Fringe Imaging Interferometer (QFII) technique for time-resolved index-of-refraction measures, and application of QFII to shock-heated air. At the same time, model advancements include aero-optic analysis of several Direct Numerical Simulation (DNS) datasets from Mach 0.5 to 14 and development of wall-modeled Large Eddy Simulation (LES) techniques for aero-optic predictions. At Mach 8 measured and predicted root mean square Optical Path Differences agree within confidence bounds but are higher than semi-empirical trends extrapolated from lower Mach conditions. Overall, results show that aero-optic effects in the hypersonic flow regime are not simple extensions from prior knowledge at lower speeds and instead reflect the added complexity of compressible hypersonic flow physics.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Flow in additively manufactured super-rough channels

Metal additive manufacturing has enabled geometrically complex internal cooling channels for turbine and heat exchanger applications, but the process gives rise to large-scale roughness whose size is comparable to the channel height (which is 500 $\mathrm {\mu }$ m). These super-rough channels pose previously unseen challenges for experimental measurements, data interpretation and roughness modelling. First, it is not clear if measurements at a particular streamwise and spanwise location still provide accurate representation of the mean (time- and plane-averaged) flow. Second, we do not know if the logarithmic layer survives. Third, it is unknown how well previously developed rough-wall models work for these large-scale roughnesses. To answer the above practical questions, we conduct direct numerical simulations of flow in additively manufactured super-rough channels. Three rough surfaces are considered, all of which are obtained from computed tomography scans of additively manufactured surfaces. The roughness’ trough to peak sizes are 0.1 $h$ , 0.3 $h$ and 0.8 $h$ , respectively, where $h$ is the intended half-channel height. Each rough surface is placed opposite a smooth wall and the other two rough surfaces, leading to six rough-wall channel configurations. Two Reynolds numbers are considered, namely $Re_\tau =180$ and $Re_\tau =395$ . We show first that measurements at one streamwise and spanwise location are insufficient due to strong mean flow inhomogeneity across the entire channel, second that the logarithmic law of the wall survives despite the mean flow inhomogeneity and third that the established roughness sheltering model remains accurate.

Altland, Samuel↗

NASA NTRS1982

On the calculation of turbulent heat and mass transport downstream from an abrupt pipe expansion

A numerical study is reported of heat/mass transfer in the separated flow region created by an abrupt pipe expansion. Computations have employed a hybrid method of central and upwind finite differencing to solve the full Navier-Stokes equations with turbulent model (k approximately equal to epsilon). The study has given its main attention to the simulation of the region in the immediate vicinity of the wall, by formulating near-wall model for the evaluation of the mean generation and destruction rate of the epsilon equation. The computed results were compared with the experimental data and they showed generally encouraging agreement with the measurements.

Amano, R. S.↗