Results for “Fluid drag” — Page 21

DOE OSTI2025

Small-scale validation tests for MFIX-Exa CFD-DEM

This report describes several bench-scale fluidization experiments that can be used to validate the CFD-DEM method as encapsulated in the MFIX-Exa code. The five cases considered are the cold-flow fluidized beds of Müller et al., Link et al. (spout-fluid), and Goldschmidt et al. (bi-disperse), the hot fluidized bed of Patil et al., and the adsorbing fluidized bed of Li et al. and Janssen. In most cases, MFIX-Exa with “standard” or “typical” CFD-DEM settings, the Gidaspow drag model, and the Gunn heat transfer provide a relatively good prediction of the quantities considered: mean void fraction profiles, mean velocity profiles, fluctuating velocity profiles, mean particle temperature and segregation index. These results, with other verification and validation tests reported elsewhere, contribute to a body of work providing confidence and credibility in CFD predictions from the MFIX-Exa code.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2003

Development of an Output-based Adaptive Method for Multi-Dimensional Euler and Navier-Stokes Simulations

The use of computational simulations in the prediction of complex aerodynamic flows is becoming increasingly prevalent in the design process within the aerospace industry. Continuing advancements in both computing technology and algorithmic development are ultimately leading to attempts at simulating ever-larger, more complex problems. However, by increasing the reliance on computational simulations in the design cycle, we must also increase the accuracy of these simulations in order to maintain or improve the reliability arid safety of the resulting aircraft. At the same time, large-scale computational simulations must be made more affordable so that their potential benefits can be fully realized within the design cycle. Thus, a continuing need exists for increasing the accuracy and efficiency of computational algorithms such that computational fluid dynamics can become a viable tool in the design of more reliable, safer aircraft. The objective of this research was the development of an error estimation and grid adaptive strategy for reducing simulation errors in integral outputs (functionals) such as lift or drag from from multi-dimensional Euler and Navier-Stokes simulations. In this final report, we summarize our work during this grant.

Darmofal, David L.↗

NASA NTRS2022

Comparison of Computational Predictions of the Mach 0.80 Transonic Truss-Braced Wing Configuration with Experimental Data

The NASA Langley and Ames Research Centers have teamed together for comparisons of computational predictions of the Boeing Mach 0.80Transonic Truss-Braced Wing (TTBW) configuration with a high-speed experimental dataset. The Mach 0.80 TTBW vehicle is a high wing, high aspect ratio configuration, designed for a high lift-to-drag ratio. System studies have predicted significant fuel burn and emissions benefits with the TTBW technology moving toward meeting NASA Subsonic Transport Systems-Level-Metrics. A 4.5% scale Mach 0.80 design TTBW model was recently tested at the NASA Ames Research Center 11-by 11-Foot Transonic Wind Tunnel(11-Ft TWT),providing a valuable dataset to validate computational tools and investigate best practices as risk reduction efforts continue for the development of the advanced TTBW vehicle. The NASA Computational Fluid Dynamics (CFD) team has computed free-air flow solutions on the Mach 0.80 design flight configuration and two wind tunnel configuration variants using the USM3D and LAVA flow solvers. Accurate modeling of the configuration tested in the wind tunnel environment is critical to validating the CFD tools, thus the team has included the internal cavity region and sting in their modeling of the configuration, similar to that tested in the 11-Ft TWT. Overall, the CFD simulations compared well and show similar trends as the corrected experimental data for lift and drag polars. The CFD predicted lift curve is shifted in angle of attack from what was observed in the experiment. The shift in lift also was seen in the pitching moment comparison plots. CFD solutions were computed at constant CL test point values and showed overall very good agreement when comparing constant spanwise cuts of pressure coefficient data on the wing and strut with experimental data. CFD cavity corrections were also investigated using the 11-Ft TWT cavity correction method, similar to that used to correct the wind tunnel data. Results showed some improvement in pitching moment coefficient predictions, and an increase in drag, shifting the data to the right in drag polars, further from the experimental data at lower lift conditions, good agreement near the design CL, and slight improvement at the higher lift conditions.

