Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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 1,027 records · Page 57

Economical Unsteady High-Fidelity Aerodynamics for Structural Optimization with a Flutter Constraint

Structural optimization with a flutter constraint for a vehicle designed to fly in the transonic regime is a particularly difficult task. In this speed range, the flutter boundary is very sensitive to aerodynamic nonlinearities, typically requiring high-fidelity Navier-Stokes simulations. However, the repeated application of unsteady computational fluid dynamics to guide an aeroelastic optimization process is very computationally expensive. This expense has motivated the development of methods that incorporate aspects of the aerodynamic nonlinearity, classical tools of flutter analysis, and more recent methods of optimization. While it is possible to use doublet lattice method aerodynamics, this paper focuses on the use of an unsteady high-fidelity aerodynamic reduced order model combined with successive transformations that allows for an economical way of utilizing high-fidelity aerodynamics in the optimization process. This approach is applied to the common research model wing structural design. As might be expected, the high-fidelity aerodynamics produces a heavier wing than that optimized with doublet lattice aerodynamics. It is found that the optimized lower skin of the wing using high-fidelity aerodynamics differs significantly from that using doublet lattice aerodynamics.

Bartels, Robert E.↗

Isolated Rotor Forward Flight Testing from One Atmosphere down to Martian Atmospheric Densities

With the recent interest in Martian exploration using Unmanned Aerial Vehicles (UAV), an experimental study was conducted to investigate rotor performance at Martian atmospheric conditions. Both simulation and testing of rotors is vital for the evaluation of performance and behavior of a rotor, especially when subjected to a Marian atmosphere. One critical test that has not been performed to date is simulated helicopter forward flight in a Martian atmosphere. To achieve this, the test must be conducted in a facility which can be evacuated to the atmospheric pressure and density of Mars. A unique 40-in diameter rotor, roughly approximating a proposed design for a Mars Helicopter (MH), was tested in forward flight at Mars atmospheric pressure at the NASA Ames Planetary Aeolian Laboratory (PAL). The goal of this experiment was to collect rotor thrust, RPM, power, torque, and acoustics measurements. Subsequently, these results are compared with simulated cases using a mid-fidelity Computational Fluid Dynamics (CFD) approach. As expected, rotor thrust and power results are drastically reduced when under low atmospheric conditions. In addition, Reynolds number effects seem to play a vital role that cannot be neglected.

Ament, Geoffrey A.↗

Hybrid CFD Engineering Model of Plume Induced Erosion and Crater Formation During Descent of Lunar Landers

With rapidly increased worldwide interest in landing on the moon, the issue of Plume Surface Interactions (PSI) is gaining attention. Hazards posed by lander plume induced dust and debris, as well as landing site deformation can be mitigated when better understood through predictive simulations. As simulation enabling computational power continues to increase, hybrid Computational Fluid Dynamics (CFD)/Engineering models provide the immediate ability to conduct parametric/trade studies driving design decisions for landers. Reduced order erosion engineering models apply correlations of the surface erosion rate to the plume induced surface forces with the correlations anchored to flight observations from Apollo LM landings. The MSFC propulsion fluid dynamics branch has developed such a hybrid model for predicting the plume induced viscous erosive regression of the Lunar surface beneath a landing vehicle. The (Descent Interpolated Gas Granular Erosion Model) DIGGEM was originally implemented as a post processing tool to calculate induced erosion rates through vehicle descent using CFD solutions of the vehicle plume at several fixed altitudes. This capability has since been advanced in the Loci/Chem-DIGGEM model to allow transient, moving vehicle, fully coupled viscous erosion modeling of vehicle descent PSI. This model has recently been used to make preflight predictions of the erosive regression of the ground beneath a Commercial Lunar Payload Services (CLPS) vehicle in support of measurements to be made by the Stereo Cameras for Lunar Plume Surface Studies (SCALPSS) instrument.

Plume Surface Interaction↗

Hybrid CFD Engineering Model of Plume Induced Erosion and Crater Formation During Descent of Lunar Landers

