Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Progress Towards a Microgravity CFD Validation Study Using the ISS SPHERES-SLOSH Experiment

Understanding, predicting, and controlling fluid slosh dynamics is critical to safety and improving performance of space missions when a significant percentage of the spacecrafts mass is a liquid. Computational fluid dynamics simulations can be used to predict the dynamics of slosh, but these programs require extensive validation. Many CFD programs have been validated by slosh experiments using various fluids in earth gravity, but prior to the ISS SPHERES-Slosh experiment1, little experimental data for long-duration, zero-gravity slosh existed. This paper presents the current status of an ongoing CFD validation study using the ISS SPHERES-Slosh experimental data.

Storey, Jedediah M.↗

Lattice Boltzmann and Navier-Stokes Cartesian CFD Approaches for Airframe Noise Predictions

Lattice Boltzmann (LB) and compressible Navier-Stokes (NS) equations based computational fluid dynamics (CFD) approaches are compared for simulating airframe noise. Both LB and NS CFD approaches are implemented within the Launch Ascent and Vehicle Aerodynamics (LAVA) framework. Both schemes utilize the same underlying Cartesian structured mesh paradigm with provision for local adaptive grid refinement and sub-cycling in time. We choose a prototypical massively separated, wake-dominated flow ideally suited for Cartesian-grid based approaches in this study - The partially-dressed, cavity-closed nose landing gear (PDCC-NLG) noise problem from AIAA's Benchmark problems for Airframe Noise Computations (BANC) series of workshops. The relative accuracy and computational efficiency of the two approaches are systematically compared. Detailed comments are made on the potential held by LB to significantly reduce time-to-solution for a desired level of accuracy within the context of modeling airframes noise from first principles.

Airframe Noise Prediction↗

CFD Predictions of N+3 Cycle Emissions for a Three-Cup Gas-Turbine Combustor

The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for a unique pre-filming type LDI-3 fuel injector for a three-cup, nineteen-element flametube configuration. All computations were performed with a consistent approach of mesh-generation, spray modeling, reduced finite-rate kinetics and turbulence-chemistry interaction, as developed for CFD analysis of single-element and multi-element LDI-3 designs with OpenNCC. Emissions and flowfield characteristics were predicted for a generic NASA N+3 engine cycle, with particular focus on the 7% and 30% ICAO power operating conditions. For both the conditions studied, the CFD analysis provided very good predictions for EINOx when compared with experimental data measured at NASA Glenn Research Center.

injectors↗

CFD Model Of The Transonic Dynamics Tunnel With Applications

This paper presents the Computational Fluid Dynamics (CFD) model of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) with some recent applications. The TDT is a continuous-flow, closed circuit, slotted-test-section wind tunnel with a 16- by 16-foot test section with cropped corners. The tunnel was originally built as the 19-ft Pressure Tunnel in 1938, but it was converted to the current transonic tunnel in the 1950s, with capabilities to use either air or heavy gas at pressures from atmosphere down to near vacuum. In this study, experimental data acquired in the empty tunnel using R-134a as the test medium was used to calibrate the computational data. Experimental data from a recent TDT test of a full-span fighter configuration in air was then selected for comparison with the numerical data. During this test, the configuration experienced a flutter event in the transonic flow regime. Numerically, the flutter event is simulated both inside the CFD model of the TDT and in a classical free-air model. The preliminary results show that the wind-tunnel walls do not affect flutter prediction.

Chwalowski, Pawel↗

Rotor CFD Analysis at Terrestrial and Martian Atmospheric Densities

Much effort has been made to enhance exploration on Mars. In addition to a rover and Mars-orbiting satellites, a Mars helicopter (MH) was proposed in order to augment planetary research. Computational Fluid Dynamics (CFD) simulations have been performed to have a better understanding of the behavior and performance of vertical lift Planetary Aerial Vehicles (PAV). Due to the large differences in atmospheric conditions between Mars and Earth, predicting and testing rotorcraft performance is a complex task. The goal of this project is to understand the capability of the mid-fidelity CFD software RotCFD to predict rotor performance in terms of thrust at 1013.25 milibar and 14 milibar corresponding to Terrestrial and Martian conditions, respectively. Also, in order to characterize the wind tunnel wall effects free field and wind tunnel simulations were performed, analyzed and compared. Different analytical tools have been used in order to aid with the design process for the future vertical lift planetary aerial vehicles. One of them includes experimental tests performed on a rotor in the Aeolian Wind Tunnel (AWT) facility at NASA Ames Research Center under different pressure conditions ranging from Terrestrial to Martian atmospheric conditions. Other software was used as well in order to capture the aerodynamic coefficients of the corresponding rotor sections based on the Mach and Reynolds numbers used for the experimental tests. The aerodynamic coefficients were input into RotCFD, and various simulations were performed under Terrestrial and Martian conditions in order to mimic the experimental test. Then, the obtained results from RotCFD were compared with the AWT collected data.

