Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes solver”

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 379 records · Page 21

Tandem Cylinder Noise Predictions

In an effort to better understand landing-gear noise sources, we have been examining a simplified configuration that still maintains some of the salient features of landing-gear flow fields. In particular, tandem cylinders have been studied because they model a variety of component level interactions. The present effort is directed at the case of two identical cylinders spatially separated in the streamwise direction by 3.7 diameters. Experimental measurements from the Basic Aerodynamic Research Tunnel (BART) and Quiet Flow Facility (QFF) at NASA Langley Research Center (LaRC) have provided steady surface pressures, detailed off-surface measurements of the flow field using Particle Image Velocimetry (PIV), hot-wire measurements in the wake of the rear cylinder, unsteady surface pressure data, and the radiated noise. The experiments were conducted at a Reynolds number of 166 105 based on the cylinder diameter. A trip was used on the upstream cylinder to insure a fully turbulent shedding process and simulate the effects of a high Reynolds number flow. The parallel computational effort uses the three-dimensional Navier-Stokes solver CFL3D with a hybrid, zonal turbulence model that turns off the turbulence production term everywhere except in a narrow ring surrounding solid surfaces. The current calculations further explore the influence of the grid resolution and spanwise extent on the flow and associated radiated noise. Extensive comparisons with the experimental data are used to assess the ability of the computations to simulate the details of the flow. The results show that the pressure fluctuations on the upstream cylinder, caused by vortex shedding, are smaller than those generated on the downstream cylinder by wake interaction. Consequently, the downstream cylinder dominates the noise radiation, producing an overall directivity pattern that is similar to that of an isolated cylinder. Only calculations based on the full length of the model span were able to capture the complete decay in the spanwise correlation, thereby producing reasonable noise radiation levels.

Lockard, David P.↗

Aeroacoustic Simulations of Tandem Cylinders with Subcritical Spacing

Tandem cylinders are being studied because they model a variety of component level interactions of landing gear. The present effort is directed at the case of two identical cylinders with their centroids separated in the streamwise direction by 1.435 diameters. Experiments in the Basic Aerodynamic Research Tunnel and Quiet Flow Facility at NASA Langley Research Center have provided an extensive experimental database of the nearfield flow and radiated noise. The measurements were conducted at a Mach number of 0.1285 and Reynolds number of 1.66x10(exp 5) based on the cylinder diameter. A trip was used on the upstream cylinder to insure a fully turbulent flow separation and, hence, to simulate a major aspect of high Reynolds number flow. The parallel computational effort uses the three-dimensional Navier-Stokes solver CFL3D with a hybrid, zonal turbulence model that turns off the turbulence production term everywhere except in a narrow ring surrounding solid surfaces. The experiments exhibited an asymmetry in the surface pressure that was persistent despite attempts to eliminate it through small changes in the configuration. To model the asymmetry, the simulations were run with the cylinder configuration at a nonzero but small angle of attack. The computed results and experiments are in general agreement that vortex shedding for the spacing studied herein is weak relative to that observed at supercritical spacings. Although the shedding was subdued in the simulations, it was still more prominent than in the experiments. Overall, the simulation comparisons with measured near-field data and the radiated acoustics are reasonable, especially if one is concerned with capturing the trends relative to larger cylinder spacings. However, the flow details of the 1.435 diameter spacing have not been captured in full even though very fine grid computations have been performed. Some of the discrepancy may be associated with the simulation s inexact representation of the experimental configuration, but numerical and flow modeling errors are also likely contributors to the observed differences.

Lockard, David P.↗

Rotor Airloads Prediction Using Unstructured Meshes and Loose CFD/CSD Coupling