With rapidly increased worldwide interest in landing on the moon, the issue of Plume Surface Interactions (PSI) is gaining attention. Hazards posed by lander plume induced dust and debris, as well as landing site deformation can be mitigated when better understood through predictive simulations. As simulation enabling computational power continues to increase, hybrid Computational Fluid Dynamics (CFD)/Engineering models provide the immediate ability to conduct parametric/trade studies driving design decisions for landers. Reduced order erosion engineering models apply correlations of the surface erosion rate to the plume induced surface forces with the correlations anchored to flight observations from Apollo LM landings. The MSFC propulsion fluid dynamics branch has developed such a hybrid model for predicting the plume induced viscous erosive regression of the Lunar surface beneath a landing vehicle. The (Descent Interpolated Gas Granular Erosion Model) DIGGEM was originally implemented as a post processing tool to calculate induced erosion rates through vehicle descent using CFD solutions of the vehicle plume at several fixed altitudes. This capability has since been advanced in the Loci/Chem-DIGGEM model to allow transient, moving vehicle, fully coupled viscous erosion modeling of vehicle descent PSI. This model has recently been used to make preflight predictions of the erosive regression of the ground beneath a Commercial Lunar Payload Services (CLPS) vehicle in support of measurements to be made by the Stereo Cameras for Lunar Plume Surface Studies (SCALPSS) instrument.

Plume Surface Interaction↗

String Potentiometers for Supersonic Split-Line Thrust Vector Amplification Factor Verification

The static fire of NASA’s Mars Ascent Vehicle (MAV) first stage Solid Rocket Motor (SRM1) Demonstration Motor-1 (DM-1) occurred at Edwards Air Force Base (AFB) on April 7, 2023. Among other objectives, this test successfully demonstrated the Thrust Vector Control (TVC) system, including the Technology Readiness Level (TRL) of the Supersonic Split-Line (SSSL) nozzle. An essential characteristic of the SSSL is the thrust vector amplification factor (TVAF) it exhibits; however, empirically quantifying this behavior is difficult due to the limitations of measuring the internal nozzle gas dynamics of a Solid Rocket Motor (SRM) during burn. Using four string potentiometers mounted inside the test chamber connected to the nozzle and motor case and a thrust measurement system (TMS), the TVAF can be derived from the known nozzle and motor geometry. Increasing the accuracy of the TVAF is essential for certain groups working on the MAV including the Guidance, Navigation, and Control (GN&C) group, who prior to DM-1, solely relied on Computational Fluid Dynamics (CFD) data in their Monte Carlo simulations for TVAF. The methodology shown in this paper can be generally applied to any test article implementing a SSSL nozzle, presuming string potentiometers are a viable option and the facility is conducive of their implementation.

Supersonic Split-Line↗

Aerodynamics of a Uranus Aerocapture System Using a Mars-Heritage Entry Vehicle

Aerodynamic characteristics of an aerocapture system intended to deliver a flagship-class orbiter and probe planetary science mission to Uranus are presented. The aeroshell of the Mars Science Laboratory and Mars 2020 entry vehicles is proposed as a baseline for this system to reduce the amount of necessary technology development. Direct Simulation Monte Carlo and Navier-Stokes computational fluid dynamics solutions are used to characterize the aerodynamic performance of the Mars-heritage vehicle for aerocapture flight at Uranus. These results are incorporated into an aerodatabase for use in six degree-of-freedom trajectory studies and mission design. Updates are made to the Mars-heritage aerodynamic uncertainty model based on observations in the Uranus-specific computational data to ensure the model is conservatively bounding for the proposed flight space. Necessary modifications to the aeroshell for system packaging are found to have minimal effect on aerodynamic performance. The resulting aerodatabase and uncertainty model are used to show the existing Mars-heritage entry vehicles have sufficient aerodynamic performance to achieve required control margin for Uranus aerocapture.

Eli R Shellabarger↗

Aerodynamics of a Uranus Aerocapture System Using a Mars-Heritage Entry Vehicle