Atmospheric Pressures↗

CFD Analysis in Advance of the NASA Juncture Flow Experiment

Outline: Experiment Motivation, Goals, Model Design; Wing Candidates; Risk reduction experiments -NASA Ames Test Cell 2 (TC2) 32 inch by 48 inch, 3 percent semispan -Virginia Tech Stability Tunnel 6 foot, 2.5 percent fullspan -NASA Langley 14 by 22 Foot Subsonic Tunnel (14 by 22) 6 percent fullspan; Results from 14 by 22 6 percent risk reduction -CFD (Computational Fluid Dynamics) Free Air -CFD with 14 by 22 WT (Wind Tunnel) walls -Risk Reduction Experiment oil flow; Observations and Upcoming Experiment.

Lee, Henry↗

CFD Evaluation of Lean-Direct Injection Combustors for Commercial Supersonics Technology

An overview is given of an effort that focused on using CFD analysis to complement design and configuration definition of Lean-Direct Injection (LDI) combustion concepts for NASA's Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for second and third generation LDI configurations at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray modeling, ignition and kinetics modeling. Emissions (EINOx) characteristics were predicted for CST cruise conditions, and compared with emissions data from experimental measurements to evaluate the fidelity of the CFD modeling approach to predict emissions changes in response to changes in supersonic cycle conditions.

nitrogen oxides injectors↗

CFD Evaluation Of Lean-Direct Injection Combustors for Commercial Supersonics Technology

An overview is given of an effort that focused on using CFD analysis to complement design and configuration definition of Lean-Direct Injection (LDI) combustion concepts for NASA's Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for second and third generation LDI configurations at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray modeling, ignition and kinetics modeling. Emissions (EINOx) characteristics were predicted for CST cruise conditions, and compared with emissions data from experimental measurements to evaluate the fidelity of the CFD modeling approach to predict emissions changes in response to changes in supersonic cycle conditions.

Ajmani, Kumud↗

CFD and Experimental Data Comparisons for Conventional and AFC-Enabled CRM High-Lift Configurations

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of a representative conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface of the flap shoulder, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-Configuration to produce the aerodynamic performance (lift) comparable to that of the CRM-HL configuration over a large angle of attack range encompassing maximum lift conditions. For comparison purposes, computations are also per-formed for the CRM-HL configuration. In the current paper, PowerFLOWR©, a CFD code based on the Lattice Boltzmann method (LBM), is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with a combination of two different types of AFC actuators. The transonic version of the PowerFLOWR© code that has been validated for high-speed flows is used to simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. This study is focused on the AFC systems and actuator arrangements that had emerged based on the parametric studies conducted previously at the nominal landing condition. Comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations over a broad angle of attack range. The numerical solutions predict the expected trends in aerodynamic forces with angle of attack variation.

Veer N Vatsa↗

CFD Simulation of an Inward-Turning Supersonic Inlet Unstart at Flight Mach Number 1.7

