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 595 records · Page 33

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗

Néel domain walls with bistable chirality in a perpendicularly magnetized ferrimagnetic insulator

Field-free spin-orbit torque-driven domain wall motion in magnetic thin films with perpendicular magnetic anisotropy (PMA) requires the domain walls to have Néel character. Conventionally, Néel domain walls are stabilized by the Dzyaloshinskii-Moriya interaction (DMI) in ultrathin films. Here, in a europium iron garnet thin film with PMA and an additional uniaxial in-plane anisotropy, we demonstrate two bistable Néel domain wall states in the absence of DMI, and the capability to toggle the wall states with an in-plane field pulse and consequently their directions of motion under a current pulse. We present a phase diagram for the bistable Néel domain wall states as a function of in-plane field pulse width and amplitude. By fitting the experimental data to an analytical model of Néel wall reversal through the nucleation and propagation of Bloch lines, we extract the length of the initial reversed domain wall segment and Bloch line nucleation energy barrier. Current-driven motion of in-plane anisotropy stabilized Néel walls is qualitatively different from that of DMI-stabilized ones owing to the different symmetry of the effective fields that stabilize the Néel configuration. Furthermore, we present a proof of principle demonstration for 2-bit random number generation based on the stochastic reversal of domain wall chirality. These results provide critical insight into the topological energy barrier of Bloch lines and identify paths towards domain wall-based memory and computing devices.

36 MATERIALS SCIENCE↗

Measurements and Computations of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

The NASA Juncture Flow experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction specifically for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This paper presents and discusses the results of a recent experiment with the juncture-flow model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of the test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model angle of incidence. Laser Doppler velocimetry (LDV) measurements were made at three model angles of incidence (0 deg: fully attached, 1 deg: incipient separation, and 5 deg: separated flow) and for each one, mean-flow and Reynolds-stress data were obtained on the fuselage and at several streamwise locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and included model and tunnel wall static pressures, tunnel wall and ceiling boundary-layer rake data, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. Comparisons between the experimental data on the test article and Reynolds-averaged Navier-Stokes CFD results are presented and discussed.

Juncture Flow↗

Measurements and Computations of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

The NASA Juncture Flow experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction specifically for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This paper will present and discuss the results of a recent experiment with the juncture-flow model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of that test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model pitch angle. Laser Doppler velocimetry (LDV) measurements were made at three model pitch angles (0 deg: fully attached, 1 deg: incipient separation, and 5 deg: separated flow) and for each one, mean-flow and Reynolds-stress data was obtained on the fuselage and at several streamwise locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and included model and tunnel wall static pressures, tunnel wall and ceiling boundary-layer rake data, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. Comparisons between the experimental data and Reynolds averaged Navier-Stokes CFD results will also be presented and discussed.

Juncture Flow↗

Measurements and Computations of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

The NASA Juncture Flow experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction specifically for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This paper will present and discuss the results of a recent experiment with the juncture-flow model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of that test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model pitch angle. Laser Doppler velocimetry (LDV) measurements were made at three model pitch angles (0 deg: fully attached, 1 deg: incipient separation, and 5 deg: separated flow) and for each one, mean-flow and Reynolds-stress data was obtained on the fuselage and at several streamwise locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and included model and tunnel wall static pressures, tunnel wall and ceiling boundary-layer rake data, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. Comparisons between the experimental data and Reynolds averaged Navier-Stokes CFD results will also be presented and discussed.

Juncture Flow↗

Effect of blockage ratio on drag and pressure distributions for bodies of revolution at transonic speeds

Experimental data were obtained in two wind tunnels for 13 models over a Mach number range from 0.70 to 1.02. Effects of increasing test-section blockage ratio in the transonic region near a Mach number of 1.0 included change in the shape of the drag curves, premature drag creep, delayed drag divergence, and a positive increment of pressures on the model afterbodies. Effects of wall interference were apparent in the data even for a change in blockage ratio from a very low 0.000343 to an even lower 0.000170. Therefore, models having values of blockage ratio of 0.0003 - an order of magnitude below the previously considered safe value of 0.0050 - had significant errors in the drag-coefficient values obtained at speeds near a Mach number of 1.0. Furthermore, the flow relief afforded by slots or perforations in test-section walls - designed according to previously accepted criteria for interference-free subsonic flow - does not appear to be sufficient to avoid significant interference of the walls with the model flow field for Mach numbers very close to 1.0.

Couch, L. M.↗

Spacecraft Internal Acoustic Environment Modeling

