Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational fluid dynamics”

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 721 records · Page 40

Prediction of BVI Noise for an Active Twist Rotor Using a Loosely Coupled CFD/CSD Method and Comparison to Experimental Data

Numerical predictions of the acoustic characteristics of an Active Twist Rotor (ATR), using two methods to compute the rotor blade aerodynamics and elastic blade motion are compared to experimental data from a wind tunnel test in the NASA Langley Transonic Dynamics Tunnel (TDT) in 2000. The first method, a loosely coupled iterative method, utilizes the Computational Fluid Dynamics (CFD) code OVERFLOW 2 and the Computational Structural Dynamics (CSD) code CAMRAD II. The second method utilizes the CAMRAD II free-wake model only. The harmonic active-twist control to the main rotor blade system is identified with three parameters - harmonic actuation frequency, actuation amplitude, and control phase angle. The resulting aerodynamics and blade motion data from the two methods are then used in the acoustics code PSU-WOPWOP to predict acoustic pressure on a spherical array of equally spaced observers surrounding the rotor. This spherical distribution of pressure is used to compute the sound power level representing baseline and actuated conditions. Sound power levels for three categories of noise are defined as - blade-vortex interaction sound power level (BVIPWL), low frequency sound power level (LFPWL), and overall sound power level, OAPWL. Comparisons with measured data indicate the CFD/CSD analysis successfully captures the trends in sound power levels and the effects of active-twist control at advance ratios of 0.14 and 0.17. The free-wake model predictions show inconsistent sound power levels relative to the trends in the experimental and CFD data. This paper presents the first ever comparison between CFD/CSD acoustic predictions for an active-twist rotor and experimental measurements.

Fogarty, David E.↗

Unstructured Finite Volume Computational Thermo-Fluid Dynamic Method for Multi-Disciplinary Analysis and Design Optimization

This paper describes a finite volume computational thermo-fluid dynamics method to solve for Navier-Stokes equations in conjunction with energy equation and thermodynamic equation of state in an unstructured coordinate system. The system of equations have been solved by a simultaneous Newton-Raphson method and compared with several benchmark solutions. Excellent agreements have been obtained in each case and the method has been found to be significantly faster than conventional Computational Fluid Dynamic(CFD) methods and therefore has the potential for implementation in Multi-Disciplinary analysis and design optimization in fluid and thermal systems. The paper also describes an algorithm of design optimization based on Newton-Raphson method which has been recently tested in a turbomachinery application.

Majumdar, Alok↗

Two-Way Coupled CFD/DEM Simulation of Electromagnetic Heating in a Fluidized Bed Reactor

A set of tools has been developed to couple a multiphase reacting computational fluid dynamics solver (MFIX) with two computational electromagnetic solvers (COMSOL and Elmer). In contrast to several recent investigations, which calculate a single static EM field prior to the CFD calculation, in the present work the electromagnetic field is periodically updated during the simulation. This enables the EM field and associated EM heating predictions to reflect the effects of evolving material properties (temperature, composition, size etc.) and the evolving distribution of material in the system. The present work uses effective medium theory (EMT) for the electromagnetic field calculations since the goal is to apply the tool for industrial scale reacting systems. Predictions using MFIX/COMSOL and MFIX/Elmer are compared to in bed thermocouple measurements from a fluidized bed experiment and recent simulation predictions and recent experiments performed at NETL. The NETL experiments also include packed bed operating conditions. Coupled MFIX/COMSOL simulations at these conditions are compared with similar simulations using only COMSOL to measure errors due to the solver coupling approach.

CFD↗

A performance comparison of the Cray-2 and the Cray X-MP

A suite of thirteen large Fortran benchmark codes were run on Cray-2 and Cray X-MP supercomputers. These codes were a mix of compute-intensive scientific application programs (mostly Computational Fluid Dynamics) and some special vectorized computation exercise programs. For the general class of programs tested on the Cray-2, most of which were not specially tuned for speed, the floating point operation rates varied under a variety of system load configurations from 40 percent up to 125 percent of X-MP performance rates. It is concluded that the Cray-2, in the original system configuration studied (without memory pseudo-banking) will run untuned Fortran code, on average, about 70 percent of X-MP speeds.

Schmickley, Ronald↗

Comparison of methodologies for describing relaxation in nonequilibrium gaseous systems

