Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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 217 records · Page 12

Computational Fluid Dynamics Analysis for the Orbiter LH2 Feedline Flowliner

In phase II, additional inducer rotations are simulated in order to understand the root cause of the flowliner crack problem. CFD results confirmed that there is a strong unsteady interaction between the backflow regions caused by the LPFTP inducer and secondary flow regions in the bellows cavity through the flowliner slots. It is observed that the swirl on the duct side of the downstream flowliner is stronger than on the duct side of the upstream flowliner. Due to this swirl, there are more significant unsteady flow interactions through the downstream slots than those observed in the upstream slots. Averaged values of the local velocity at the slots were provided to the NESC-ITA flow physics acoustics team to guide them in designing the acoustics experiment. A parametric study was performed to compare the flow field in the flowliner area when one upstream slot and one corresponding downstream slot were enlarged. No significant differences were observed between the flow field obtained from the enlarged slot configuration when compared with the original configuration. More cases must be analyzed with various enlarged slot configurations to generalize this observation. The flow through the A1 test stand and the flow through the orbiter fuel feedline manifold were simulated without the LPFTP. It was observed that incoming flow to the flowliner and inducer was more uniform in the A1 test stand then in the orbiter manifold. Additionally, each engine LPFTP in the orbiter receives significantly different velocity distributions. Because of the differences observed in the computed results, it is not possible for the A1 test stand to represent the three different engine feedlines simultaneously.

Kiris, Cetin C.↗

Three-Dimensional Structure and Modeling of a Normal Bifurcated Shock Train from Experimental Measurements

An asymmetric, normal bifurcated shock train in a Mach 2 constant area, rectangular duct is investigated in the University of Michigan Direct Connect Isolator facility. High-speed schlieren imaging, wall static pressure measurements, surface oil flow visualization, and particle image velocimetry of the shock train are synthesized into a three-dimensional representation. This visualization is then used to inform the underlying flow physics of the distributed fluid dynamical processes along the shock train. The detailed, three-dimensional morphology of the flow profile entering the shock train is shown to have a significant impact on the separated flow morphology of the shock train. This results in a skewed supersonic core flow with alternating separation bubbles in the surrounding boundary layers. Finally, the efficacy of the prevailing pseudo-shock models in the literature is analyzed with the support of the available flow measurements. None are found to model both the pressure and streamwise Mach number profile of the UMDCI pseudo-shock. Note: This presentation is accompanied by the video presentation that can be downloaded for viewing and is formatted as an mp4 with a runtime of 10 min 15 secs.

Hypersonics↗

Detailed Velocity, Temperature, and Heat Flux Measurements on a Large Scale Film Cooling Model

To investigate the flow physics of turbine film cooling, detailed mean and fluctuating velocity and temperature surveys were made in the turbulent jet flow field behind a row of large scale film cooling holes. Additionally, surface heat transfer, film effectiveness, and turbulent heat flux values were determined. Measurements were made in the NASA Engine Research Building (ERB) SW-6 test facility, which consists of an 8.15-in2 open-inlet wind tunnel with a temperature-controlled secondary flow system. The film-cooling test plate fit into the wind tunnel such that the heat transfer surface formed the wind tunnel floor. This plate consisted of a three-hole array of film cooling holes that were fed from a plenum of different temperature.

Fluid Mechanics↗

Forebody and Inlet Design for the HIFiRE 2 Flight Test

A forebody and inlet have been designed for the HIFiRE 2 scramjet flight test. The test will explore the operating, performance, and stability characteristics of a simple hydrocarbon-fueled scramjet combustor as it transitions from dual-mode to scramjet-mode operation and during supersonic combustion at Mach 8+ flight conditions. Requirements for the compression system were derived from inlet starting and combustor inflow requirements as well as physical size constraints. The design process is described. A planar, fixed geometry, mixed compression concept was used to produce laterally uniform flow at the inlet entrance and a conservative amount of internal contraction with respect to inlet starting. A grid sensitivity study was performed so that important flow physics caused by three-dimensional shock boundary layer interactions could be captured with confidence. Results from low Mach number operability studies, nominal trajectory cases, and high dynamic pressure heat load cases are discussed. The forebody and inlet solutions provide information for on-going combustor calculations, mass capture across the trajectory for fuel system design, and surface heating rates for thermal/structural analysis. The design has a one freestream Mach number margin for inlet starting, exceeds the high Mach number combustor entrance pressure requirement, produces high quality flow at the inlet exit for all Mach numbers and vehicle attitudes in the design space, and fits inside the booster shroud.