The FUN3D unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids has been modified to allow prediction of trimmed rotorcraft airloads. The trim of the rotorcraft and the aeroelastic deformation of the rotor blades are accounted for via loose coupling with the CAMRAD II rotorcraft computational structural dynamics code. The set of codes is used to analyze the HART-II Baseline, Minimum Noise and Minimum Vibration test conditions. The loose coupling approach is found to be stable and convergent for the cases considered. Comparison of the resulting airloads and structural deformations with experimentally measured data is presented. The effect of grid resolution and temporal accuracy is examined. Rotorcraft airloads prediction presents a very substantial challenge for Computational Fluid Dynamics (CFD). Not only must the unsteady nature of the flow be accurately modeled, but since most rotorcraft blades are not structurally stiff, an accurate simulation must account for the blade structural dynamics. In addition, trim of the rotorcraft to desired thrust and moment targets depends on both aerodynamic loads and structural deformation, and vice versa. Further, interaction of the fuselage with the rotor flow field can be important, so that relative motion between the blades and the fuselage must be accommodated. Thus a complete simulation requires coupled aerodynamics, structures and trim, with the ability to model geometrically complex configurations. NASA has recently initiated a Subsonic Rotary Wing (SRW) Project under the overall Fundamental Aeronautics Program. Within the context of SRW are efforts aimed at furthering the state of the art of high-fidelity rotorcraft flow simulations, using both structured and unstructured meshes. Structured-mesh solvers have an advantage in computation speed, but even though remarkably complex configurations may be accommodated using the overset grid approach, generation of complex structured-mesh systems can require months to set up. As a result, many rotorcraft simulations using structured-grid CFD neglect the fuselage. On the other hand, unstructured-mesh solvers are easily able to handle complex geometries, but suffer from slower execution speed. However, advances in both computer hardware and CFD algorithms have made previously state-of-the-art computations routine for unstructured-mesh solvers, so that rotorcraft simulations using unstructured grids are now viable. The aim of the present work is to develop a first principles rotorcraft simulation tool based on an unstructured CFD solver.

Biedron, Robert T.↗

A Numerical Study of Anti-Vortex Film Cooling Designs at High Blowing Ratio

A concept for mitigating the adverse effects of jet vorticity and liftoff at high blowing ratios for turbine film cooling flows has been developed and studied at NASA Glenn Research Center. This "anti-vortex" film cooling concept proposes the addition of two branched holes from each primary hole in order to produce a vorticity counter to the detrimental kidney vortices from the main jet. These vortices typically entrain hot freestream gas and are associated with jet separation from the turbine blade surface. The anti-vortex design is unique in that it requires only easily machinable round holes, unlike shaped film cooling holes and other advanced concepts. The anti-vortex film cooling hole concept has been modeled computationally for a single row of 30deg angled holes on a flat surface using the 3D Navier-Stokes solver Glenn-HT. A modification of the anti-vortex concept whereby the branched holes exit adjacent to the main hole has been studied computationally for blowing ratios of 1.0 and 2.0 and at density ratios of 1.0 and 2.0. This modified concept was selected because it has shown the most promise in recent experimental studies. The computational results show that the modified design improves the film cooling effectiveness relative to the round hole baseline and previous anti-vortex cases, in confirmation of the experimental studies.

Heidmann, James D.↗

Multidisciplinary Optimization for Aerospace Using Genetic Optimization

In support of the ARMD guidelines NASA's Dryden Flight Research Center is developing a multidisciplinary design and optimization tool This tool will leverage existing tools and practices, and allow the easy integration and adoption of new state-of-the-art software. Optimization has made its way into many mainstream applications. For example NASTRAN(TradeMark) has its solution sequence 200 for Design Optimization, and MATLAB(TradeMark) has an Optimization Tool box. Other packages, such as ZAERO(TradeMark) aeroelastic panel code and the CFL3D(TradeMark) Navier-Stokes solver have no built in optimizer. The goal of the tool development is to generate a central executive capable of using disparate software packages ina cross platform network environment so as to quickly perform optimization and design tasks in a cohesive streamlined manner. A provided figure (Figure 1) shows a typical set of tools and their relation to the central executive. Optimization can take place within each individual too, or in a loop between the executive and the tool, or both.

Pak, Chan-gi↗

Performance Analysis of Scientific and Engineering Applications Using MPInside and TAU

