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 55 records · Page 3

Spanwise Wing Loads on the Space Shuttle Orbiter during Roll Maneuver

Spanwise aerodynamic loads for the low-Mach, high-attitude portion of ascent for the Space Shuttle Orbiter are presented. In this Mach 0.3 flight regime, also called the roll maneuver, pre-stall and post-stall distributions of aerodynamic wing shear force, bending moment, and torsion moment were obtained from wind tunnel test data and computational fluid dynamics simulations of the Space Shuttle Launch Vehicle. The spanwise loads were computed by integration of surface pressure data. The existing historical operational database of spanwise wing loads for the Orbiter does not cover this low-Mach, high-attitude condition, however for Mach 0.6 low-attitude conditions the experimental and computational results compare well with the operational data which has been validated by past flight measurements. Spanwise load distributions exhibit typical delta-wing characteristics. The computational results capture well the peak loading condition in the pre-stall case, but show more load relief for the post-stall case than was observed in the wind tunnel test data.

Doggett, Glen P.↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Grid Sensitivity Study for Slat Noise Simulations

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations in conjunction with a Ffowcs Williams-Hawkings acoustics solver. Many previous simulations have been performed for the configuration, and the case was introduced as a new category for the Second AIAA workshop on Benchmark problems for Airframe Noise Configurations (BANC-II). However, the cost of the simulations has restricted the study of grid resolution effects to a baseline grid and coarser meshes. In the present study, two different approaches are being used to investigate the effect of finer resolution of near-field unsteady structures. First, a standard grid refinement by a factor of two is used, and the calculations are performed by using the same CFL3D solver employed in the majority of the previous simulations. Second, the OVERFLOW code is applied to the baseline grid, but with a 5th-order upwind spatial discretization as compared with the second-order discretization used in the CFL3D simulations. In general, the fine grid CFL3D simulation and OVERFLOW calculation are in very good agreement and exhibit the lowest levels of both surface pressure fluctuations and radiated noise. Although the smaller scales resolved by these simulations increase the velocity fluctuation levels, they appear to mitigate the influence of the larger scales on the surface pressure. These new simulations are used to investigate the influence of the grid on unsteady high-lift simulations and to gain a better understanding of the physics responsible for the noise generation and radiation.

Lockard, David P.↗

Influence of Spanwise Boundary Conditions on Slat Noise Simulations

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations with the OVERFLOW code in conjunction with a Ffowcs Williams-Hawkings acoustics solver. In the present study, two different spanwise grids are being used to investigate the effect of the spanwise extent and periodicity on the near-field unsteady structures and radiated noise. The baseline grid with periodic boundary conditions has a short span equal to 1/9th of the stowed chord, whereas the other, longer span grid adds stretched grids on both sides of the core, baseline grid to allow inviscid surface boundary conditions at both ends. The results indicate that the near-field mean statistics obtained using the two grids are similar to each other, as are the directivity and spectral shapes of the radiated noise. However, periodicity forces all acoustic waves with less than one wavelength across the span to be two-dimensional, without any variation in the span. The spanwise coherence of the acoustic waves is what is needed to make estimates of the noise that would be radiated from realistic span lengths. Simulations with periodic conditions need spans of at least six slat chords to allow spanwise variation in the low-frequencies associated with the peak of broadband slat noise. Even then, the full influence of the periodicity is unclear, so employing grids with a fine, central region and highly stretched meshes that go to slip walls may be a more efficient means of capturing the spanwise decorrelation of low-frequency acoustic phenomena.

Lockard, David P.↗

Observations of Keyhole Porosity and Comparisons to Analytical Models for Ti-6Al-4V Powder Bed Fusion

Keyhole porosity defects are a common concern in powder bed fusion (PBF) processing. Keyhole porosity prediction models have generally fallen into two categories – high fidelity computational fluid dynamics simulations and low fidelity analytical or empirical models based on keyhole vibration dynamics. The low computational cost of low fidelity models allows them to better approach part-scale predictions. However, studies on low fidelity modeling techniques are limited by the lack of experimental data to assess the validity of the calibration over a wide range of processing conditions. This work extracts and quantifies keyhole porosity across 14 laser velocities, two laser powers, and six+ repetitions for a total of 176 independent trials. The measured porosity data are compared to low fidelity keyhole models to assess their success in predictive porosity occurrence. This work impacts the field by providing independent validation of keyhole porosity models for PBF for use in part-scale defect prediction.

powder bed fusion↗

Extension of Miles Equation for Ring Baffle Damping Predictions to Small Slosh Amplitudes and Large Baffle Widths

