SEARCH · Engineering Papers

Results for “CFD simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

DOE OSTI2021

Forecasting Dynamic Line Rating with Spatial Variation Considerations

Dynamic line rating (DLR) is a technology that allows the ampacity of an electrical conductor to be calculated using real-time or forecasted weather conditions. Historically, the ampacity of a conductor has been determined using a static line rating method which assumes conservative weather assumptions. Therefore, not only can DLR give a more accurate measurement of the true ampacity of a conductor, but it can also increase its ampacity during weather conditions with greater thermal mitigations. The two primary cooling factors in the ampacity calculations are wind speed and direction. In complex terrain, wind speed and direction can have large variations over short distances. Therefore, accurately identifying the limiting span of a transmission line requires high spatial resolution of the wind along its path. One solution is to install dense weather stations along their path, though this can become costly over long distances. Therefore, researchers have investigated the use of Computation Fluid Dynamic (CFD) simulations to accurately compute the wind field along the path of a transmission line and use these results to identify the limiting section of the conductor. This work presents a case study that evaluates the coupling of CFD simulations and forecasted weather simulations using the High-Resolution Rapid Refresh (HRRR) model points over a 2-year span within a region in south eastern Idaho. The primary goal of the work is the evaluation of the number of HRRR model points used, i.e., weather stations, along the path of the line and the accuracy of the resulting ampacity. This was done using 4, 10, 17, 26, and 35 HRRR model points along two transmission line paths. The results indicate that as the number of model points are increased, the DLR ampacity of the lines decrease, yet converge as more points are added and demonstrate little change with additional HRRR points. It is expected that these results can help transmission line operators identify the number of weather stations that must be installed when coupled with CFD simulations and DLR ampacity to ensure accurate ratings and safe operations.

17 WIND ENERGY↗

DOE OSTI2021

Forecasting Dynamic Line Rating with Spatial Variation Considerations

Dynamic line rating (DLR) is a technology that allows the ampacity of an electrical conductor to be calculated using real-time or forecasted weather conditions. Historically, the ampacity of a conductor has been determined using a static line rating method which assumes conservative weather assumptions. Therefore, not only can DLR give a more accurate measurement of the true ampacity of a conductor, but it can also increase its ampacity during weather conditions with greater thermal mitigations. The two primary cooling factors in the ampacity calculations are wind speed and direction. In complex terrain, wind speed and direction can have large variations over short distances. Therefore, accurately identifying the limiting span of a transmission line requires high spatial resolution of the wind along its path. One solution is to install dense weather stations along their path, though this can become costly over long distances. Therefore, researchers have investigated the use of Computation Fluid Dynamic (CFD) simulations to accurately compute the wind field along the path of a transmission line and use these results to identify the limiting section of the conductor. This work presents a case study that evaluates the coupling of CFD simulations and forecasted weather simulations using the High-Resolution Rapid Refresh (HRRR) model points over a 2-year span within a region in south eastern Idaho. The primary goal of the work is the evaluation of the number of HRRR model points used, i.e., weather stations, along the path of the line and the accuracy of the resulting DLR ampacity. This was done using 4, 10, 17, 26, and 35 HRRR model points along two transmission line paths. The results indicate that as the number of model points are increased, the DLR ampacity of the lines decrease, yet converge as more points are added and demonstrate little change with additional HRRR points. It is expected that these results can help transmission line operators identify the number of weather stations that must be installed when coupled with CFD simulations and DLR ampacity to ensure accurate ratings and safe operations.

17 WIND ENERGY↗

NASA NTRS2009

Improvements in Modeling 90 degree Bleed Holes for Supersonic Inlets

The modeling of porous bleed regions as boundary conditions in computational fluid dynamics (CFD) simulations of supersonic inlet flows has been improved through a scaling of sonic flow coefficient data for 90deg bleed holes. The scaling removed the Mach number as a factor in computing the sonic flow coefficient and allowed the data to be fitted with a quadratic equation, with the only factor being the ratio of the plenum static pressure to the surface static pressure. The implementation of the bleed model into the Wind-US CFD flow solver was simplified by no longer requiring the evaluation of the flow properties at the boundary-layer edge. The quadratic equation can be extrapolated to allow the modeling of small amounts of blowing, which can exist when recirculation of the bleed flow occurs within the bleed region. The improved accuracy of the bleed model was demonstrated through CFD simulations of bleed regions on a flat plate in supersonic flow with and without an impinging oblique shock. The bleed model demonstrated good agreement with experimental data and three-dimensional CFD simulations of bleed holes.

Slater, John W.↗

DOE OSTI2020