Ferlemann, Paul G.↗

CFD for hypersonic airbreathing aircraft

A general discussion is given on the use of advanced computational fluid dynamics (CFD) in analyzing the hypersonic flow field around an airbreathing aircraft. Unique features of the hypersonic flow physics are presented and an assessment is given of the current algorithms in terms of their capability to model hypersonic flows. Several examples of advanced CFD applications are then presented.

Kumar, Ajay↗

CFD assessment of the pollutant environment from RD-170 propulsion system testing

Computational Fluid Dynamics (CFD) technology has been used to assess the exhaust plume pollutant environment of the RD-170 engine hot-firing on the F1 Test Stand at Marshall Space Flight Center. Researchers know that rocket engine hot-firing has the potential for forming thermal nitric oxides (NO(x)), as well as producing carbon monoxide (CO) when hydrocarbon fuels are used. Because of the complicated physics involved, however, little attempt has been made to predict the pollutant emissions from ground-based engine testing, except for simplified methods which can grossly underpredict and/or overpredict the pollutant formations in a test environment. The objective of this work, therefore, has been to develop a technology using CFD to describe the underlying pollutant emission physics from ground-based rocket engine testing. This resultant technology is based on a three-dimensional (3D), viscous flow, pressure-based CFD formulation, where wet CO and thermal NO finite-rate chemistry mechanisms are solved with a Penalty Function method. A nominal hot-firing of a RD-170 engine on the F1 stand has been computed. Pertinent test stand flow physics such as the multiple-nozzle clustered engine plume interaction, air aspiration from base and aspirator, plume mixing with entrained air that resulted in contaminant dilution and afterburning, counter-afterburning due to flame bucket water-quenching, plume impingement on the flame bucket, and restricted multiple-plume expansion and turning have been captured. The predicted total emission rates compared reasonably well with those of the existing hydrocarbon engine hot-firing test data.

Wang, Ten-See↗

Aeroheating Design Issues for Reusable Launch Vehicles: A Perspective

An overview of basic aeroheating design issues for Reusable Launch Vehicles (RLV), which addresses the application of hypersonic ground-based testing, and computational fluid dynamic (CFD) and engineering codes, is presented. Challenges inherent to the prediction of aeroheating environments required for the successful design of the RLV Thermal Protection System (TPS) are discussed in conjunction with the importance of employing appropriate experimental/computational tools. The impact of the information garnered by using these tools in the resulting analyses, ultimately enhancing the RLV TPS design is illustrated. A wide range of topics is presented in this overview; e.g. the impact of flow physics issues such as boundary-layer transition, including effects of distributed and discrete roughness, shock-shock interactions, and flow separation/reattachment. Also, the benefit of integrating experimental and computational studies to gain an improved understanding of flow phenomena is illustrated. From computational studies, the effect of low-density conditions and of uncertainties in material surface properties on the computed heating rates a r e highlighted as well as the significant role of CFD in improving the Outer Mold Line (OML) definition to reduce aeroheating while maintaining aerodynamic performance. Appropriate selection of the TPS design trajectories and trajectory shaping to mitigate aeroheating levels and loads are discussed. Lastly, an illustration of an aeroheating design process is presented whereby data from hypersonic wind-tunnel tests are integrated with predictions from CFD codes and engineering methods to provide heating environments along an entry trajectory as required for TPS design.

Zoby, E. Vincent↗

Adaptive refinement tools for tetrahedral unstructured grids

An exemplary embodiment providing one or more improvements includes software which is robust, efficient, and has a very fast run time for user directed grid enrichment and flow solution adaptive grid refinement. All user selectable options (e.g., the choice of functions, the choice of thresholds, etc.), other than a pre-marked cell list, can be entered on the command line. The ease of application is an asset for flow physics research and preliminary design CFD analysis where fast grid modification is often needed to deal with unanticipated development of flow details.

Pao, S. Paul↗

Unsteady Loss in the Stator Due to the Incoming Rotor Wake in a Highly-Loaded Transonic Compressor