Transonic Truss-Braced Wing↗

NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗

NASA NTRS2013

A Non-Cut Cell Immersed Boundary Method for Use in Icing Simulations

This paper describes a computational fluid dynamic method used for modelling changes in aircraft geometry due to icing. While an aircraft undergoes icing, the accumulated ice results in a geometric alteration of the aerodynamic surfaces. In computational simulations for icing, it is necessary that the corresponding geometric change is taken into consideration. The method used, herein, for the representation of the geometric change due to icing is a non-cut cell Immersed Boundary Method (IBM). Computational cells that are in a body fitted grid of a clean aerodynamic geometry that are inside a predicted ice formation are identified. An IBM is then used to change these cells from being active computational cells to having properties of viscous solid bodies. This method has been implemented in the NASA developed node centered, finite volume computational fluid dynamics code, FUN3D. The presented capability is tested for two-dimensional airfoils including a clean airfoil, an iced airfoil, and an airfoil in harmonic pitching motion about its quarter chord. For these simulations velocity contours, pressure distributions, coefficients of lift, coefficients of drag, and coefficients of pitching moment about the airfoil's quarter chord are computed and used for comparison against experimental results, a higher order panel method code with viscous effects, XFOIL, and the results from FUN3D's original solution process. The results of the IBM simulations show that the accuracy of the IBM compares satisfactorily with the experimental results, XFOIL results, and the results from FUN3D's original solution process.

Sarofeen, Christian M.↗

NASA NTRS2015

Application of Pinniped Vibrissae to Aeropropulsion

Vibrissae of Phoca Vitulina (Harbor Seal) and Mirounga Angustirostris (Elephant Seal) possess undulations along their length. Harbor Seal Vibrissae were shown to reduce vortex induced vibrations and reduce drag compared to appropriately scaled cylinders and ellipses. Samples of Harbor Seal vibrissae, Elephant Seal vibrissae and California Sea Lion vibrissae were collected from the Marine Mammal Center in California. CT scanning, microscopy and 3D scanning techniques were utilized to characterize the whiskers. Computational fluid dynamics simulations of the whiskers were carried out to compare them to an ellipse and a cylinder. Leading edge parameters from the whiskers were used to create a 3D profile based on a modern power turbine blade. The NASA SW-2 facility was used to perform wind tunnel cascade testing on the 'Seal Blades'. Computational Fluid Dynamics simulations were used to study incidence angles from -37 to +10 degrees on the aerodynamic performance of the Seal Blade. The tests and simulations were conducted at a Reynolds number of 100,000. The Seal Blades showed consistent performance improvements over the baseline configuration. It was determined that a fuel burn reduction of approximately 5 could be achieved for a fixed wing aircraft. Noise reduction potential is also explored.

wakes↗

NASA NTRSDate not supplied

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

NASA NTRSDate not supplied

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

NASA NTRSDate not supplied

Sensitivity Study of Dust-Induced Surface Erosion during Martian Planetary Entry

Large, global dust storms that periodically form in the Martian atmosphere pose a threat to NASA missions due to a lack of understanding of how they will affect the heat loading and mass loss of the entry vehicle’s thermal protection system (TPS). Previous studies have laid out viable methodologies for estimating the amount of surface recession from a significant dust storm, but rely on underlying modeling assumptions that are difficult to validate. The current work utilizes the DUST library, which is a Lagrangian particle solver under development at NASA Ames Research Center. DUST includes variable-fidelity fluid-particle coupling and has been tested for a wide range of problem types. The DUST library is used in conjunction with the US3D Naiver-Stokes solver. Using these tools, a sensitivity study is performed to determine the influence of the underlying dust modeling (drag coefficient, Nusselt number, and size distribution) on the surface erosion through a hypersonic Martian entry trajectory. The surface erosion was found to be insensitive to the drag and heating models, but was affected significantly by the size distributions parameters. The total recession was comparable to the thermochemical ablation predicted by material response codes, but the lack of a particle erosion model for PICA limits the predictive capability of the current work.

