Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANS CFD”

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 127 records · Page 7

Investigating the Effects of Chemical Mechanism on Soot Formation Under High-Pressure Fuel Pyrolysis

We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.

02 PETROLEUM↗

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

This paper will describe three new supersonic jet cases for the NASA Turbulence Modeling Resource (TMR). These cases were taken from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test problem, which utilized data taken at NASA GRC consisting of non-intrusive measurements of jet plume velocities and temperatures, including both mean values and turbulence statistics. The three jet cases all used a Mach 1.63 axisymmetric nozzle: A heated jet for each at on-design conditions, a temperature-matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle. Computational fluid dynamics (CFD) solutions are obtained using three established Reynolds-averaged Navier-Stokes (RANS) codes: Wind-US, FUN3D, and Vulcan, all using the Menter Shear Stress Transport (SST) k-ω turbulence model. The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, that is - providing solutions that may not necessarily provide close comparisons with experimental data, but solve the posed CFD equations as intended.

nozzle↗

Validation Experiments of Incipient Turbulent Separation over an Axisymmetric Afterbody

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges in validation between experimental and numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires the support of advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed information of the flow field over a smooth body, undergoing adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states from fully attached,incipient separation, to fully separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configuration of the axisymmetric model features a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation. Experiments include steady pressure measurements to serve as preliminary comparisons to simulation studies, which examine the effect of variable grid domains and RANS turbulence models, in an effort to understand and evaluate the critical variability between solutions for the present configuration. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and inconsistent reference parameters. Ongoing work will focus on the experimental campaigns to obtain surface flow visualizations and high-resolution, off-body flow field measurements using Stereoscopic Particle Image Velocimetry (SPIV) and Laser Doppler Velocimetry (LDV) techniques to provide a detailed benchmark dataset to aid turbulence modelers.

axisymmetric afterbody↗

Development of the Glenn Heat-Transfer (Glenn-HT) Computer Code to Enable Time-Filtered Navier-Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations that are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminar/turbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes that take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-Heat-Transfer (Glenn-HT) code and applied to film-cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30deg holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and two blowing ratios of 0.5 and 1.0 are shown. Flow features under those conditions are also described.

Ameri, Ali↗

Development of the Glenn-HT Computer Code to Enable Time-Filtered Navier-Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations which are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminarturbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes which take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-HT code and applied to film cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30 holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and four blowing ratios of 0.5, 1.0, 1.5 and 2.0 are shown. Flow features under those conditions are also described.

infrared↗

Development of the Glenn-Heat-Transfer (Glenn-HT) Computer Code to Enable Time-Filtered Navier Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations that are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminar/turbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes that take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-Heat-Transfer (Glenn-HT) code and applied to film-cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30deg holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and two blowing ratios of 0.5 and 1.0 are shown. Flow features under those conditions are also described.

heat transfer↗

Dynamic Aeroelastic Flight Dynamic Modeling of Mach 0.745 Transonic Truss-Braced Wing

This paper presents a modeling approach for dynamic aeroelastic flight dynamic analysis of the Mach0.745 Transonic Truss-Braced Wing. The modeling approach is based on a transonic correction method to correct the Theodorsen’s theory for transonic flow. CFD unsteady Reynolds-averaged Navier-Stokes equations (RANS) simulations are conducted using FUN3D solver for a series of wing sections from the Mach 0.745Transonic Truss-Braced Wing in pitch and plunge oscillations. Unsteady lift and pitching moment coefficients are obtained and used to develop the correction terms in the Theodorsen’s theory to account for transonic aerodynamics. The unsteady lift and pitching moment derivatives with respect to the unsteady angle of attack are obtained as functions of the reduced frequency. These derivatives are used to compute the unsteady lift and pitching moment contributions by the angle of attack, pitch rate, and roll rate. They are then approximated using a frequency-domain regression to obtain the dynamic stability derivatives for the Mach 0.745 Transonic Truss-Braced Wing. The structural dynamic mode shapes of the Mach 0.745 TTBW are extracted from a NAS-TRAN finite-element model. These mode shapes are used to compute the generalized unsteady aerodynamic forces. The aerodynamic mass, damping, and stiffness and the aerodynamic lag states are constructed to couple the dynamic aeroelastic contribution to the flight dynamic model of the Mach 0.745 Transonic Truss-Braced Wing. The coupled dynamic aeroelastic flight dynamic equations of motion are formulated. The eigenvalues of the coupled system are computed. All the flight dynamic modes and structural dynamic modes are stable at Match 0.745. The effect of transonic aerodynamics generally causes all the dynamic modes to have lower damping values