In this paper, we present performance analysis of two NASA applications using performance tools like Tuning and Analysis Utilities (TAU) and SGI MPInside. MITgcmUV and OVERFLOW are two production-quality applications used extensively by scientists and engineers at NASA. MITgcmUV is a global ocean simulation model, developed by the Estimating the Circulation and Climate of the Ocean (ECCO) Consortium, for solving the fluid equations of motion using the hydrostatic approximation. OVERFLOW is a general-purpose Navier-Stokes solver for computational fluid dynamics (CFD) problems. Using these tools, we analyze the MPI functions (MPI_Sendrecv, MPI_Bcast, MPI_Reduce, MPI_Allreduce, MPI_Barrier, etc.) with respect to message size of each rank, time consumed by each function, and how ranks communicate. MPI communication is further analyzed by studying the performance of MPI functions used in these two applications as a function of message size and number of cores. Finally, we present the compute time, communication time, and I/O time as a function of the number of cores.

Saini, Subhash↗

Combustor Operability and Performance Verification for HIFiRE Flight 2

As part of the Hypersonic International Flight Research Experimentation (HIFiRE) Direct-Connect Rig (HDCR) test and analysis activity, three-dimensional computational fluid dynamics (CFD) simulations were performed using two Reynolds-Averaged Navier Stokes solvers. Measurements obtained from ground testing in the NASA Langley Arc-Heated Scramjet Test Facility (AHSTF) were used to specify inflow conditions for the simulations and combustor data from four representative tests were used as benchmarks. Test cases at simulated flight enthalpies of Mach 5.84, 6.5, 7.5, and 8.0 were analyzed. Modeling parameters (e.g., turbulent Schmidt number and compressibility treatment) were tuned such that the CFD results closely matched the experimental results. The tuned modeling parameters were used to establish a standard practice in HIFiRE combustor analysis. Combustor performance and operating mode were examined and were found to meet or exceed the objectives of the HIFiRE Flight 2 experiment. In addition, the calibrated CFD tools were then applied to make predictions of combustor operation and performance for the flight configuration and to aid in understanding the impacts of ground and flight uncertainties on combustor operation.

Storch, Andrea M.↗

Computation of UH-60A Airloads Using CFD/CSD Coupling on Unstructured Meshes

An unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids is used to compute the rotor airloads on the UH-60A helicopter at high-speed and high thrust conditions. The flow solver is coupled to a rotorcraft comprehensive code in order to account for trim and aeroelastic deflections. Simulations are performed both with and without the fuselage, and the effects of grid resolution, temporal resolution and turbulence model are examined. Computed airloads are compared to flight data.

Biedron, Robert T.↗

CFD Computations for a Generic High-Lift Configuration Using TetrUSS

Assessment of the accuracy of computational results for a generic high-lift trapezoidal wing with a single slotted flap and slat is presented. The paper is closely aligned with the focus of the 1st AIAA CFD High Lift Prediction Workshop (HiLiftPW-1) which was to assess the accuracy of CFD methods for multi-element high-lift configurations. The unstructured grid Reynolds-Averaged Navier-Stokes solver TetrUSS/USM3D is used for the computational results. USM3D results are obtained assuming fully turbulent flow using the Spalart-Allmaras (SA) and Shear Stress Transport (SST) turbulence models. Computed solutions have been obtained at seven different angles-of-attack ranging from 6 -37 . Three grids providing progressively higher grid resolution are used to quantify the effect of grid resolution on the lift, drag, pitching moment, surface pressure and stall angle. SA results, as compared to SST results, exhibit better agreement with the measured data. However, both turbulence models under-predict upper surface pressures near the wing tip region.

Pandya, Mohagna J.↗

Low Boom Configuration Analysis with FUN3D Adjoint Simulation Framework