Optimization of Window Positions for Wind-Driven Natural Ventilation Performance

This paper optimizes opening positions on building facades to maximize the natural ventilation’s potential for ventilation and cooling purposes. The paper demonstrates how to apply computational fluid dynamics (CFD) simulation results to architectural design processes, and how the CFD-driven decisions impact ventilation and cooling: (1) background: A CFD helps predict the natural ventilation’s potential, the integration of CFD results into design decision-making has not been actively practiced; (2) methods: Pressure data on building facades were obtained from CFD simulations and mapped into the 3D modeling environment, which were then used to identify optimal positions of two openings of a zone. The effect of the selected opening positions was validated with building energy simulations; (3) results: The cross-comparison study of different window positions based on different geographical locations quantified the impact on natural ventilation effectiveness; and (4) conclusions: The optimized window position was shown to be effective, and some optimal solutions contradicted the typical cross-ventilation strategy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

A Computational Investigation of Engine Heat Transfer with Ducted Fuel Injection

Ducted fuel injection (DFI) is an innovative method that curtails or prevents soot formation in direct-injection compression-ignition engines. DFI uses a simple duct, positioned outside each injector hole, facilitating the fuel/charge gas mixing before ignition. This reduces the equivalence ratio below two, in the autoignition zone, which in turn decreases soot formation. But this method also reduces fuel-conversion efficiency. This study investigates the effects of DFI on in-cylinder heat transfer. Experiments with conventional diesel combustion (CDC) and DFI were performed at four different dilution levels. Computational fluid dynamics (CFD) simulations were carried out at conditions matching those of the experiments, and the simulations were validated by the experimental data. The CFD simulations enabled to examine of in-cylinder heat release and temperature distributions. The heat transfer to the piston, head, and cylinder was investigated. The results show that DFI increased the heat transfer to the walls compared to CDC under the same conditions. This could help explain why DFI has been observed to reduce fuel-conversion efficiency by approximately 1% (absolute) relative to CDC under certain conditions. The efficiency loss typically decreases with dilution, such that DFI can improve fuel-conversion efficiencies relative to CDC at higher dilution levels.

Şener, Ramazan↗

NASA NTRS2018

Automated Boundary Conditions for Wind Tunnel Simulations

Computational fluid dynamic (CFD) simulations of models tested in wind tunnels require a high level of fidelity and accuracy particularly for the purposes of CFD validation efforts. Considerable effort is required to ensure the proper characterization of both the physical geometry of the wind tunnel and recreating the correct flow conditions inside the wind tunnel. The typical trial-and-error effort used for determining the boundary condition values for a particular tunnel configuration are time and computer resource intensive. This paper describes a method for calculating and updating the back pressure boundary condition in wind tunnel simulations by using a proportional-integral-derivative controller. The controller methodology and equations are discussed, and simulations using the controller to set a tunnel Mach number in the NASA Langley 14- by 22-Foot Subsonic Tunnel are demonstrated.

Carlson, Jan-Renee↗

DOE OSTI2021

Enabling Load Following Capability in the Transatomic Power MSR

This project is dedicated towards designing a fuel processing system that enables liquid-fueled molten salt reactors (MSR) to load follow by removing the dissolved xenon in the fuel salt. As one of the Gen-IV nuclear reactor concepts, the molten salt reactor receives increasing development interests in the recent years. One distinguishing feature of the liquid-fueled molten salt reactor is its improved ability to operate in a load-following mode by including the unique online fission product removal system. Load-following means that the reactor changes its power output based on the demand on the grid. Most of the current operating nuclear reactors have limited load-following ability and operate as the base load on the grid. Due to the rapid increase of solar energy, the requirement on load-following capacity is significantly increased because of the varying power output of the solar panels, yet the traditional load-following capacity is expected to decrease as the decarbonization of the grid continues. Therefore, the ability to perform load-following operation for the nuclear reactors will greatly enhance the resilience of the grid and make nuclear energy more economically competitive. This load-following feature is included in many commercial molten salt reactor designs, such as the designs by Transatomic Power, Terrestrial Energy, and Flibe Energy. Unfortunately, detailed analysis of the fuel processing system for commercial scale MSRs is still lacking, as well as how the fuel processing quantitively impacts the load-following operation. Moreover, experimental data for many of the underlying physics of fuel processing is limited. This project aims to pave the way for the fuel processing technology to advance to the commercial stage by performing combined experimental and simulation research. During the project period, four interconnected aspects of the development of the fuel processing system in liquid-fueled molten salt reactors are investigated. These aspects are the simulation and analysis of the fission product removal system, the fuel cycle simulation, the coupled reactor neutronics and thermal hydraulics transient simulation, and the gaseous fission product removal experiment. Multiphase CFD simulations are performed for components of the processing systems, and simplified air-water experiments are carried out to provide validation data. It is concluded that the CFD simulation can satisfactorily predict the system level performance of the components, and engineering models are constructed based on this success. Fuel cycle analysis is performed for two representative MSR design, the MSBR and the Transatomic Power MSR. Open-source code SaltProc is developed to incorporate the unique fuel processing system of the MSRs. It is concluded that the removal of xenon is essential for load-following operation in thermal spectrum MSR and Molten Salt Breeder Reactor. For the Transatomic Power MSR, the xenon poisoning effect is negligible due to its relatively fast neutron spectrum, though the overall fuel cycle economics still benefits from the removal of xenon. Coupled reactor neutronics and thermal hydraulics transient simulation is performed specifically for the Transatomic Power MSR. It is concluded that the reactor core design could perform power ramping fast enough to satisfy load-following operation. Combining the findings from each aspect, it is concluded that the load-following operation of a thermal neutron MSR is dependent upon the removal of xenon, which could be achieved for a commercial sized reactor using continuous inert gas sparging in a separate system with reasonable dimensions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2006