Aerodynamic characteristics of an aerocapture system intended to deliver a flagship-class orbiter and probe planetary science mission to Uranus are presented. The aeroshell of the Mars Science Laboratory and Mars 2020 entry vehicles is proposed as a baseline for this system to reduce the amount of necessary technology development. Direct Simulation Monte Carlo and Navier-Stokes computational fluid dynamics solutions are used to characterize the aerodynamic performance of the Mars-heritage vehicle for aerocapture flight at Uranus. These results are incorporated into an aerodatabase for use in six degree-of-freedom trajectory studies and mission design. Updates are made to the Mars-heritage aerodynamic uncertainty model based on observations in the Uranus-specific computational data to ensure the model is conservatively bounding for the proposed flight space. Necessary modifications to the aeroshell for system packaging are found to have minimal effect on aerodynamic performance. The resulting aerodatabase and uncertainty model are used to show the existing Mars-heritage entry vehicles have sufficient aerodynamic performance to achieve required control margin for Uranus aerocapture.

Eli Shellabarger↗

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗

High-Fidelity Multidisciplinary Design Optimization Methodology with Application to Rotor Blades

A multidisciplinary design optimization procedure has been developed and applied to rotorcraft simulations involving tightly-coupled, high-fidelity computational fluid dynamics and comprehensive analysis. A discretely-consistent, adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows on unstructured, dynamic, overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system for high-fidelity rotorcraft analysis is verified; simulation results exhibit good agreement with established solutions. A constrained gradient-based design optimization for a HART-II rotorcraft configuration is demonstrated. The computational cost for individual components of the multidisciplinary sensitivity analysis is assessed and improved.

Wang, Li↗

Tensoral: A system for post-processing turbulence simulation data

Many computer simulations in engineering and science -- and especially in computational fluid dynamics (CFD) -- produce huge quantities of numerical data. These data are often so large as to make even relatively simple post-processing of this data unwieldy. The data, once computed and quality-assured, is most likely analyzed by only a few people. As a result, much useful numerical data is under-utilized. Since future state-of-the-art simulations will produce even larger datasets, will use more complex flow geometries, and will be performed on more complex supercomputers, data management issues will become increasingly cumbersome. My goal is to provide software which will automate the present and future task of managing and post-processing large turbulence datasets. My research has focused on the development of these software tools -- specifically, through the development of a very high-level language called 'Tensoral'. The ultimate goal of Tensoral is to convert high-level mathematical expressions (tensor algebra, calculus, and statistics) into efficient low-level programs which numerically calculate these expressions given simulation datasets. This approach to the database and post-processing problem has several advantages. Using Tensoral the numerical and data management details of a simulation are shielded from the concerns of the end user. This shielding is carried out without sacrificing post-processor efficiency and robustness. Another advantage of Tensoral is that its very high-level nature lends itself to portability across a wide variety of computing (and supercomputing) platforms. This is especially important considering the rapidity of changes in supercomputing hardware.

Dresselhaus, Eliot↗

Computational Characterization and Boundary Condition Models of a Jet Interaction Fluidic Oscillator

Fluidically oscillating jet actuators that generate sweeping jets when supplied with a pressurized fluid have been used to mitigate separation and reduce drag in a range of flow control applications. Their implementation in future flight vehicles would require fundamental understanding and accurate predictive techniques of the physics of their internal flow and jet formation. The present investigations focus on high-fidelity, time-accurate simulations to characterize the flow physics of the actuator in quiescent conditions. An important element of the present simulations is to demonstrate the ability of the computational fluid dynamics (CFD) solver to predict the jet characteristics and provide a basis for the development of improved boundary conditions (BC) without entirely resolving the geometrical features of the fluidic device. The CFD-predicted oscillation frequencies of the engendered jets were found to be in excellent agreement with experiments, even on two-dimensional meshes. The study revealed that three-dimensional simulations are required to capture some of the flow features of the sweeping jet such as the double peak in time-averaged velocity distributions downstream of the actuator’s orifice that were measured in experiments. Furthermore, several approaches for modeling the actuator were implemented and assessed in quiescent conditions. The evaluation of a boundary condition at the device throat, based on the phase-averaged flow variables, provides the basis for devising surface-based boundary conditions. The influence and necessity of including turbulent characteristics as part of the boundary conditions have also been identified.

42 ENGINEERING↗