Flight Dynamics↗

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

Several wall-modeled large eddy simulation (WMLES) methods are tested by simulating an unsteady Mach 2.0 compression corner geometry in Mississippi State’s Loci/CHEM solver. This study is conducted to evaluate the usage and requirements of these WMLES methods for applications regarding fluctuating pressure environments on launch vehicles with computational fluid dynamics (CFD). Two hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulation (LES) methods, Dynamic Hybrid RANS-LES (DHRL) and Improved Delayed Detached Eddy Simulation (IDDES), and one wall-stress-model, the Algebraic Wall Model for Wall-Modeled LES (AWMLES), are tested on varying grid and timestep refinement levels. These grid and timestep sizes are chosen to test the minimum requirements for successfully running these WMLES methods. The simulations are evaluated based on turbulent boundary layer properties in the developed boundary layer as well as unsteady quantities relating to fluctuating pressure environments in the region of the compression corner. The DHRL method shows good agreement with the comparison wind tunnel data and shows good grid and timestep convergence. The results from the IDDES and WMLES simulations show good agreement for several quantities with some discrepancies regarding others. The results presented in this paper will be used to inform further studies in predicting unsteady environments on higher-complexity geometries.

Bryson Frank↗

Advancing Thermal-Hydraulic Modeling Capabilities for Molten Salt Reactors in Pronghorn

This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.

42 ENGINEERING↗

A CFD Database for Airfoils and Wings at Post-Stall Angles of Attack

This paper presents selected results from an ongoing effort to develop an aerodynamic database from Reynolds-Averaged Navier-Stokes (RANS) computational analysis of airfoils and wings at stall and post-stall angles of attack. The data obtained from this effort will be used for validation and refinement of a low-order post-stall prediction method developed at NCSU, and to fill existing gaps in high angle of attack data in the literature. Such data could have potential applications in post-stall flight dynamics, helicopter aerodynamics and wind turbine aerodynamics. An overview of the NASA TetrUSS CFD package used for the RANS computational approach is presented. Detailed results for three airfoils are presented to compare their stall and post-stall behavior. The results for finite wings at stall and post-stall conditions focus on the effects of taper-ratio and sweep angle, with particular attention to whether the sectional flows can be approximated using two-dimensional flow over a stalled airfoil. While this approximation seems reasonable for unswept wings even at post-stall conditions, significant spanwise flow on stalled swept wings preclude the use of two-dimensional data to model sectional flows on swept wings. Thus, further effort is needed in low-order aerodynamic modeling of swept wings at stalled conditions.

Petrilli, Justin↗

Computational Study of Additively Manufactured Internally Cooled Airfoils for Industrial Gas Turbine Applications

Internal cooling features such as pin-fins, impingement jets, and rib-turbulators are necessary to keep turbine components cool, but if sufficiently advanced can potentially also eliminate the need for film cooling on turbine blades particularly in industrial gas turbines where temperatures are not extreme. Furthermore, by leveraging additive manufacturing, other advanced designs such as lattice and incremental impingement configurations are possible and have recently been experimentally tested. While the performance of such configurations has been quantified through means of overall cooling effectiveness, it is not as clear why certain designs were better than others, or what the mechanisms were behind the observed external cooling patterns. The purpose of this study was to computationally analyze different advanced turbine blade internal cooling designs previously tested by the National Energy Technology Laboratory.

CFD↗

Supersonic Combusting Jet Experiments for Code Development and Validation

Computational fluid dynamics (CFD) methods based on the Reynolds averaged Navier-Stokes (RANS) equations are extensively employed in the design of hypersonic airbreathing engines. A fundamental weakness in these methods is the accurate mathematical modeling of turbulence and turbulence-combustion interactions. This paper is one of a companion pair of papers which describe the use of the combined dual pump CARS and newly developed interferometric Rayleigh scattering optical system to acquire time-accurate turbulence and mixing data with which to aid calibration and verification of the current turbulence models. This test program consisted of a pilot laboratory-scale experiment used to develop the optical system and data acquisition techniques, and a significantly larger experiment to acquire the necessary turbulence data. The flow is an axially-symmetric, supersonic, combusting, free jet that provides good optical access, consisting of a central jet of hot "vitiated air" and a coflow jet of hydrogen or ethylene fuel. In this paper the development of both experiments is described. Facility and flow visualization data are presented for various types of flames, including the flames selected for detailed study with the CARS/Rayleigh optical techniques.