Thermal Protection Systems↗

NASA NTRSDate not supplied

Fluid-Structure Interaction Simulations Of Supersonic Parachute Inflation: Model Sensitivities

Fluid-structure interaction simulations have the potential to accelerate the certification process of parachute systems for interplanetary robotic exploration missions while helping to mitigate risks by covering a wider range of scenarios than can be flight tested. The Launch, Ascent, and Vehicle Aerodynamics team is working toward actualizing this potential by further developing the capability to perform such fluid-structure interaction simulations and validating it with the best supersonic parachute flight test data available. We perform several simulations of the Advanced Supersonic Parachute Inflation Research Experiments’ third flight test (SR03) where we independently vary several parameters to determine their sensitivity on the predicted drag. We investigate sensitivity of the drag to the added mass of the confluence fitting, to the effective porosity, to how the gore seams are modeled, to viscous effects, and to deceleration. Modeling the added mass of the confluence fitting has virtually no impact on the drag, but is worthwhile to include to better match the overall system’s mass, which is important for the simulations with deceleration. Parachute drag is only marginally sensitive to porosity and viscous effects. Maximum parachute drag is most sensitive to how the gore seams are modeled, causing a nearly 20% drop when they are modeled as stiffer radials, and to deceleration, causing an additional 10% reduction. Results show improved agreement with flight test measurements when the FSI model is closest to the as-built configuration and physics experienced in flight: with the confluence fitting, viscous effects, the gore seams modeled as stiff radials, and including deceleration.

ESM↗

NASA NTRSDate not supplied

High-Lift Prediction Workshop 5: Summary of Reynolds-Averaged Navier–Stokes Technology Focus Group Summary

The fifth High-Lift Prediction Workshop (HLPW-5), which involved the high-lift version of the NASA common research model in several configurations, assessed various computational fluid dynamics methods, including Reynolds-averaged Navier-Stokes (RANS) and hybrid large-eddy simulations. This paper summarizes RANS solutions computed on fixed grids. Case 1, a verification case, considered a simple wing-body configuration and focused on grid convergence of lift, drag, and pitching moment coefficients. Case 2 is a configuration buildup case that focused on predicting the effects of increasing geometric complexity. For buildup configurations with slats, flaps, and a nacelle/pylon (configurations 2.2, 2.3, and 2.4), wind-tunnel data were provided by ONERA. Case 3 focused on the reference landing configuration at four Reynolds-number conditions. For Case 1, grid-converged RANS solutions were achieved using the standard Spalart-Allmaras (SA) turbulence model and the SA model with a quadratic constitutive relation and a rotation correction. Agreement between RANS solutions was observed for the simplest configuration 2.1 of Case 2 with the standard SA model. For other configurations, agreement between RANS solutions was hampered by insufficient iterative and grid convergence, especially at high angles of attack. In comparison with the experiment, RANS solutions qualitatively showed the correct configuration buildup trend but underpredicted lift and overpredicted both drag and pitching moment at high angles of attack.

Boris Diskin↗

NASA NTRS1993

Interpretation of waverider performance data using computational fluid dynamics