Off-body pressure, forces, and moments for the Gulfstream Low Boom Model are computed with a Reynolds Averaged Navier Stokes solver coupled with the Spalart-Allmaras (SA) turbulence model. This is the first application of viscous output-based adaptation to reduce estimated discretization errors in off-body pressure for a wing body configuration. The output adaptation approach is compared to an a priori grid adaptation technique designed to resolve the signature on the centerline by stretching and aligning the grid to the freestream Mach angle. The output-based approach produced good predictions of centerline and off-centerline measurements. Eddy viscosity predicted by the SA turbulence model increased significantly with grid adaptation. Computed lift as a function of drag compares well with wind tunnel measurements for positive lift, but predicted lift, drag, and pitching moment as a function of angle of attack has significant differences from the measured data. The sensitivity of longitudinal forces and moment to grid refinement is much smaller than the differences between the computed and measured data.

Park, Michael A.↗

Continuing Validation of Computational Fluid Dynamics for Supersonic Retropropulsion

A large step in the validation of Computational Fluid Dynamics (CFD) for Supersonic Retropropulsion (SRP) is shown through the comparison of three Navier-Stokes solvers (DPLR, FUN3D, and OVERFLOW) and wind tunnel test results. The test was designed specifically for CFD validation and was conducted in the Langley supersonic 4 x4 Unitary Plan Wind Tunnel and includes variations in the number of nozzles, Mach and Reynolds numbers, thrust coefficient, and angles of orientation. Code-to-code and code-to-test comparisons are encouraging and possible error sources are discussed.

Schauerhamer, Daniel Guy↗

An Examination of Unsteady Airloads on a UH-60A Rotor: Computation Versus Measurement

An unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids is used to simulate the flow over a UH-60A rotor. Traditionally, the computed pressure and shear stresses are integrated on the computational mesh at selected radial stations and compared to measured airloads. However, the corresponding integration of experimental data uses only the pressure contribution, and the set of integration points (pressure taps) is modest compared to the computational mesh resolution. This paper examines the difference between the traditional integration of computed airloads and an integration consistent with that used for the experimental data. In addition, a comparison of chordwise pressure distributions between computation and measurement is made. Examination of this unsteady pressure data provides new opportunities to understand differences between computation and flight measurement.

Biedron, Robert T.↗

FUN3D Analyses in Support of the First Aeroelastic Prediction Workshop

This paper presents the computational aeroelastic results generated in support of the first Aeroelastic Prediction Workshop for the Benchmark Supercritical Wing (BSCW) and the HIgh REynolds Number AeroStructural Dynamics (HIRENASD) configurations and compares them to the experimental data. The computational results are obtained using FUN3D, an unstructured grid Reynolds-averaged Navier-Stokes solver developed at NASA Langley Research Center. The analysis results for both configurations include aerodynamic coefficients and surface pressures obtained for steady-state or static aeroelastic equilibrium (BSCW and HIRENASD, respectively) and for unsteady flow due to a pitching wing (BSCW) or modally-excited wing (HIRENASD). Frequency response functions of the pressure coefficients with respect to displacement are computed and compared with the experimental data. For the BSCW, the shock location is computed aft of the experimentally-located shock position. The pressure distribution upstream of this shock is in excellent agreement with the experimental data, but the pressure downstream of the shock in the separated flow region does not match as well. For HIRENASD, very good agreement between the numerical results and the experimental data is observed at the mid-span wing locations.

Chwalowski, Pawel↗

FUN3D Airload Predictions for the Full-Scale UH-60A Airloads Rotor in a Wind Tunnel