Cutler, Andrew D.↗

Computational Study of Flow Interactions in Coaxial Rotors

Although the first idea of coaxial rotors appeared more than 150 years ago, most helicopters have used single main-rotor/tail-rotor combination. Since reactive moments of coaxial rotors are canceled by contra-rotation, no tail rotor is required to counter the torque generated by the main rotor. Unlike the single main rotor design that distributes power to both main and tail rotors, all of the power for coaxial rotors is used for vertical thrust. Thus, no power is wasted for anti-torque or directional control. The saved power helps coaxial rotors reach a higher hover ceiling than single rotor helicopters. Another advantage of coaxial rotors is that the overall rotor diameter can be reduced for a given vehicle gross weight because each rotor provides a maximum contribution to vertical thrust to overcome vehicle weight. However, increased mechanical complexity of the hub has been one of the challenges for manufacturing coaxial rotorcraft. Only the Kamov Design Bureau of Russia had been notably successful in production of coaxial helicopters until Sikorsky built X2, an experimental compound helicopter. Recent developments in unmanned aircraft systems and high-speed rotorcraft have renewed interest in the coaxial configuration. Multi-rotors are frequently used for small electric unmanned rotorcraft partly due to mechanical simplicity. The use of multiple motors provides redundancy as well as cost-efficiency. The multi-rotor concept has rarely been used until recently because of its inherent stability and control problems. However, advances in inexpensive electronic flight control systems have opened the floodgates for small drones using multirotors. Coaxial rotors have started to appear in some multi-rotor configurations. Small coaxial rotors have often been designed using a hundred year old approach that is "sketch, build, fly, and iterate." In that approach, there is no systematic way to explore trade-offs or determine logical next steps. It is neither possible to account for multiple real-world constraints up front in design nor possible to know what performance is possible with a given design. Since unmanned vehicles are sized and optimized for the particular mission, a modern low-fidelity conceptual design and sizing tool that has been used for the design of large helicopters can be used for design of small coaxial rotorcraft. However, unlike most helicopters with single main rotor, the interactions between the upper and lower rotors emerge as an important factor to consider in design because an increase in performance of a multi-rotor system is not proportional to the number of rotors. Interference losses and differences in thrusts between the upper and lower rotors were investigated by theoretical methods as well as a computational fluid dynamics (CFD) method using the Reynolds-Averaged Navier-Stokes (RANS) equations. In this work, hybrid turbulence models are used to investigate the physics of interactions between coaxial rotors and a fuselage that are not well understood. Present study covers not only small-scale drones but also large-scale coaxial rotors for heavy-lifting missions. Considering the recently proposed FAA drone rules that require the flight only in visual line-of-sight, a large multirotor might be used as an airborne carrier for launch and recovery of unmanned aircraft systems with a human operator onboard. For applications to civil operations, their aerodynamic performance and noise levels need to be assessed. Noise is one of the largest limiting factors to rotorcraft operations in urban area. Since the high-frequency noise of multi-rotors may increase the annoyance, noise may turn out to be a key issue that must be addressed for market acceptability. One of the objectives of the present work is to study the effects of inter-rotor spacing and collectives on the performance, efficiency, and acoustics of coaxial rotor systems.

Study↗

Comparative Results from a CFD Challenge Over a 2D Three-Element High-Lift Airfoil