Simulations of SSLV Ascent and Debris Transport

A viewgraph presentation on Computational Fluid Dynamic (CFD) Simulation of Space Shuttle Launch Vehicle (SSLV) ascent and debris transport analysis is shown. The topics include: 1) CFD simulations of the Space Shuttle Launch Vehicle ascent; 2) Debris transport analysis; 3) Debris aerodynamic modeling; and 4) Other applications.

Rogers, Stuart↗

DOE OSTI2024

Performance Testing of a Moving-Bed Gasifier Using Coal, Biomass, and Waste Plastic Blends with Washed and Unwashed Legacy Coals and Other Waste Fuels to Generate White Hydrogen

The objective of this effort, primarily funded by the United States Department of Energy (DOE), and led by the Electric Power Research Institute, Inc. (EPRI), with support by Hamilton Maurer International (HMI) and Sotacarbo S.p.A. (Sotacarbo), has been to qualify coal, biomass, and plastic waste blends based on performance testing of selected fuel pellet compositions in a pilot-scale updraft moving-bed (UDMB) gasifier. The testing provided relevant data to advance the commercial-scale design of the moving-bed gasifier to be able to successfully use these feedstocks to produce hydrogen. In particular, the effects of waste plastics on feedstock development (i.e., blending and pelletizing) and the resulting products (i.e., syngas compositions, organic condensate production, and ash characteristics) are the focus. The gasifier used for testing is HMI’s moving-bed gasifier, which has been proven capable of gasifying nearly all coal ranks. It has also shown the ability in prior testing work to gasify wood chips (biomass). However, mixtures of these fuels with plastic wastes have not been prepared and gasified together. The three feedstocks were densified and pelletized by California Pellet Mill (CPM) to meet the feedstock size required by Sotacarbo’s 30mm ID UDMB gasifier, under contract to HMI. The technical tasks and results from this two-year research project included: (1) Feed Procurement and Preparation: Nine different tri-fuel pellets were prepared from varying compositions of fresh mined PRB coal, corn stover biomass, and car fluff waste plastics. Tri-fuel pellets were produced by CPM and shipped to Sotacarbo’s test facility in Carbonia, Sardinia, Italy. (2) Test Plan Development: A test plan was created to define the test runs to be performed. The test plan detailed the different UDMB gasification tests to be performed in Sotacarbo’s 12-inch ID pilot scale gasifier, the process monitoring instrumentation used, and the extractive samples recovered for analysis of the total gasification process mass and energy balance. (3) Gasifier Testing: Nine different gasification runs were performed in the pilot-scale gasifier at Sotacarbo using nine different fuel feedstock compositions generated from varying mixtures of PRB coal, biomass, and plastic wastes. The testing generated performance data on gasification reaction efficiency and performance, yielding relevant data for models used to scale up the gasifier design. This task also included work to refurbish and reassemble the pilot gasifier at Sotacarbo and perform a baseline 100% PRB coal run. (4) Data Analysis and Reporting: Review of the data, determination of figures of merit, and interpretation of the results are reported in the project’s final report, published in March 2024. The results show that all tri-fuel pellets gasified well and maintained structural integrity throughout the gasification process. The syngas generated can be shifted to hydrogen by using commercial syngas shifting technologies. (5) High Fidelity computational fluid dynamics (CFD) Simulation: The National Energy Technology Laboratory (NETL) team performed CFD simulations of the UDMB gasifier for two of the tri-fuel pellets gasified in Sotacarbo’s pilot scale gasifier. The kinetic mechanisms for the pyrolysis of each constituent, PRB coal, corn stover biomass, and waste plastics are based on thermogravimetric analysis performed by Sotacarbo. The gasification model was validated by comparing the predicted syngas composition at the exit of the gasifier with the measured syngas composition. In addition, the reactor’s measured internal temperature profile agreed well with the predicted internal reactor temperature profile. These results validate that the model can be used to predict the performance of the updraft moving bed gasifier for different feedstocks and operating conditions. This paper summarizes the results of the completed work in which the pelletizing procedure was validated to ensure the viability of the tri-fuel pellets for the gasification runs performed at Sotacarbo’s 30 mm UDMB gasifier. The gasification performance data from this series of nine runs will enable modeling of a full-scale HMI industrial scale gasifier supporting both combined heat and power, and Hydrogen production from coal (both fresh mined and legacy) combined with various biomass and waste plastics. Additionally, plans and progress on a follow-up project, being executed by the same project team, will be presented. In this project, a total of twenty (20) different feedstocks are being prepared from varying compositions of biomass (both woody biomass and corn stover) with a mixture of legacy coal waste, plastic waste, and refuse-derived fuel (RDF). The testing will provide information on gasification reaction efficiency/performance, yielding relevant data for models used to scale up the gasifier design to 50 megawatt electric (MWe) (equivalent hydrogen production). Tests will also be performed on a bench-scale fluidized-bed gasifier for comparison purposes. The results of this testing will be used to specify the range of feedstock blends that can be successfully gasified as well as quantify gasifier outputs based on specific blends.