A computational study was conducted to better understand experimental results obtained from wind tunnel tests of a Mach 4 waverider model and a comparative reference configuration. The experimental results showed that the performance of the reference configuration was slightly better than that of the waverider model. These results contradict waverider design theory, which suggests that a waverider optimized for maximum lift-to-drag should provide better performance than any other non-waverider configuration at a given design point, especially at hypersonic speeds. The computational results showed that the predicted surface pressure values and the integrated lift and drag coefficients from the pressure distributions were much lower for the reference model than for the flat-top model, due to the reference model bottom surface having a slight expansion. The lift-to-drag ratios for the flat-top model were higher due to a relatively low drag for the same amount of lift. These results indicate that the performance advantage of the reference model was due to the shape of the bottom surface and not due to the flat top surface. The results also showed that the reference model exhibited the same shock attachment characteristics as the waverider because the planform shapes were identical. CFD predictions show that the planform shape gives the waverider an advantage in performance over conventional hypersonic vehicles and that altering the bottom surface of a waverider does not cause significant performance degradation.

Cockrell, Charles E., Jr.↗

NASA NTRS1993

Interpretation of Waverider Performance Data Using Computational Fluid Dynamics

A computational study was conducted to better understand experimental results obtained from wind tunnel tests of a Mach 4 waverider model and a comparative reference configuration. The experimental results showed that the performance of the reference configuration was slightly better than that of the waverider model. These results contradict waverider design theory, which suggests that a waverider optimized for maximum lift-to-drag should provide better performance than any other non-waverider configuration at a given design point, especially at hypersonic speeds. The computational results showed that the predicted surface pressure values and the integrated lift and drag coefficients from the pressure distributions were much lower for the reference model than for the flat-top model, due to the reference model bottom surface having a slight expansion. The lift-to-drag ratios for the flat-top model were higher due to a relatively low drag for the same amount of lift. These results indicate that the performance advantage of the reference model was due to the shape of the bottom surface and not due to the flat top surface. The results also showed that the reference model exhibited the same shock attachment characteristics as the waverider because the planform shapes were identical. CFD predictions show that the planform shape gives the waverider an advantage in performance over conventional hypersonic vehicles and that altering the bottom surface of a waverider does not cause significant performance degradation.

Cockrell, Charles E., Jr.↗

NASA NTRS1987

Typical testing experience in cryogenic wind tunnels

Airfoil testing as a major activity of many wind tunnels is considered, and emphasis is placed on airfoil testing with sidewall boundary-layer removal, oscillating airfoil testing, and three-dimensional wake topology. Buffeting as the structural response of an aircraft or test model to the aerodynamic excitation produced by separated flows is analyzed, as well as flutter problems and testing techniques. Strain-gage balance studies, optical methods and visualization techniques, and fluid mechanics are examined. Focus is placed on a two-spot laser as a boundary-layer probe, hot-film gages, fluctuating-pressure measurements, skin-friction balance, and nonadiabatic airfoil testing. Boattail afterbody pressure tests are covered, and Space Shuttle Orbiter base drag and launch configuration are discussed.

Kilgore, Robert A.↗

NASA NTRS2011

Flow Simulation of Supersonic Inlet with Bypass Annular Duct

A relaxed isentropic compression supersonic inlet is a new concept that produces smaller cowl drag than a conventional inlet, but incurs lower total pressure recovery and increased flow distortion in the (radially) outer flowpath. A supersonic inlet comprising a bypass annulus to the relaxed isentropic compression inlet dumps out airflow of low quality through the bypass duct. A reliable computational fluid dynamics solution can provide considerable useful information to ascertain quantitatively relative merits of the concept, and further provide a basis for optimizing the design. For a fast and reliable performance evaluation of the inlet performance, an equivalent axisymmetric model whose area changes accounts for geometric and physical (blockage) effects resulting from the original complex three-dimensional configuration is proposed. In addition, full three-dimensional calculations are conducted for studying flow phenomena and verifying the validity of the equivalent model. The inlet-engine coupling is carried out by embedding numerical propulsion system simulation engine data into the flow solver for interactive boundary conditions at the engine fan face and exhaust plane. It was found that the blockage resulting from complex three-dimensional geometries in the bypass duct causes significant degradation of inlet performance by pushing the terminal normal shock upstream.

Kim, HyoungJin↗

NASA NTRSDate not supplied

Pumping Loss of Shrouded Meshed Spur Gears