Interest in inward-turning supersonic inlets renewed during NASA’s Commercial Supersonic Technology Project, as a means of reducing disturbances caused by propulsion integration. Inward-turning designs have the potential to greatly reduce the strength of shock waves emanating outward from the nacelle, and therefore their contribution to the overall sonic boom signature. For top-mounted, or over-wing installations, disturbances to aircraft aerodynamics are reduced. The nature of the inward-turning geometry results in a boundary-layer thickening effect, opposite that of the familiar thinning of an outward-turning cone boundary-layer. The application of a small amount of boundary-layer bleed at the throat has been shown to alleviate the negative effects of the boundary-layer on subsonic diffuser performance. Boundary-layer bleed, however, comes at the expense of increased complexity, and external disturbances at bleed exits that would partially defeat the purpose of the inward-turning scheme. Vortex generators were proposed as an alternative to bleed, and an initial CFD study of their effectiveness in this application forms the basis of the present work. A source-term model in the Wind-US CFD code was used to model the effects of an array of vortex generators placed at the entrance to the subsonic diffuser. Baseline results with no vortex generators compared favorably to previously reported results with no bleed at back-pressures near the design point, but stability margin in the present case was limited by a curious non-linear phenomena similar to unstart in mixed-compression inlets. While some improvement was noted with vortex generators, their effectiveness was diminished by the limited stability margin. A time-accurate calculation was subsequently performed to better understand this process. Results of this calculation are presented and instantaneous solutions compared with steady-state, subcritical results previously reported. The stability margin limitation is attributed to the use of a constant-pressure outflow boundary condition, which differed from that used in the prior simulations.

Pre-Compressor Cooling, Water Injection, Hypersoni↗

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can either be used to improve the model through calibration or quantify the model-form uncertainty. This work focuses on providing the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given uncertainty in both the simulation and experimental data. A procedure will be presented that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs. The work will also assess the impact between using different high-fidelity solvers and different turbulence models.The goal of performing this comparison is to provide a quantifiable measurement for the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis will also be performed to get a model-form uncertainty for the lower fidelity analysis. The two high-fidelty CFD solvers used are VULCAN-CFD and FUN3D and the low-fidelity results will come from Cart3D. The experimental data came from the 20-inch Mach 6Tunnel located at NASA Langley Research Center.

Laura M. White↗

Evaluation of CFD as a Surrogate for Wind Tunnel Testing - Experimental Uncertainty Quantification

A series of wind tunnel tests is being performed at the Unitary Plan Wind Tunnel at Langley Research Center to assess the validity of using computational fluid dynamics (CFD) as a surrogate for wind tunnel testing. In order to make proper comparisons, uncertainties in CFD results and experimental data must be well understood. The material presented will highlight the methods, assumptions, and elemental inputs used to achieve experimental uncertainty estimates for several variables of interest.The work performed to date has focused on gaining insight into random uncertainty via statistical analysis of repeat data and systematic uncertainty via Monte Carlo propagation analysis. These methods were combined using a second-order Monte Carlo propagation, resulting in a probability box (P-box) for each variable and at all conditions of interest. Figure 1 provides an example of one such P-box, showing the experimental uncertainty in a locally determined dynamic pressure (QC10) for one of many tunnel conditions evaluated during the Flow Survey experiment. This particular condition (Condition 27) has nominal set points of Mach number = 3.85, Reynolds number = 3x10(exp 6) ft(exp -1), total pressure = 5160 psfa, and total temperature = 150°F. The uncertainty in dynamic pressure at this condition is now well defined by this plot. Armed with this experimental data uncertainty, meaningful comparisons can be made with computational results, when their associated uncertainties are also considered.

Heather P Houlden↗

CFD Modeling of Phase Change and Pressure Drop during Violent Sloshing of Cryogenic Fluid in a Small-Scale Tank

This paper presents the development of a two-phase CFD model that is used to study sloshing of a cryogenic fluid (LN2) undergoing phase change generated by lateral acceleration and its effects on the ensuing heat and mass transfer and pressure drop in a small-scale tank. In this model the interface is captured using the volume-of-fluid (VOF) method [1]. Kinetic-theory based Schrage relation [2] is used for calculating phase change mass transfer at the liquid-vapor interface. Computational results are compared to the data provided by a non-isothermal sloshing experiment with phase change conducted by the University of Tokyo and JAXA in 2018 using liquid Nitrogen in a transparent tank. The effect of different thicknesses of a liquid temperature stratification layer created at the interface on the pressure drop during sloshing was studied experimentally. CFD results for interface movement and tank pressure are presented and compared in this paper to experimental data for the cases with the smallest and largest liquid temperature stratification thickness.