A high-lift workshop was held in May of 1993 at NASA Langley Research Center. A major part of the workshop centered on a blind test of various computational fluid dynamics (CFD) methods in which the flow about a two- dimensional (2D) three-element airfoil was computed without prior knowledge of the experimental data. The results of this 'blind' test revealed: (1) The Reynolds Averaged Navier-Stokes (RANS) methods generally showed less variability among codes than did potential/Euler solvers coupled with boundary-layer solution techniques. However, some of the coupled methods still provided excellent predictions. (2) Drag prediction using coupled methods agreed more closely with experiment than the RANS methods. Lift was more accurately predicted than drag for both methods. (3) The CFD methods did well in predicting lift and drag changes due to changes in Reynolds number, however, they did not perform as well when predicting lift and drag increments due to changing flap gap, (4) Pressures and skin friction compared favorably with experiment for most of the codes. (5) There was a large variability in most of the velocity profile predictions. Computational results predict a stronger siat wake than measured suggesting a missing component in turbulence modeling, perhaps curvature effects.

Klausmeyer, Steven M.↗

Flow Control Analysis on the Hump Model with RANS Tools

A concerted effort is underway at NASA Langley Research Center to create a benchmark for Computational Fluid Dynamic (CFD) codes. both unstructured and structured, against a data set for the hump model with actuation. The hump model was tested in the NASA Langley 0.3-m Transonic Cryogenic Tunnel. The CFD codes used for the analyses are the FUN2D (Full Unstructured Navier-Stokes 2-Dimensional) code, the structured TLNS3D (Thin-Layer Navier-Stokes 3-Dimensional) code, and the structured CFL3D code, all developed at NASA Langley. The current investigation uses the time-accurate Reynolds-Averaged Navier-Stokes (RANS) approach to predict aerodynamic performance of the active flow control experimental database for the hump model. Two-dimensional computational results verified that steady blowing and suction and oscillatory suction/blowing can be used to significantly reduce the separated flow region on the model. Discrepancies do exist between the CFD results and experimental data in the region downstream of the slot with the largest differences in the oscillatory cases. Overall, the structured CFD codes exhibited similar behavior with each other for a wide range of control conditions, with the unstructured FUN2D code showing moderately different results in the separated flow region for the suction and oscillatory cases.

Viken, Sally A.↗

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗

Simulations of Mach 0.8 Transonic Truss-Braced Wing Aircraft Aerodynamics at High Angles of Attack

This paper presents numerical simulations of the Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft at high angle of attack using CFD solver FUN3D. Three different approaches - steady Reynolds Averaged Navier-Stokes (RANS), Unsteady Reynolds Averaged Navier-Stokes (URANS), and Delayed Detached Eddy Simulation (DDES) are used to simulation the TTBW aircraft at Mach number 0.8 and altitude 40, 000 ft. The pitch break phenomenon is observed at this flight condition. The simulation results show that the pitch break angle of attack occurs at 2.75˚ by DDES, 2.81◦ by URANS, and 2.89˚ by RANS. The DDES unsteady pressure begins to grow at angle of attack of 3.5˚, which might indicate the buffet onset, and rises one order of magnitude larger at angle of attack of 4.0˚. However, the URANS unsteady pressure begins to grow at angle of attack of 4.0˚. Both unsteady simulation results show the pitch break occurs before buffet onset.

TTBW↗

Recent Improvements in Pronghorn for Advanced Reactor Modeling

Pronghorn is a thermal-hydraulics computational tool developed using the Idaho National Laboratory's Multiphysics Object-Oriented Simulation Environment (MOOSE). It is designed to support Computational Fluid Dynamics (CFD) modeling, ranging from subchannel and porous media analysis to Reynolds Averaged Navier-Stokes (RANS) turbulence modeling. As an integral part of the MOOSE-based suite of tools, Pronghorn seamlessly couples with other MOOSE-based applications to simulate a variety of physical phenomena. This article highlights recent significant enhancements to Pronghorn's CFD modeling capabilities and demonstrates their application to advanced nuclear reactor designs. The recent improvements in Pronghorn primarily focus on modifications to its turbulence modeling capabilities, near-wall corrections and numerical schemes. In terms of turbulence modeling, the two-equation $k-\epsilon$ and $k-\omega$ SST models have been implemented and validated with both equilibrium and non-equilibrium wall treatments. Additionally, corrections for wall roughness, and curvature, and wall-channeling in pebble beds have been introduced in the near-wall modeling. These developments enable more accurate simulations of advanced nuclear reactors. Two case studies are presented in this work: a pool-type Molten Chloride Reactor and a salt-cooled Pebble-Bed High Temperature Reactor. In both cases, the previous models in Pronghorn are compared with the new implementations, demonstrating the improved accuracy achieved with the updated models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