The Miles equation has long been used to predict slosh damping in liquid propellant tanks due to ring baffles. The original work by Miles identifies defined limits to its range of application. Recent evaluations of the Space Launch System identified that the Core Stage baffle designs resulted in violating the limits of the application of the Miles equation. This paper describes the work conducted by NASA/MSFC to develop methods to predict slosh damping from ring baffles for conditions for which Miles equation is not applicable. For asymptotically small slosh amplitudes or conversely large baffle widths, an asymptotic expression for slosh damping was developed and calibrated using historical experimental sub-scale slosh damping data. For the parameter space that lies between region of applicability of the asymptotic expression and the Miles equation, Computational Fluid Dynamics simulations of slosh damping were used to develop an expression for slosh damping. The combined multi-regime slosh prediction methodology is shown to be smooth at regime boundaries and consistent with both sub-scale experimental slosh damping data and the results of validated Computational Fluid Dynamics predictions of slosh damping due to ring baffles.

West, Jeff↗

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

Extraction of Dynamic Inflow Models for Coaxial and Tandem Rotors from CFD Simulations

The dynamic inflow coupling with rotor/body dynamics is crucial in the analysis of stability and control law design for helicopters. Over the past several decades, finite-state inflow models for single rotor configurations in hover, forward flight, and maneuver have developed (Ref.1-3). By capturing the interference effects between rotors, the extension of pressure potential finite state inflow model has promising result for coaxial rotor configuration (Ref.4-6). Recently, the focus of the dynamic inflow modeling has shifted to tandem rotor configurations (Ref.7, 8). The development of the dynamic inflow models for tandem rotor configuration still have some limitations due to the lack of knowledge of rotor-to-rotor interference, and rotor-wake interference. Experimental methods, and computational fluid dynamics methods are commonly used to understand the rotor performance and rotor airload variations, and measure or predict inflow velocity distributions at the rotor desk. The inflow distributions are subsequently used to improve the dynamic inflow models. Tandem rotor configurations have been studied experimentally and computationally for several decades (Ref.9-12). Sweet (Ref.10) observed that a tandem rotor with 76-percent-radius overlap required 14% more induced power at hovering condition, relative to an isolated rotor of equivalent disk area. Sweet also found that, above a shaft-to-shaft distance of 1.03 diameter, the performance of the tandem rotor was nearly the same as two isolated rotors. The objective of the present study is to apply computational fluid dynamics simulations of tandem rotors for the extraction of dynamic inflow models. The extended methodology is first validated by comparing the computed induced power against test data. Subsequently inflow distributions and wake structures are analyzed.

Extraction of Dynamic Inflow Models↗

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

Goekcen, Tahir↗

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

CFD↗

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

Visualization of a Sweeping Jet by Laser Speckle Retro-reflective Background Oriented Schlieren

Sweeping jet actuators are currently being studied as a method of active flow control to increase the performance of wings and tails. Computational Fluid Dynamics simulations show that the jet does not sweep at a constant rate, instead spending more time at the minimum and maximum of its extent. As part of a test done to provide dynamic flow measurements for Computation Fluid Dynamics validation in NASA Ames Fluid Mechanics Lab, Particle Image Velocimetry (PIV) was performed on a single sweeping jet actuator designed by NASA and Boeing that was inset in a small airfoil. However, even at low jet pressures, the PIV could not be acquired at high enough frequencies to both spatially and temporally resolve the stickiness of the jet motion. Therefore, a small portion of the PIV laser was picked off and used to acquire flow visualization at higher frequencies through a new technique, laser speckle Retroreflective Background Oriented Schlieren (RBOS).

sweeping jet actuator↗

Computational Simulations of the 10-MW TP3 Arc-Jet Facility Flow

This paper reports computational simulations and analysis in support of calibration tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using stagnation calorimeters and two different blunted wedge models with calibration plates at a wide range of conditions in the NASA Ames 10-megawatt TP3 (Test Position 3) facility. Data were obtained using four different conical nozzles with the same test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot and null-point heat flux probes were also performed at several arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. These computational simulations provide estimates of the arc-jet test environment parameters that are not measured but are needed to evaluate the performance of thermal protection system materials, along with further valuable insights into the arc-jet testing environment. Simulation results are used to estimate centerline total enthalpy, surface shear, boundary layer thickness, and boundary layer edge Mach number and to verify that specific test requirements from the Orion program are met.

CFD↗

CFD Simulations of the Truncated IHF Semi-Elliptical Nozzle Flow: Convective Calibration Plate

This paper reports computational simulations of a high-enthalpy panel test facility at NASA Ames Research Center. The recently developed facility configuration uses a truncated version of an existing semi-elliptical nozzle in the 60-MW Interaction Heating Facility, and it provides testing capability for panel test articles as large as 43.2 cm x 43.2 cm (17 in x 17 in). A water-cooled calibration plate, which is attached to the semi-elliptical nozzle bottom surface and instrumented with twenty Gardon-gage heat flux calorimeters and twelve pressure gages, is used to characterize the convective heating environment over the panel test articles. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, and flowfield over panel test articles. Both laminar and turbulent simulations are performed, and the computational results are compared with the calibration plate measurements. Comparisons of computed and measured heat fluxes indicate that the boundary layer in the semi-elliptical nozzle is transitional at certain conditions (relatively higher arc-heater pressures or mass flow rates).

Arc-jets↗