phase change↗

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can be used to either improve the model through calibration or quantify the model-form uncertainty. This work focuses on estimating the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given the uncertainty in both the simulation and experimental data. A procedure that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs is outlined. The work also assesses the impact of using different high-fidelity solvers and different turbulence models. The goal of performing this comparison is to provide a quantifiable measurement of the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis is also performed to get a model-form uncertainty for this type of analysis. The two high-fidelity CFD solvers used here are VULCAN-CFD and FUN3D, and the low-fidelity results comes from Cart3D.The experimental data comes from the 20-inch Mach 6 Tunnel located at NASA Langley Research Center. The work here shows that the models used tend to under-predict the drag and moment aerodynamic coefficients for this type of design leading to a model-form uncertainty estimate which can account for up to an 86% of the total uncertainty in the model predictions.

Laura M White↗

Multi-point Jig Twist Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft and High-Fidelity CFD Validation

This paper presents a multi-point jig twist optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft using an in-house developed aero-structural analysis solver VSPAERO coupled to BEAM3D. A vortex-lattice model of the TTBW model is developed, and a transonic small disturbance /integral boundary layer correction method is implemented in the VSPAERO model to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. Also, a structural dynamic finite-element model of the TTBW aircraft is developed. This finite-element model includes the geometric nonlinear effect due to the tension in the struts which causes a deflection-dependent nonlinear stiffness. The VSPAERO model is coupled to the corresponding finite-element model to provide a rapid aero-structural analysis. The VSPAERO model has been validated against wind tunnel test data and FUN3D CFD simulation data. Three different flight conditions corresponding to Mach 0.745 at 40,000 ft are selected for the TTBW aircraft jig twist optimization to reduce the drag coefficients at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized jig twist achieves a drag reduction of about 9 counts, 7 counts, and8 counts corresponding to the lift coefficients 0.68, 0.73, and 0.78, respectively. A high-fidelity CFD solver FUN3D is used to validate the jig twist optimization.

Multi-point↗

Atmospheric Boundary-Layer and Flutter Computations Using CFD Model of the Transonic Dynamics Tunnel

This paper presents two Computational Fluid Dynamics (CFD) models of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The TDT is a continuous-flow, closed circuit wind tunnel with a 16- by 16-foot slotted test section with cropped corners. The tunnel was originally built as the 19-ft Pressure Tunnel in 1938, but it was converted to the current transonic tunnel in the 1950s, with capabilities to use either air or heavy gas as the test medium. The first computational model describes the generation of an atmospheric-boundary-layer (ABL) profile inside the tunnel. An ABL, which includes both a wind profile and turbulence content, is one of the aerodynamic characteristics affecting the occurrence of wind-induced oscillations for a launch vehicle sitting on a pad. The challenging part of this analysis was modeling the turbulent flow inside the tunnel. This is due to the special ABL-generating hardware that was installed at the entrance of the TDT test section in order to change the downstream velocity profile and to introduce velocity fluctuations into the flow. The second CFD model builds on the computational aeroelastic results that were generated in support of the second Aeroelastic Prediction Workshop (AePW) for the NASA Benchmark Supercritical Wing (BSCW) configuration. During the AePW, the wing-only configuration (classical free-air model) was analyzed. In the current study, the flutter computations were conducted on the configuration as it was mounted in the TDT during the experiment. This includes the wing attached to the splitter plate that was attached to the wind-tunnel walls. The preliminary results show that the wind-tunnel walls marginally affect flutter prediction.

Pawel Chwalowski↗

CFD-Based Frequency Domain Method for Dynamic Stability Derivative Estimation with Application to Transonic Truss-Braced Wing

This paper presents a dynamic stability estimation technique obtained from high-fidelity CFD simulations of the Mach 0.8 Transonic Truss-Braced Wing (TTBW). A series of unsteady RANS CFD simulations in FUN3D is performed on the TTBW in pitch and plunge oscillations at various reduced frequencies. The time-domain data are transformed into the frequency-domain data by Fourier series. Transfer functions of the dynamic stability derivatives are then estimated by a frequency-domain regression. The dynamic stability derivatives with respect to the angle of attack are determined by the regression of the unsteady aerodynamic coefficients for the plunge motion. The dynamic stability derivatives with respect to the pitch rate are deter-mined by the regression of the differential unsteady aerodynamic coefficients for the pitch motion upon the removal of the angle of attack contribution by the plunge motion. The steady-state stability derivatives are then compared to the results obtained from a stability analysis code VSPAERO as well as steady-state FUN3D simulations. The comparison of the steady-state stability derivatives shows excellent agreement.

Dynamic Stability Derivatives↗

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