The heat transfer process in hypervelocity vehicles is dominated by nonequilibrium gas dynamics. One model used in computational fluid dynamics (CFD) codes to predict hypervelocity heat transfer is the 'two-temperature' model. An analysis has been made to test the validity of the two-temperature model for predicting another nonequilibrium phenomenon, sound absorption and deviation of signal speed in a high temperature gas. It is found that the two temperature model's prediction capabilities degenerate with increasing temperature. These results are felt to have significance concerning the two-temperature's ability to predict heat transfer in hypervelocity flows.

Schreiber, Willard C.↗

Parallel CFD design on network-based computer

Combining multiple engineering workstations into a network-based heterogeneous parallel computer allows application of aerodynamic optimization with advanced computational fluid dynamics codes, which can be computationally expensive on mainframe supercomputers. This paper introduces a nonlinear quasi-Newton optimizer designed for this network-based heterogeneous parallel computing environment utilizing a software called Parallel Virtual Machine. This paper will introduce the methodology behind coupling a Parabolized Navier-Stokes flow solver to the nonlinear optimizer. This parallel optimization package is applied to reduce the wave drag of a body of revolution and a wing/body configuration with results of 5% to 6% drag reduction.

Cheung, Samson↗

Numerical Stability and Control Analysis Towards Falling-Leaf Prediction Capabilities of Splitflow for Two Generic High-Performance Aircraft Models

Aerodynamic analysis are performed using the Lockheed-Martin Tactical Aircraft Systems (LMTAS) Splitflow computational fluid dynamics code to investigate the computational prediction capabilities for vortex-dominated flow fields of two different tailless aircraft models at large angles of attack and sideslip. These computations are performed with the goal of providing useful stability and control data to designers of high performance aircraft. Appropriate metrics for accuracy, time, and ease of use are determined in consultations with both the LMTAS Advanced Design and Stability and Control groups. Results are obtained and compared to wind-tunnel data for all six components of forces and moments. Moment data is combined to form a "falling leaf" stability analysis. Finally, a handful of viscous simulations were also performed to further investigate nonlinearities and possible viscous effects in the differences between the accumulated inviscid computational and experimental data.

Charlton, Eric F.↗

Assessment of Various Flow Solvers Used to Predict the Thermal Environment inside Space Shuttle Solid Rocket Motor Joints

It is very important to accurately predict the gas pressure, gas and solid temperature, as well as the amount of O-ring erosion inside the space shuttle Reusable Solid Rocket Motor (RSRM) joints in the event of a leak path. The scenarios considered are typically hot combustion gas rapid pressurization events of small volumes through narrow and restricted flow paths. The ideal method for this prediction is a transient three-dimensional computational fluid dynamics (CFD) simulation with a computational domain including both combustion gas and surrounding solid regions. However, this has not yet been demonstrated to be economical for this application due to the enormous amount of CPU time and memory resulting from the relatively long fill time as well as the large pressure and temperature rising rate. Consequently, all CFD applications in RSRM joints so far are steady-state simulations with solid regions being excluded from the computational domain by assuming either a constant wall temperature or no heat transfer between the hot combustion gas and cool solid walls.

Wang, Qun-Zhen↗

Rotor Airloads Prediction Using Loose Aerodynamic Structural Coupling

This work couples a computational fluid dynamics (CFD) code and rotorcraft computational structural dynamics (CSD) code to calculate helicopter rotor airloads across a range of flight conditions. An iterative loose (weak) coupling methodology is used to couple the CFD and CSD codes on a per revolution, periodic basis. The CFD uses a high fidelity, Navier-Stokes, overset grid methodology with first principles-based wake capturing. Modifications are made to the CFD code for aeroelastic analysis. For a UH-60A Blackhawk helicopter, four challenging level flight conditions are computed: 1) low speed (u = 0.15) with blade-vortex interaction, 2) high speed (u = 0.37) with advancing blade negative lift, 3) high thrust with dynamic stall (u = 0.24), and 4) hover. Results are compared with UH-60A Airloads Program fight test data. Most importantly, for all cases the loose coupling methodology is shown to be stable, convergent, and robust with full coupling of normal force, pitching moment, and chord force. In comparison with flight test data, normal force and pitching moment magnitudes are in good agreement. For the high speed and dynamic stall cases a phase lag in comparison with the data is seen, nonetheless, the shapes of the curves are very good. Overall, the results are noteworthy improvement over lifting line aerodynamics used in rotorcraft comprehensive codes.

Potsdam, Mark↗