The present paper reports an investigation of unsteady loss generation in the stator due to the incoming rotor wake in an advanced GE transonic compressor design with a high-fidelity numerical method. This advanced compressor with high reaction and high stage loading has been investigated both experimentally and analytically in the past. The measured efficiency in this advanced compressor is significantly lower than the design intention/goal. The general understanding is that the current generation of compressor design/analysis tools miss some important flow physics in this modern compressor design. To pinpoint the source of the efficiency miss, an advanced test with a detailed flow traverse was performed for the front one and a half stage at the NASA Glenn Research Center. Detailed data-match analysis by GE identified an unexpected high loss generation in the pressure side of the stator passage. Higher total temperature and lower total pressure are measured near the pressure side of the stator. Various analyses based on the RANS and URANS of the compressor stage do not calculate the measured higher total temperature and lower total pressure on the pressure side of the stator. In the present paper, a Large Eddy Simulation (LES) is applied to find the fundamental mechanism of this unsteady loss generation in the stator due to the incoming rotor wake. The results from the LES were first compared with the NASA test results and the GE interpretation of the test data. LES calculates lower total pressure and higher total temperature on the pressure side of the stator, as the measured data showed, resulting in large loss generation on the pressure side of the stator. Detailed examination of the unsteady flow field from LES shows that the rotor wake, which has higher total temperature and higher total pressure relative to the free stream, interacts quite differently with the pressure side of the blade compared to the suction side of the blade. The higher temperature in the wake remains high as the wake passes through the pressure side of the blade. On the other hand, the total temperature diffuses as it passes through near the suction surface. For the presently investigated compressor, the classical intra-stator wake transport to the pressure side of the blade by the slip velocity in the wake seems to be minor. The main causes of this phenomenon are three-dimensional unsteady vortex interactions near the blade surface. The stabilizing effect of the concave curvature on the suction side keeps the rotor wake thin. On the other hand, the destabilizing effect of the convex curvature of the pressure side makes the rotor wake thicker, which results in a higher total temperature measurement at the stator exit. Additionally, wake stretching through the stator seems to contribute to the redistribution of the total temperature and the loss generation.

Efficiency↗

Unsteady Loss in the Stator Due to the Incoming Rotor Wake in a Highly-Loaded Transonic Compressor

The present paper reports an investigation of unsteady loss generation in the stator due to the incoming rotor wake in an advanced GE transonic compressor design with a high-fidelity numerical method. This advanced compressor with high reaction and high stage loading has been investigated both experimentally and analytically in the past. The measured efficiency in this advanced compressor is significantly lower than the design intention goal. The general understanding is that the current generation of compressor design analysis tools miss some important flow physics in this modern compressor design. To pinpoint the source of the efficiency miss, an advanced test with a detailed flow traverse was performed for the front one and a half stage at the NASA Glenn Research Center.

Compressor↗

Wavenumber-frequency Spectra of Pressure Fluctuations on a Generic Space Vehicle Measured via Unsteady Pressure-Sensitive Paint

Time histories of pressure fluctuations on a generic, hammerhead space vehicle model were measured using unsteady Pressure-Sensitive Paint (uPSP). The test was conducted in the 11-foot transonic wind tunnel of NASA Ames Research Center over a Mach number range of 0.6 M 1.2, and angles of attack of -4 4. The model was coated with a porous binder and PtTFPP-based porous polymer paint. An elaborate system of four high-speed cameras, and forty LED lamps was used for image acquisition. Various steps for image registration, reduction of shot noise, photogrammetry procedure to map images from the four cameras on a grid for the model, and finally a calibration procedure to convert the measured fluctuations in light intensity to fluctuating pressure, are discussed in the paper. The calibration process using a set of unsteady pressure sensors mounted on the model, was found to overcome some of the inherent problems of the fast response paint, such as rapid photo-degradation, non-linearity in pressure response, and significant temperature sensitivity. Comparison of spectra of pressure fluctuations between UPSP and pressure sensors demonstrated the ability of the paint to faithfully follow fluctuations up to 10 kHz, the maximum attempted. It was also found that the camera bit-depth and the illumination level limited the lowest measurable levels of pressure fluctuations to around 140dB. The large data set exposed various critical transonic flow physics not seen before, such as a coupling of the shock motion on the Payload Fairing (PF) with the separated flow region on the upper stage of the launch vehicle, and upstream convection of pressure fluctuation on PF at certain Mach numbers. The data also confirmed the expectation of a general lowering of the coefficient of pressure fluctuation with Mach number. The availability of the data set on a dense, regularly-spaced, surface grid allowed for the calculation of wavenumber-frequency (k-) spectra via straightforward applications of Fourier transform. The k- spectra were compared for the separated flow regions on the Second Stage, and the shock-boundary layer interactions on PF. The former showed self-similarity with Mach number while the latter was distinctly different, and confirmed the upstream propagation of pressure fluctuations. The k- spectra were dominated by the convected fluctuations; the acoustic domain was not discernable. These data, valuable for the vibro-acoustics analysis of aerospace vehicles, are believed to be the first obtained for the transonic flight regime, and pave the path for application on production models of aerospace vehicles.