08 HYDROGEN↗

NASA NTRS2017

Validation Process for LEWICE by Use of a Navier-Stokes Solver

A research project is underway at NASA Glenn to produce computer software that can accurately predict ice growth under any meteorological conditions for any aircraft surface. This report will present results from the latest LEWICE release, version 3.5. This program differs from previous releases in its ability to model mixed phase and ice crystal conditions such as those encountered inside an engine. It also has expanded capability to use structured grids and a new capability to use results from unstructured grid flow solvers. A quantitative comparison of the results against a database of ice shapes that have been generated in the NASA Glenn Icing Research Tunnel (IRT) has also been performed. This paper will extend the comparison of ice shapes between LEWICE 3.5 and experimental data from a previous paper. Comparisons of lift and drag are made between experimentally collected data from experimentally obtained ice shapes and simulated (CFD) data on simulated (LEWICE) ice shapes. Comparisons are also made between experimentally collected and simulated performance data on select experimental ice shapes to ensure the CFD solver, FUN3D, is valid within the flight regime. The results show that the predicted results are within the accuracy limits of the experimental data for the majority of cases.

large-eddy simulations↗

NASA NTRS2017

A Revised Validation Process for Ice Accretion Codes

turbulent heat flux↗

DOE OSTI2021

MFIX-Exa: A path toward exascale CFD-DEM simulations

MFIX-Exa is a computational fluid dynamics–discrete element model (CFD-DEM) code designed to run efficiently on current and next-generation supercomputing architectures. MFIX-Exa combines the CFD-DEM expertise embodied in the MFIX code—which was developed at NETL and is used widely in academia and industry—with the modern software framework, AMReX, developed at LBNL. The fundamental physics models follow those of the original MFIX, but the combination of new algorithmic approaches and a new software infrastructure will enable MFIX-Exa to leverage future exascale machines to optimize the modeling and design of multiphase chemical reactors.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

OC7 phase I: CFD investigation of viscous forces on rectangular members of semisubmersibles

In Phase I of the Offshore Code Comparison Collaboration 7 (OC7) project, we aim to address key questions related to the modeling of viscous hydrodynamic loads on floating platforms. Here, this paper presents results from a collaborative effort within Phase I of the OC7 project. We investigate viscous forces acting on the rectangular members of semisubmersible platforms using high-fidelity computational fluid dynamics (CFD) simulations. The parametric CFD analyses include two-and-a-half-dimensional (2.5D) simulations of pontoon sections of two different semisubmersible platform designs (the VolturnUS-S and the INO WINDMOOR) and full three-dimensional (3D) structure simulations of the VolturnUS-S platform. We extract both the added mass and drag coefficients that are typically used in engineering-level hydrodynamic models from the 2.5D forced oscillation simulations with various Keulegan-Carpenter numbers and compare them to the coefficients from the 3D forced oscillation simulations. We observe consistent variation of the coefficients along the pontoon legs due to 3D platform and free-surface effects. With the breadth of the simulations and component-level detailed analyses, this work establishes a useful database and several empirical formulas for selecting hydrodynamic coefficients when modeling semisubmersible platforms in engineering tools.