High speed rotorcraft transmissions are subject to load-independent power losses consisting of drag loss and pumping loss. Tightly conforming shrouds enclosing the transmission gears are often incorporated to reduce the drag component of the total load-independent losses. However, tightly conforming axial shrouds can result in an increase in the pumping loss component. Quantifying the pumping loss of shrouded gear transmissions has been the subject of many studies. This study presents a new approach for estimating pumping loss based on the concept of swept volume and examines the applicability of the approach to various shroud configurations.The drag loss and pumping loss of a shrouded spur gear pair has been determined through testing using the NASAGlenn Research Center (GRC) Gear Windage Test Facility.The results from this testing have been compared to theoretical results using the formulations presented in this study. In addition, computational fluid dynamic (CFD) analysis has been conducted for the various shroud configurations tested at NASA GRC. The results from the CFD analysis confirm the theoretical and empirical results and provide insight into the applicability of the swept volume approach for estimating pumping power loss of shrouded gear transmissions.

Michael J. Hurrell↗

NASA NTRSDate not supplied

Eigenvalue Sensitivity Computations for Linear Stability Theory

To realize the drag reduction benefit of boundary-layer transition control strategies, it is crucial to integrate transition prediction into the vehicle design through an optimization process. The integration of transition prediction based on linear stability analysis into adjoint d design optimization requires coupling an adjoint enabled computational fluid dynamics (CFD) solver with an adjoint enabled linear stability code. In particular, the boundary-layer transition location is often predicted using the N-factor method based on linear stability theory (LST). Thus, sensitivity of the linear-stability eigenvalues constitute an essential building block for optimizing the laminar flow performance. The present paper describes an implementation of LST eigenvalue sensitivity analysis that can be easily coupled with a CFD solver. Specifically, we describe a discrete adjoint formulation for the transition location prediction based on the N-factor method. The verification of this formulation is carried out by comparing the adjoint-based sensitivity of the local growth rate of a given instability mode with respect to the disturbance frequency, and the adjoint-based sensitivity of the transition location with respect to spanwise wavenumber with those sensitivities computed using a finite-difference approximation. Finally, the adjoint LST formulation is applied to flat-plate boundary-layer flows at transonic, supersonic, and hypersonic conditions, to determine the behavior and sensitivities of the transition location with respect to a range of disturbance spanwise wavenumbers.

Boundary Layer Transition↗

DOE OSTI2022

Role of Three-Dimensional Swirl in Forced Convection Heat Transfer Enhancement in Wavy-Plate-Fin Channels

The influence of wall-corrugation-induced swirl flow on enhanced forced convection in wavy-plate-fin cores has been investigated. Three-dimensional computational simulations were carried out for steady-state periodically developed air flow (Pr ~ 0.71; 50 ≤ Re ≤ 4000) with channel walls subject to constant-uniform temperature conditions. The recirculation that develops in the wall troughs and grows to have an axially helical character is scaled by the Swirl number Sw. As Sw increases with higher flowrate and/or corrugation severity, tornado-shaped vortices appear in the wave trough region midway of the interfin channel height, then extend longitudinally to encompass majority of the flow channel. The local wall-shear and heat transfer coefficient variations indicate that boundary-layer thinning upstream of the wave peak aids in intensifying momentum and heat transfer. However, the flow recirculation in wall trough impedes heat transfer at low Sw due to flow stagnation but promotes it at high Sw because of the vortices-induced augmented fluid mixing. The effects of this secondary flow are quantified by Φf(or j), which is seen to increase log-linearly as fin corrugation aspect ratio γ and/or fin spacing ratio ζ increases; the influence of cross section aspect ratio α is marginal. Moreover, the pressure drag penalty due to swirl critically affects overall pressure loss, and its proportion remains nearly constant when α varies, but grows as Sw, γ, and/or ζ increases and can be as much as 80% of the total pressure drop.

Engineering↗