Theoretical and Experimental Unsteady Aerodynamics Compared for a Linear Oscillating Cascade With a Supersonic Leading-Edge Locus

Experimental data were obtained to help validate analytical and computational fluid dynamics (CFD) codes used to compute unsteady cascade aerodynamics in a supersonicaxial- flow regime. Results from two analytical codes and one CFD code were compared with experimental data. One analytical code did not account for airfoil thickness or camber; another, using piston theory (piston code), accounted for thickness and camber upstream of the first shockwave/airfoil impingement locations. The Euler CFD code accounted fully for airfoil shape.

Ramsey, John K.↗

Computational Aeroelastic Analysis of the Ares Launch Vehicle During Ascent

This paper presents the static and dynamic computational aeroelastic (CAE) analyses of the Ares crew launch vehicle (CLV) during atmospheric ascent. The influence of launch vehicle flexibility on the static aerodynamic loading and integrated aerodynamic force and moment coefficients is discussed. The ultimate purpose of this analysis is to assess the aeroelastic stability of the launch vehicle along the ascent trajectory. A comparison of analysis results for several versions of the Ares CLV will be made. Flexible static and dynamic analyses based on rigid computational fluid dynamic (CFD) data are compared with a fully coupled aeroelastic time marching CFD analysis of the launch vehicle.

Bartels, Robert E.↗

An Assessment of Comprehensive Code Prediction State-of-the-Art Using the HART II International Workshop Data

Despite significant advancements in computational fluid dynamics and their coupling with computational structural dynamics (= CSD, or comprehensive codes) for rotorcraft applications, CSD codes with their engineering level of modeling the rotor blade dynamics, the unsteady sectional aerodynamics and the vortical wake are still the workhorse for the majority of applications. This is especially true when a large number of parameter variations is to be performed and their impact on performance, structural loads, vibration and noise is to be judged in an approximate yet reliable and as accurate as possible manner. In this paper, the capabilities of such codes are evaluated using the HART II Inter- national Workshop data base, focusing on a typical descent operating condition which includes strong blade-vortex interactions. Three cases are of interest: the baseline case and two cases with 3/rev higher harmonic blade root pitch control (HHC) with different control phases employed. One setting is for minimum blade-vortex interaction noise radiation and the other one for minimum vibration generation. The challenge is to correctly predict the wake physics - especially for the cases with HHC - and all the dynamics, aerodynamics, modifications of the wake structure and the aero-acoustics coming with it. It is observed that the comprehensive codes used today have a surprisingly good predictive capability when they appropriately account for all of the physics involved. The minimum requirements to obtain these results are outlined.

vanderWall, Berend G.↗

Aerodynamic Loading on Crew Access Arm

An aerodynamic database for the crew access arm during deployment is generated. While initial attempts used data derived from wind tunnel testing, limitations lead to the use of computational fluid dynamics to improve analysis. This computational data compared favorably to the experimental data increasing confidence in both. Further analysis lead to simplification showing the ability to build the database with a single relative wind variable which takes into account both arm deployment angle and wind direction. Eventually geometric simplifications and conservative results lead to the use of simulations of the CAA alone to build the database.Finally these static conditions are then compared to a dynamic simulation where the CAA rotates into position to test the quasi-steady state assumptions used to generate the database. The dynamic case showed an increase in the moment of interest but within acceptable bounds for now.

SLS↗

Ground Wind Loads on the Space Launch System’s Mobile Launcher Crew Access Arm

An aerodynamic database for the crew access arm during deployment is generated. While initial attempts used data derived from wind tunnel testing, limitations lead to the use of computational fluid dynamics to improve analysis. This computational data compared favorably to the experimental data increasing confidence in both. Further analysis lead to simplification showing the ability to build the database with a single relative wind variable which takes into account both arm deployment angle and wind direction. Eventually geometric simplifications and conservative results lead to the use of simulations of the CAA alone to build the database.Finally these static conditions are then compared to a dynamic simulation where the CAA rotates into position to test the quasi-steady state assumptions used to generate the database. The dynamic case showed an increase in the moment of interest but within acceptable bounds for now.

SLS↗

High-Fidelity Analysis of Lift+Cruise VTOL Urban Air Mobility Concept Aircraft