13 HYDRO ENERGY↗

NASA NTRSDate not supplied

Pretest Simulations of a Supersonic Mixing and Combustion Validation Experiment to Assess Sensitivities

The reliance on CFD simulations to develop, design, and optimize scramjet systems (or components) has become commonplace. This reliance inevitably hinges on the ability of the computational analyst to quantify the level of confidence in their computational results. Unfortunately, nearly all the measured data available for this assessment comes from antiquated experimental datasets, or from tests that focused on the extraction of scramjet system (or component) performance. The objective of a CFD validation experiment is to quantify the predictive accuracy of one or more of the CFD physics submodels, implying that other uncertainties related to replicating the facility flow environment (e.g., knowledge of boundary conditions) must be minimized to the extent possible. This inevitably places stringent requirements on the quality and quantity of measurements taken to accurately specify inflow, outflow, and surface conditions for the CFD simulations; in addition to the measurements taken for the validation of physics submodels. This places additional demands on the experimental process above and beyond those for test article performance assessment. A recent high speed code credibility workshop series sponsored by AFRL identified a gap in existing validation data for fundamental assessments of turbulent mixing and combustion CFD closure models at scramjet engine relevant conditions. To address this gap, engineers at AFRL have designed a coaxial jet flame configuration that will be tested at two facilities (Research Cell 19 at the Air Force Research Lab, and at Purdue University). The effort described here documents pretest simulations of this validation experiment with the goal of fleshing out the extent of the facility flowpath that must be included to adequately reproduce the facility test section flow environment. The findings indicate that the flow around the support structure for the fuel injection centerbody upstream of the facility nozzle generates disturbances that persist throughout the nozzle expansion process; corrupting the azimuthal symmetry that was desired in the fuel/air mixing region of the test section. Simulations without this support structure maintained a high degree of azimuthal symmetry up until the fuel injection plane. However, even in this scenario the azimuthal symmetry was not maintained once the centerbody boundary layer transitioned to a wake flow downstream of the fuel injection plane.

CFD↗

NASA NTRSDate not supplied

Pretest Simulations of a Supersonic Mixing and Combustion Validation Experiment to Assess Sensitivities

CFD↗

NASA NTRS2018

Generation of Performance Model for the Aeolian Wind Tunnel (AWT) Rotor at Reduced Pressure

The NASA Jet Propulsion Laboratory (JPL) designed the Mars Helicopter (MH) in collaboration with AeroVironment Inc., NASA Ames Research Center, and NASA Langley Research Center to explore the possibility of a vertical takeoff and landing (VTOL) Unmanned Aerial Vehicle (UAV) for flight on Mars. A 40-inch-diameter Aeolian Wind Tunnel (AWT) rotor, roughly approximating the proposed MH design by JPL, was tested in forward flight at Mars atmospheric pressure at the NASA Ames Planetary Aeolian Laboratory (PAL) in support of MH research efforts. This report describes the generation of the rotor model used to correlate with that experimental effort as reported by Ament and Koning. The 40-inch-diameter rotor was 3D-scanned and transformed into an airfoil deck. The scanned rotor airfoil sections are analyzed using C81 Generator (C81Gen) to generate the sectional aerodynamic coefficients for comprehensive analyses. A mid-fidelity computational fluid dynamics (CFD) simulation using Rotorcraft CFD (RotCFD) is pursued to efficiently estimate rotor hover and forward flight performance. Simulations at two pressures, 7 mbar (approximate Martian atmospheric pressure) and 1018 mbar (1 atmosphere), are performed to gain an understanding of the performance differences and Reynolds number effects observed. Experimental 1-atmosphere thrust for single- and dual-rotor isolated hover cases correlate well with the modeled rotor. Performance results at reduced pressure (7 mbar) show a drastic decrease in lift for equivalent RPMs tested at 1 atmosphere. Although this is primarily due to pressure reduction, Reynolds number effects also contribute to this decrease, as airfoil lift and drag coefficients are affected when compared with 1-atmosphere results. Further, simulated rotor power coefficient shows drastic increases at reduced pressures, attributed to laminar boundary layer separation, as described in Koning et al. for the MH rotor analysis. PAL experimental Martian Surface Wind Tunnel (MARSWIT) results are presented in the paper by Ament and Koning. The very low Reynolds number range is currently not well understood and presents various challenges for both experimentation and simulation.

Koning, Witold J. F.↗

NASA NTRSDate not supplied

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