buffet↗

History and Status of the CDISC Aerodynamic Design Method

This paper will review the development and application of the CDISC aerodynamic design method that began in the early 1980s and continues in active use today. The method uses an iterative, knowledge-based approach that has been coupled with numerous flow solvers utilizing a variety of grid types and levels of flow physics. It has been applied to configurations varying in complexity from 2-D airfoils to full 3-D aircraft at flow conditions ranging from low-speed, high-lift to supersonic cruise. The knowledge-based approach provides a very rapid design process, with designs typically completed in approximately 1-2x the time required for an initial converged analysis of the baseline. The evolution of the method as it tracked the expansive growth of grid and flow solver technology, along with some major applications over 5 decades, are given, followed by a description of the philosophy and process used in design. Finally, an example of its use in a preliminary design environment is given.

CDISC↗

History and Status of the CDISC Aerodynamic Design Method

This paper will review the development and application of the CDISC aerodynamic design method that began in the early 1980s and continues in active use today. The method uses an iterative, knowledge-based approach that has been coupled with numerous flow solvers utilizing a variety of grid types and levels of flow physics. It has been applied to configurations varying in complexity from 2-D airfoils to full 3-D aircraft at flow conditions ranging from low-speed, high-lift to supersonic cruise. The knowledge-based approach provides a very rapid design process, with designs typically completed in approximately 1-2x the time required for an initial converged analysis of the baseline. The evolution of the method as it tracked the expansive growth of grid and flow solver technology, along with some major applications over 5 decades, are given, followed by a description of the philosophy and process used in design. Finally, an example of its use in a preliminary design environment is given.

CDISC↗

PeleMP: The Multiphysics Solver for the Combustion Pele Adaptive Mesh Refinement Code Suite

Combustion encompasses multiscale, multiphase reacting flow physics spanning a wide range of scales from the molecular scales, where chemical reactions occur, to the device scales, where the turbulent flow is affected by the geometry of the combustor. This scale disparity and the limited measurement capabilities from experiments make modeling combustion a significant challenge. Recent advancements in high-performance computing (HPC), particularly with the Department of Energy's Exascale Computing Project (ECP), have enabled high-fidelity simulations of practical applications to be performed. The major physics submodels, including chemical reactions, turbulence, sprays, soot, and thermal radiation, exhibit distinctive computational characteristics that need to be examined separately to ensure efficient utilization of computational resources. This paper presents the multiphysics solver for the Pele code suite, called PeleMP, which consists of models for spray, soot, and thermal radiation. Here, the mathematical and algorithmic aspects of the model implementations are described in detail as well as the verification process. The computational performance of these models is benchmarked on multiple supercomputers, including Frontier, an exascale machine. Results are presented from production simulations of a turbulent sooting ethylene flame and a bluff-body swirl stabilized spray flame with sustainable aviation fuels to demonstrate the capability of the Pele codes for modeling practical combustion problems with multiphysics. This work is an important step toward the exascale computing era for high-fidelity combustion simulations providing physical insights and data for predictive modeling of real-world devices.

42 ENGINEERING↗

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental Investigation of the Near Wall Flow Structure of a Low Reynolds Number 3-D Turbulent Boundary Layer

Laser Doppler velocimetry (LDV) measurements and hydrogen bubble flow visualization techniques were used to examine the near-wall flow structure of 2D and 3D turbulent boundary layers (TBLs) over a range of low Reynolds numbers. The goals of this research were (1) an increased understanding of the flow physics in the near wall region of turbulent boundary layers,(2) to observe and quantify differences between 2D and 3D TBL flow structures, and (3) to document Reynolds number effects for 3D TBLs. The LDV data have provided results detailing the turbulence structure of the 2D and 3D TBLs. These results include mean Reynolds stress distributions, flow skewing results, and U and V spectra. Effects of Reynolds number for the 3D flow were also examined. Comparison to results with the same 3D flow geometry but at a significantly higher Reynolds number provided unique insight into the structure of 3D TBLs. While the 3D mean and fluctuating velocities were found to be highly dependent on Reynolds number, a previously defined shear stress parameter was discovered to be invariant with Reynolds number. The hydrogen bubble technique was used as a flow visualization tool to examine the near-wall flow structure of 2D and 3D TBLs. Both the quantitative and qualitative results displayed larger turbulent fluctuations with more highly concentrated vorticity regions for the 2D flow.

FLOW VISUALIZATION↗

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗

Aerodynamic shape optimization using control theory

Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Reuther, James↗