Acoustic modeling can be used to identify key noise sources, determine/analyze sub-allocated requirements, keep track of the accumulation of minor noise sources, and to predict vehicle noise levels at various stages in vehicle development, first with estimates of noise sources, later with experimental data. This paper describes the implementation of acoustic modeling for design purposes by incrementally increasing model fidelity and validating the accuracy of the model while predicting the noise of sources under various conditions. During FY 07, a simple-geometry Statistical Energy Analysis (SEA) model was developed and validated using a physical mockup and acoustic measurements. A process for modeling the effects of absorptive wall treatments and the resulting reverberation environment were developed. During FY 08, a model with more complex and representative geometry of the Orion Crew Module (CM) interior was built, and noise predictions based on input noise sources were made. A corresponding physical mockup was also built. Measurements were made inside this mockup, and comparisons were made with the model and showed excellent agreement. During FY 09, the fidelity of the mockup and corresponding model were increased incrementally by including a simple ventilation system. The airborne noise contribution of the fans was measured using a sound intensity technique, since the sound power levels were not known beforehand. This is opposed to earlier studies where Reference Sound Sources (RSS) with known sound power level were used. Comparisons of the modeling result with the measurements in the mockup showed excellent results. During FY 10, the fidelity of the mockup and the model were further increased by including an ECLSS (Environmental Control and Life Support System) wall, associated closeout panels, and the gap between ECLSS wall and mockup wall. The effect of sealing the gap and adding sound absorptive treatment to ECLSS wall were also modeled and validated.

Chu, Shao-Sheng R.↗

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Turbulence closure modeling near rigid boundaries

The near-wall region plays an essential role in turbulent boundary layers: it is a region of high shear; the peak rate of production and peak intensity of turbulence occurs there; and the peak rate of dissipation occurs right at the wall. Nevertheless, this region has received less attention from modelers than have more nearly homogeneous flows. One reason for this is that when the boundary layer is near equilibrium, experimental data can be used to prescribe the flow in the wall layer. Another reason is that most turbulence models are developed under assumptions of near homogeneity. This is a poor approximation in the wall region. A single-point moment closure model for the strongly non-homogeneous A turbulent flow near a rigid boundary is developed.

Durbin, Paul A.↗

Progress Towards Fuselage Drag Reduction via Active Flow Control: A Combined CFD and Experimental Effort

A combined computational and experimental effort has been undertaken to study fuselage drag reduction on a generic, non-proprietary rotorcraft fuselage by the application of active ow control. Fuselage drag reduction is an area of research interest to both the United States and France and this area is being worked collaboratively as a task under the United States/France Memorandum of Agreement on Helicopter Aeromechanics. In the first half of this task, emphasis is placed on the US generic fuselage, the ROBIN-mod7, with the experimental work being conducted on the US side and complementary US and French CFD analysis of the baseline and controlled cases. Fuselage simulations were made using Reynolds-averaged Navier-Stokes ow solvers and with multiple turbulence models. Comparisons were made to experimental data for numerical simulations of the isolated fuselage and for the fuselage as installed in the tunnel, which includes modeling of the tunnel contraction, walls, and support fairing. The numerical simulations show that comparisons to the experimental data are in good agreement when the tunnel and model support are included. The isolated fuselage simulations compare well to each other, however, there is a positive shift in the centerline pressure when compared to the experiment. The computed flow separation locations on the rear ramp region had only slight differences with and without the tunnel walls and model support. For the simulations, the flow control slots were placed at several locations around the flow separation lines as a series of eight slots that formed a nearly continuous U-shape. Results from the numerical simulations resulted in an estimated 35% fuselage drag reduction from a steady blowing flow control configuration and a 26% drag reduction for unsteady zero-net-mass flow control configuration. Simulations with steady blowing show a delayed flow separation at the rear ramp of the fuselage that increases the surface pressure acting on the ramp, thus decreasing the overall fuselage pressure drag.

Schaeffler, Norman W.↗

Space Shuttle main engine nozzle-steerhorn dynamics

On two occasions during the Space Shuttle main engine development, the LH2 feedline (called the steerhorn, because of its shape) failed during the cutoff transient. A dynamic test was undertaken, and an analytical model was developed and correlated to the dynamic test. Detailed models of the tube bundle were required to obtain the equivalent shell coefficients. All-shell models of the nozzle wall were found better than beam-shell models. The most difficult part of the structure to simulate was the felt-metal pad between the feedline and its mount, which introduced nonlinear stiffness and damping and led to the use of separate low amplitude and high amplitude models. The total structure was found to have 400 modes in the frequency range of interest, 0 to 500 Hz. Good test analysis correlation was obtained and a modified feedline configuration was found to demonstrate a 40% reduction of response stress from the original configuration.

Kiefling, L.↗

The present state and future direction of second order closure models for compressible flows

The topics are presented in viewgraph form and include: (1) Reynolds stress closure models; (2) Favre averages and governing equations; (3) the model for the deviatoric part of the pressure-strain rate correlation; (4) the SSG pressure-strain correlation model; (5) a compressible turbulent dissipation rate model; (6) variable viscosity effects; (7) near-wall stiffness problems; (8) models of the Reynolds mass and heat flux; and (9) a numerical solution of the compressible turbulent transport equation.

Gatski, Thomas B.↗