CFD code calibration and inlet-fairing effects on a 3D hypersonic powered-simulation model

A three-dimensional (3D) computational study has been performed addressing issues related to the wind tunnel testing of a hypersonic powered-simulation model. The study consisted of three objectives. The first objective was to calibrate a state-of-the-art computational fluid dynamics (CFD) code in its ability to predict hypersonic powered-simulation flows by comparing CFD solutions with experimental surface pressure dam. Aftbody lower surface pressures were well predicted, but lower surface wing pressures were less accurately predicted. The second objective was to determine the 3D effects on the aftbody created by fairing over the inlet; this was accomplished by comparing the CFD solutions of two closed-inlet powered configurations with a flowing-inlet powered configuration. Although results at four freestream Mach numbers indicate that the exhaust plume tends to isolate the aftbody surface from most forebody flowfield differences, a smooth inlet fairing provides the least aftbody force and moment variation compared to a flowing inlet. The final objective was to predict and understand the 3D characteristics of exhaust plume development at selected points on a representative flight path. Results showed a dramatic effect of plume expansion onto the wings as the freestream Mach number and corresponding nozzle pressure ratio are increased.

Huebner, Lawrence D.↗

CFD Code Calibration and Inlet-Fairing Effects On a 3D Hypersonic Powered-Simulation Model

A three-dimensional (3D) computational study has been performed addressing issues related to the wind tunnel testing of a hypersonic powered-simulation model. The study consisted of three objectives. The first objective was to calibrate a state-of-the-art computational fluid dynamics (CFD) code in its ability to predict hypersonic powered-simulation flows by comparing CFD solutions with experimental surface pressure data. Aftbody lower surface pressures were well predicted, but lower surface wing pressures were less accurately predicted. The second objective was to determine the 3D effects on the aftbody created by fairing over the inlet; this was accomplished by comparing the CFD solutions of two closed-inlet powered configurations with a flowing- inlet powered configuration. Although results at four freestream Mach numbers indicate that the exhaust plume tends to isolate the aftbody surface from most forebody flow- field differences, a smooth inlet fairing provides the least aftbody force and moment variation compared to a flowing inlet. The final objective was to predict and understand the 3D characteristics of exhaust plume development at selected points on a representative flight path. Results showed a dramatic effect of plume expansion onto the wings as the freestream Mach number and corresponding nozzle pressure ratio are increased.

Huebner, Lawrence D.↗

Visual Computing Environment

The Visual Computing Environment (VCE) is a NASA Lewis Research Center project to develop a framework for intercomponent and multidisciplinary computational simulations. Many current engineering analysis codes simulate various aspects of aircraft engine operation. For example, existing computational fluid dynamics (CFD) codes can model the airflow through individual engine components such as the inlet, compressor, combustor, turbine, or nozzle. Currently, these codes are run in isolation, making intercomponent and complete system simulations very difficult to perform. In addition, management and utilization of these engineering codes for coupled component simulations is a complex, laborious task, requiring substantial experience and effort. To facilitate multicomponent aircraft engine analysis, the CFD Research Corporation (CFDRC) is developing the VCE system. This system, which is part of NASA's Numerical Propulsion Simulation System (NPSS) program, can couple various engineering disciplines, such as CFD, structural analysis, and thermal analysis. The objectives of VCE are to (1) develop a visual computing environment for controlling the execution of individual simulation codes that are running in parallel and are distributed on heterogeneous host machines in a networked environment, (2) develop numerical coupling algorithms for interchanging boundary conditions between codes with arbitrary grid matching and different levels of dimensionality, (3) provide a graphical interface for simulation setup and control, and (4) provide tools for online visualization and plotting. VCE was designed to provide a distributed, object-oriented environment. Mechanisms are provided for creating and manipulating objects, such as grids, boundary conditions, and solution data. This environment includes parallel virtual machine (PVM) for distributed processing. Users can interactively select and couple any set of codes that have been modified to run in a parallel distributed fashion on a cluster of heterogeneous workstations. A scripting facility allows users to dictate the sequence of events that make up the particular simulation.

Lawrence, Charles↗

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