This paper presents a high-fidelity multidisciplinary analysis for the NASA lift+cruise vertical takeoff and landing urban air mobility concept aircraft. The reported simulations couple comprehensive rotorcraft aeromechanics and high-fidelity computational fluid dynamics (CFD). Aerodynamic solutions are computed on dynamic, deforming, unstructured, overset grid systems by an unsteady CFD solver, FUN3D, developed at the NASA Langley Research Center. An integrated overset-grid assembler, Yoga, is used for communications between component grids. Two turbulence models are compared for loose-coupling simulations in hover and low-speed forward flight conditions. One model is the negative variant of the one-equation Spalart-Allmaras model (SA-neg), and the other model is the SA-neg-R model adding a simple rotation correction term to the SA-neg model. It demonstrated that the rotation correction in the SA-neg-R turbulence model significantly improves resolution of secondary vortices and wake interactions for rotorcraft simulations.

CFD↗

High-Fidelity Analysis of Lift+Cruise VTOL Urban Air Mobility Concept Aircraft

This paper presents a high-fidelity multidisciplinary analysis for the NASA Lift+Cruise vertical takeoff and landing urban air mobility concept aircraft. The reported simulations couple comprehensive rotorcraft aeromechanics and high-fidelity computational fluid dynamics (CFD). Aerodynamic solutions are computed on dynamic, deforming, unstructured, overset grid systems by an unsteady CFD solver, FUN3D, developed at the NASA Langley Research Center. An integrated overset-grid assembler, Yoga, is used for communications between component grids. Two turbulence models are compared for loose-coupling simulations in hover and low-speed forward flight conditions. One model is the negative variant of the one-equation Spalart-Allmaras model (SA-neg), and the other model is the SA-neg-R model adding a simple rotation correction term to the SA-neg model. It demonstrated that the rotation correction in the SA-neg-R turbulence model significantly improves resolution of secondary vortices and wake interactions for rotorcraft simulations.

CFD↗

Two-Way Coupled CFD/DEM Simulation of Electromagnetic Heating in a Fluidized Bed Reactor - Extended Abstract

A set of tools has been developed to couple a multiphase reacting computational fluid dynamics solver (MFIX) with two computational electromagnetic solvers (COMSOL and Elmer). In contrast to several recent investigations, which calculate a single static EM field prior to the CFD calculation, in the present work the electromagnetic field is periodically updated during the simulation. This enables the EM field and associated EM heating predictions to reflect the effects of evolving material properties (temperature, composition, size etc.) and the evolving distribution of material in the system. The present work uses effective medium theory (EMT) for the electromagnetic field calculations since the goal is to apply the tool for industrial scale reacting systems. Predictions using MFIX/COMSOL and MFIX/Elmer are compared to in bed thermocouple measurements from a fluidized bed experiment and recent simulation predictions and recent experiments performed at NETL. The NETL experiments also include packed bed operating conditions. Coupled MFIX/COMSOL simulations at these conditions are compared with similar simulations using only COMSOL to measure errors due to the solver coupling approach.

CFD↗

Transition Modeling Based on the Dual N-factor Method for the CRM-NLF Wind Tunnel Configuration

The dual N-factor method is used to model the boundary-layer transition over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. The paper presents a systematic methodology for transition prediction in the presence of a dual shock system and extends the prior capability for iteratively coupled computational fluid dynamics (CFD) predictions to incorporate three-dimensional, transonic wings. The method employs stability computations based on the linear parabolized stability equations (PSE), along with a dual N-factor criterion. The iterative process begins with the fully turbulent Reynolds-averaged-Navier-Stokes (RANS) mean flow solution. For the first iteration, a mean flow solution is calculated with an imposed transition front that aligns with the shock front from the fully turbulent solution. Subsequently, stability computations are performed along a set of streamlines across the wing to calculate the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) modes. The transition criterion based on the dual N-factor method is used to infer the updated transition front and the process is successively repeated until convergence of the solution. Within three iterations, the predicted fronts for angles of attack of 1.45, 1.98, 2.46 and 2.94 degrees and a mean-aerodynamic-chord Reynolds number equal to 15 million, approach visual convergence in most regions of the studied cases, and the resulting predictions are in good agreement with the transition fronts deduced from measurements of temperature-sensitive paint. Even though surface pressure measurements based on fully-turbulent flow agree well with the measured pressure coefficient distributions, strong viscous-inviscid interaction effects cause significant shifts in the shock locations based on the imposed transition front, underscoring the intrusive nature of static pressure measurements using surface mounted ports on the CRM-NLF configuration.

Boundary Layer Transition↗