Heat transfer in rotating serpentine passages with trips skewed to the flow

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass, heat transfer model with both radially inward and outward flow. Trip strips, skewed at 45 deg to the flow direction, were machined on the leading and trailing surfaces of the radial coolant passages. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature, rotation number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from similar stationary and rotating models with smooth walls and with trip strips normal to the flow direction. The heat transfer coefficients on surfaces, where the heat transfer decreased with rotation and buoyancy, decreased to as low as 40 percent of the value without rotation. However, the maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels previously obtained with the smooth wall models. It was concluded that (1) both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips, (2) the effects of rotation are markedly different depending upon the flow direction, and (3) the heat transfer with skewed trip strips is less sensitive to buoyancy than the heat transfer in models with either smooth or normal trips. Therefore, skewed trip strips rather than normal trip strips are recommended and geometry-specific tests are required for accurate design information.

Johnson, B. V.↗

Empower Wall: Active insulation system leveraging additive manufacturing and model predictive control

Buildings are one of the largest energy consumers worldwide, using large amounts of energy during their construction and for climate control during operation. Active insulation systems (AIS) have been shown to reduce the energy needed for climate control in buildings by dynamically regulating the heat transferred between a building’s interior and exterior. Infrastructure-scale additive manufacturing (AM) has the potential to reduce the resources needed for building construction. Combining these two technologies into a single building envelope would create a path towards more sustainable buildings. A test was conducted for the Federal Energy Management Program (FEMP) Energy Exchange training and trade show, in August 2021, to investigate a new building envelope design, termed the Empower Wall, that utilized an AIS and was constructed using AM. Model predictive control was implemented to manage operation of the Empower Wall in concert with the existing HVAC system. Finally, the prototype system demonstrated that the Empower Wall lowered total energy consumption and reduced the cost of energy used.

36 MATERIALS SCIENCE↗

Internal energy balance and aerodynamic heating predictions for hypersonic turbulent boundary layers

The elemental equation governing heat transfer in aerodynamic flows is the internal energy equation. For a boundary layer flow, a double integration of the Reynolds-averaged form of this equation provides an expression of the wall heat flux in terms of the integrated effects, over the boundary layer, of various physical processes: turbulent dissipation, mean dissipation, turbulent heat flux, etc. Recently available direct numerical simulation data for a Mach 11 cold-wall turbulent boundary layer allows a comparison of the exact contributions of these terms in the energy equation to the wall heat flux with their counterparts modeled in the Reynolds-averaged Navier-Stokes (RANS) framework. Various approximations involved in RANS, both closure models as well as approximations involved in adapting incompressible RANS models to a compressible form, are assessed through examination of the internal energy balance. There are a number of potentially problematic assumptions and terms identified through this analysis. Here, the effect of compressibility corrections of the dilatational dissipation type is explored, as is the role of the modeled turbulent dissipation, in the context of wall heat flux predictions. The results indicate several potential avenues for RANS model improvement for hypersonic cold-wall boundary-layer flows.

42 ENGINEERING↗

Discontinuous Galerkin Methods for Turbulence Simulation

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Collis, S. Scott↗

Efficient Parametric Uncertainty Analysis of an Earth Entry Vehicle Concept Using Least Angle Regression

The objective of this work was to outline and apply an efficient and accurate parametric un-certainty propagation approach to the analysis of convective heating on an Earth entry vehicle concept. The described approach was based on Least Angle Regression used to solve a sparse and underdetermined linear system in the point-collocation non-intrusive polynomial chaos surrogate method. This approach involved an iterative process to computing the non-zero terms of the underlying polynomial chaos model using only enough samples to converge uncertainty interval predictions and Sobol index values based global nonlinear sensitivity estimates. The Earth entry vehicle was analyzed at three points along a representative trajectory for a Mars return mission. 329 sources of uncertainty were identified in the computational fluid dynamics model used to predict the forebody convective heating. These included uncertainty in flow field chemical rates, collision integrals, heats of formation, surface finite rate char model reaction rates, wall roughness height, and the turbulent Schmidt number. Results from this study showed that convective heating uncertainty as high as 50% of the nominal was predicted with only about 50 evaluations of the computational model. This was far fewer than would be required for a sampling-based approach or a full basis polynomial chaos model, which would have required over 50,000 samples. Additionally, results showed that over 90% of the total convective heating uncertainty was due to uncertainty in the N2catalytic rate on the surface, while the remainder of the uncertainty was attributed to the turbulent Schmidt number and the wall roughness uncertainties.

Thomas K West IV↗

Heat transfer in rotating serpentine passages with smooth walls

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass, smooth-wall heat transfer model with both radially inward and outward flow. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages (coolant-to-wall temperature ratio, Rossby number, Reynolds number and radius-to-passage hydraulic diameter ratio). These four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. It was found that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs and that the effect of rotation on the heat transfer coefficients was markedly different depending on the flow direction. Local heat transfer coefficients were found to decrease by as much as 60 percent and increase by 250 percent from no rotation levels. Comparisons with a pioneering stationary vertical tube buoyancy experiment showed reasonably good agreement. Correlation of the data is achieved employing dimensionless parameters derived from the governing flow equations.

Wagner, J. H.↗