An unsteady Reynolds-Averaged Navier-Stokes solver for unstructured grids, FUN3D, is used to compute the rotor performance and airloads of the UH-60A Airloads Rotor in the National Full-Scale Aerodynamic Complex (NFAC) 40- by 80-foot Wind Tunnel. The flow solver is loosely coupled to a rotorcraft comprehensive code, CAMRAD-II, to account for trim and aeroelastic deflections. Computations are made for the 1-g level flight speed-sweep test conditions with the airloads rotor installed on the NFAC Large Rotor Test Apparatus (LRTA) and in the 40- by 80-ft wind tunnel to determine the influence of the test stand and wind-tunnel walls on the rotor performance and airloads. Detailed comparisons are made between the results of the CFD/CSD simulations and the wind tunnel measurements. The computed trends in solidity-weighted propulsive force and power coefficient match the experimental trends over the range of advance ratios and are comparable to previously published results. Rotor performance and sectional airloads show little sensitivity to the modeling of the wind-tunnel walls, which indicates that the rotor shaft-angle correction adequately compensates for the wall influence up to an advance ratio of 0.37. Sensitivity of the rotor performance and sectional airloads to the modeling of the rotor with the LRTA body/hub increases with advance ratio. The inclusion of the LRTA in the simulation slightly improves the comparison of rotor propulsive force between the computation and wind tunnel data but does not resolve the difference in the rotor power predictions at mu = 0.37. Despite a more precise knowledge of the rotor trim loads and flight condition, the level of comparison between the computed and measured sectional airloads/pressures at an advance ratio of 0.37 is comparable to the results previously published for the high-speed flight test condition.

Lee-Rausch, Elizabeth M.↗

Blade Displacement Predictions for the Full-Scale UH-60A Airloads Rotor

An unsteady Reynolds-Averaged Navier-Stokes solver for unstructured grids is loosely coupled to a rotorcraft comprehensive code and used to simulate two different test conditions from a wind-tunnel test of a full-scale UH-60A rotor. Performance data and sectional airloads from the simulation are compared with corresponding tunnel data to assess the level of fidelity of the aerodynamic aspects of the simulation. The focus then turns to a comparison of the blade displacements, both rigid (blade root) and elastic. Comparisons of computed root motions are made with data from three independent measurement systems. Finally, comparisons are made between computed elastic bending and elastic twist, and the corresponding measurements obtained from a photogrammetry system. Overall the correlation between computed and measured displacements was good, especially for the root pitch and lag motions and the elastic bending deformation. The correlation of root lead-lag motion and elastic twist deformation was less favorable.

Bledron, Robert T.↗

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three dimensional unsteady results were produced and then time averaged for the collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Collection efficiencies were generated for three spherical particle sizes, 100, 20, and 5 micron in diameter, using the codes LEWICE3D and LEWICE2D. The free stream Mach number was 0.27, representing a velocity of approximately 86 ms. It was observed that a reduction in velocity of about 15-20 occurred as the flow entered the shroud of the probe.Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the effect of the protective shroud.

Aircraft Icing↗

Towards Full Aircraft Airframe Noise Prediction: Detached Eddy Simulations

Results from a computational study on the aeroacoustic characteristics of an 18%-scale, semi-span Gulf-stream aircraft model are presented in this paper. NASA's FUN3D unstructured compressible Navier-Stokes solver was used to perform steady and unsteady simulations of the flow field associated with this high-fidelity aircraft model. Solutions were obtained for free-air at a Mach number of 0.2 with the flap deflected at 39 deg, with the main gear off and on (the two baseline configurations). Initially, the study focused on accurately predicting the prominent noise sources at both flap tips for the baseline configuration with deployed flap only. Building upon the experience gained from this initial effort, subsequent work involved the full landing configuration with both flap and main landing gear deployed. For the unsteady computations, we capitalized on the Detached Eddy Simulation capability of FUN3D to capture the complex time-dependent flow features associated with the flap and main gear. To resolve the noise sources over a broad frequency range, the tailored grid was very dense near the flap inboard and outboard tips and the region surrounding the gear. Extensive comparison of the computed steady and unsteady surface pressures with wind tunnel measurements showed good agreement for the global aerodynamic characteristics and the local flow field at the flap inboard tip. However, the computed pressure coefficients indicated that a zone of separated flow that forms in the vicinity of the outboard tip is larger in extent along the flap span and chord than measurements suggest. Computed farfield acoustic characteristics from a FW-H integral approach that used the simulated pressures on the model solid surface were in excellent agreement with corresponding measurements.

Khorrami, Mehdi R.↗

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Water Content↗