Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

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 1,153 records · Page 64

A Framework for Integrated Component and System Analyses of Instabilities

Instabilities associated with fluid handling and operation in liquid rocket propulsion systems and test facilities usually manifest themselves as structural vibrations or some form of structural damage. While the source of the instability is directly related to the performance of a component such as a turbopump, valve or a flow control element, the associated pressure fluctuations as they propagate through the system have the potential to amplify and resonate with natural modes of the structural elements and components of the system. In this paper, the authors have developed an innovative multi-level approach that involves analysis at the component and systems level. The primary source of the unsteadiness is modeled with a high-fidelity hybrid RANS/LES based CFD methodology that has been previously used to study instabilities in feed systems. This high fidelity approach is used to quantify the instability and understand the physics associated with the instability. System response to the driving instability is determined through a transfer matrix approach wherein the incoming and outgoing pressure and velocity fluctuations are related through a transfer (or transmission) matrix. The coefficients of the transfer matrix for each component (i.e. valve, pipe, orifice etc.) are individually derived from the flow physics associated with the component. A demonstration case representing a test loop/test facility comprised of a network of elements is constructed with the transfer matrix approach and the amplification of modes analyzed as the instability propagates through the test loop.

Ahuja, Vineet↗

Efficiency Considerations in Low Pressure Turbines

Issues & Topics Discussed: a) Aviation Week reported shortfall In LPT efficiency due to the application of "high lift airfoils". b) Progress in the design technologies in LPTs during the last 20 years: 1) Application of RANS based CFD codes. 2) Integration of recent experimental data and modeling of LPT airfoil specific flows into design methods. c) Opportunities to further enhance LPT efficiency for commercial aviation and military transport application and to impact emissions, noise, weight & cost.

Source record↗

Reynolds-Averaged Navier-Stokes Simulation of a 2D Circulation Control Wind Tunnel Experiment

Numerical simulations are performed using a Reynolds-averaged Navier-Stokes (RANS) flow solver for a circulation control airfoil. 2D and 3D simulation results are compared to a circulation control wind tunnel test conducted at the NASA Langley Basic Aerodynamics Research Tunnel (BART). The RANS simulations are compared to a low blowing case with a jet momentum coefficient, C(sub u), of 0:047 and a higher blowing case of 0.115. Three dimensional simulations of the model and tunnel walls show wall effects on the lift and airfoil surface pressures. These wall effects include a 4% decrease of the midspan sectional lift for the C(sub u) 0.115 blowing condition. Simulations comparing the performance of the Spalart Allmaras (SA) and Shear Stress Transport (SST) turbulence models are also made, showing the SST model compares best to the experimental data. A Rotational/Curvature Correction (RCC) to the turbulence model is also evaluated demonstrating an improvement in the CFD predictions.

Allan, Brian G.↗

Time Accurate CFD Simulations of the Orion Launch Abort Vehicle in the Transonic Regime

Significant asymmetries in the fluid dynamics were calculated for some cases in the CFD simulations of the Orion Launch Abort Vehicle through its abort trajectories. The CFD simulations were performed steady state and in three dimensions with symmetric geometries, no freestream sideslip angle, and motors firing. The trajectory points at issue were in the transonic regime, at 0 and +/- 5 angles of attack with the Abort Motors with and without the Attitude Control Motors (ACM) firing. In some of the cases the asymmetric fluid dynamics resulted in aerodynamic side forces that were large enough that would overcome the control authority of the ACMs. MSFC's Fluid Dynamics Group supported the investigation into the cause of the flow asymmetries with time accurate CFD simulations, utilizing a hybrid RANS-LES turbulence model. The results show that the flow over the vehicle and the subsequent interaction with the AB and ACM motor plumes were unsteady. The resulting instantaneous aerodynamic forces were oscillatory with fairly large magnitudes. Time averaged aerodynamic forces were essentially symmetric.

Rojahn, Josh↗

Time Accurate CFD Simulations of the Orion Launch Abort Vehicle in the Transonic Regime

Significant asymmetries in the fluid dynamics were calculated for some cases in the CFD simulations of the Orion Launch Abort Vehicle through its abort trajectories. The CFD simulations were performed steady state with symmetric boundary conditions and geometries. The trajectory points at issue were in the transonic regime, at 0 and 5 angles of attack with the Abort Motors with and without the Attitude Control Motors (ACM) firing. In some of the cases the asymmetric fluid dynamics resulted in aerodynamic side forces that were large enough that would overcome the control authority of the ACMs. MSFC s Fluid Dynamics Group supported the investigation into the cause of the flow asymmetries with time accurate CFD simulations, utilizing a hybrid RANS-LES turbulence model. The results show that the flow over the vehicle and the subsequent interaction with the AB and ACM motor plumes were unsteady. The resulting instantaneous aerodynamic forces were oscillatory with fairly large magnitudes. Time averaged aerodynamic forces were essentially symmetric.

Ruf, Joseph↗

Uncertainty Quantification of the FUN3D-Predicted NASA CRM Flutter Boundary

A nonintrusive point collocation method is used to propagate parametric uncertainties of the flexible Common Research Model, a generic transport configuration, through the unsteady aeroelastic CFD solver FUN3D. A range of random input variables are considered, including atmospheric flow variables, structural variables, and inertial (lumped mass) variables. UQ results are explored for a range of output metrics (with a focus on dynamic flutter stability), for both subsonic and transonic Mach numbers, for two different CFD mesh refinements. A particular focus is placed on computing failure probabilities: the probability that the wing will flutter within the flight envelope.

Stanford, Bret K.↗

Preliminary Computational Assessment of Disk Rotating Detonation Engine Configurations

The pressure gain combustion (PGC) community is currently investigating rotating detonation engine (RDE) configurations where the flow direction is predominantly radial while the detonation travels circumferentially. These configurations are sometimes referred to as disk rotating detonation engines (DRDE) due to their nominal appearance as two disks in parallel with a gap between them. Having radial flow between disks, as opposed to the conventional RDE with axial flow in an annulus, may have profound effects on both the flow field and the performance. It may also yield extraordinarily compact devices which are well suited to particular propulsion and power applications. This presentation describes a preliminary effort to model the DRDE using a modified computational fluid dynamics (CFD) code originally written for analyzing ordinary RDE's. The quasi-two-dimensional code modifications are described, and some simple test flows are analyzed to insure that the modifications are functioning as envisioned. The code is then used to examine several DRDE scenarios such as radially inward and radially outward devices to see if stable operation is possible and if so, to assess the performance in terms of pressure gain. It is found that several flow scenarios are not only stable, but show superior performance to the ordinary RDE.

Paxson, Daniel E.↗

An Experimental Investigation of a Wing-Fuselage Junction Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel

Current turbulence models, such as those employed in Reynolds-averaged Navier-Stokes CFD, are unable to reliably predict the onset and extent of the three-dimensional separated flow that typically occurs in wing-fuselage junctions. To critically assess, as well as to improve upon, existing turbulence models, experimental validation-quality flow-field data in the junction region is needed. In this report, we present an overview of experimental measurements on a wing-fuselage junction model that addresses this need. The experimental measurements were performed in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The model was a full-span wing-fuselage body that was configured with truncated DLR-F6 wings, both with and without leading-edge extensions at the wing root. The model was tested at a fixed chord Reynolds number of 2.4 million, and angles-of-attack ranging from -10 degrees to +10 degrees were considered. Flow-field measurements were performed with a pair of miniature laser Doppler velocimetry (LDV) probes that were housed inside the model and attached to three-axis traverse systems. One LDV probe was used to measure the separated flow field in the trailing-edge junction region. The other LDV probe was alternately used to measure the flow field in the leading-edge region of the wing and to measure the incoming fuselage boundary layer well upstream of the leading edge. Both LDV probes provided measurements from which all three mean velocity components, all six independent components of the Reynolds-stress tensor, and all ten independent components of the velocity triple products were calculated. In addition to the flow-field measurements, static and dynamic pressures were measured at selected locations on the wings and fuselage of the model, infrared imaging was used to characterize boundary-layer transition, oil-flow visualization was used to visualize the separated flow in the leading- and trailing-edge regions of the wing, and unsteady shear stress was measured at limited locations using capacitive shear-stress sensors. Sample results from the measurement techniques employed during the test are presented and discussed.

Kegerise, Michael A.↗

March of the Interns

Summer 2019 brought 70 interns — from Puerto Rico, 44 continental US states, Sweden, and the Republic of Trinidad and Tobago — into the NASA Ames Aeromechanics Branch. Some were high school students learning engineering for the first time, while most were mechanical and aerospace engineering students, though two physics majors and a math major jumped into the mix to provide balance. This year’s interns completed work on 30 different projects with a collective of 30,000 dedicated hours of work. The projects this year were on urban air mobility (UAM), search and rescue vehicle design, all things Mars, Titan exploration concept designs, 3D modeling (i.e. computer aided design, CAD), computational fluid dynamics (CFD), and much more. Last year we addressed the question of what to do with an army of interns. This year, we are following its path to change the world.

Dromiack, Hannah↗

TPSAS-NF1676L-31264-DND

Nonlinear unsteady aerodynamic models (TSD, Euler/Navier-Stokes (RANS), LES) are complex, expensive solutions. CFD used by industry mostly for steady/static loads cases. Limit Cycle Oscillation (LCO) is an equilibrium point of the nonlinear aeroelastic system induced by nonlinear unsteady aerodynamics (shock dynamics, separated flow) and/or nonlinear structural dynamics (geometric and material nonlinearities). While Euler solutions can generate LCOs, viscosity yields improved accuracy (shock strength, position). Most common aeroelastic CFD codes currently consist of nonlinear unsteady aerodynamics with linear structural dynamics. Van der Pol system is a simple example of LCOs.

Walter A Silva↗

Improved Prediction of Losses with Large Eddy Simulation in a Low-Pressure Turbine

There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of freestream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Reex=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.

CFD↗

Flow Reconstruction in A Transonic Turbine Cascade Using Physics-Informed Neural Networks (PINNS)

This paper investigates the application of Physics-Informed Neural Networks (PINNs) for the analysis of turbine blades in a transonic cascade. The 2-D flow field in a transonic turbine cascade is reconstructed in three ways: the traditional forward approach (PINN not trained on experimental data), by training the PINN using discrete sets of experimentally measured pressure at midspan, and in the inverse sense where no inlet or outlet pressure boundary conditions are applied. Comparisons between the PINN solutions to measured data are made. This is repeated for three different turbine blades with distinct loading characteristics. Good agreement is shown between a CFD calculation of the CMC7 blade, and the PINN model trained with all data. The PINN is trained utilizing all available data, half the available data, data from only the leading edge region, and data from only the trailing edge region. The forward problem results deviate the most from experimental data but show promise. Solutions from the assisted training cases show that the PINN can reconstruct the flow field with acceptable accuracy when trained on measurements along the entire blade. In the inverse case, it is shown that to simultaneously achieve acceptable errors for inlet Mach number and outlet isentropic Mach number, the PINN must be trained on the static pressure data along the entire blade.

Machine Learning↗

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

The validation and application of numerical modelling to supersonic mixing and reacting flows

This study demonstrates that the computational fluid dynamics (CFD) code GASP can be used to predict the features for flowfields representative of scramjet combustors. The validation program undertaken was composed of two parts: a low enthalpy mixing study and a simplified autoignition study. The low enthalpy mixing study was done using low angle helium injection into a Mach 6 airstream over a flat plate. Both matched and overpressurized injection cases were numerically modelled. The agreement with experimental data for the farfield bulk mixing in both cases was reasonable. The combustion studies modelled premixed hydrogen-air reaction with the results being evaluated qualitatively. GASP proved to be highly sensitive to the chemical species present for flow near the ignition temperature, but provided reasonable qualitative results for a subsonic flameholder model.

Vitt, Paul H.↗

CFD Sensitivity Analysis of a Modern Civil Transport Near Buffet-Onset Conditions

A computational fluid dynamics (CFD) sensitivity analysis is conducted for a modern civil transport at several conditions ranging from mostly attached flow to flow with substantial separation. Two different Navier-Stokes computer codes and four different turbulence models are utilized, and results are compared both to wind tunnel data at flight Reynolds number and flight data. In-depth CFD sensitivities to grid, code, spatial differencing method, aeroelastic shape, and turbulence model are described for conditions near buffet onset (a condition at which significant separation exists). In summary, given a grid of sufficient density for a given aeroelastic wing shape, the combined approximate error band in CFD at conditions near buffet onset due to code, spatial differencing method, and turbulence model is: 6% in lift, 7% in drag, and 16% in moment. The biggest two contributers to this uncertainty are turbulence model and code. Computed results agree well with wind tunnel surface pressure measurements both for an overspeed 'cruise' case as well as a case with small trailing edge separation. At and beyond buffet onset, computed results agree well over the inner half of the wing, but shock location is predicted too far aft at some of the outboard stations. Lift, drag, and moment curves are predicted in good agreement with experimental results from the wind tunnel.

Rumsey, Christopher L.↗

Gas Phase Effects on Slosh Dynamics

Gas phase effects on slosh dynamics were quantified using computational fluid dynamics (CFD) simulation for a range of propellant and ullage gas combinations. Historical slosh modeling using potential flow solutions typically neglects gas phase effects. Regardless, the results have been shown to compare well with slosh ground tests typically performed with water and air at standard temperature and pressure. CFD analysis reveals that as the liquid-to-gas density ratio decreases, slosh dynamics change due to the relative increase in gas inertia and thus influence on liquid motion. The result is a profound impact on slosh dynamics over certain parameter spaces particularly for liquid hydrogen. Gas phase effects on slosh dynamics should be considered in slosh models especially for liquid hydrogen propellant tanks.

Jacob M Brodnick↗

Three-Dimensional Effects on Multi-Element High Lift Computations

In an effort to discover the causes for disagreement between previous 2-D computations and nominally 2-D experiment for flow over the 3-clement McDonnell Douglas 30P-30N airfoil configuration at high lift, a combined experimental/CFD investigation is described. The experiment explores several different side-wall boundary layer control venting patterns, document's venting mass flow rates, and looks at corner surface flow patterns. The experimental angle of attack at maximum lift is found to be sensitive to the side wall venting pattern: a particular pattern increases the angle of attack at maximum lift by at least 2 deg. A significant amount of spanwise pressure variation is present at angles of attack near maximum lift. A CFD study using 3-D structured-grid computations, which includes the modeling of side-wall venting, is employed to investigate 3-D effects of the flow. Side-wall suction strength is found to affect the angle at which maximum lift is predicted. Maximum lift in the CFD is shown to be limited by the growth of all off-body corner flow vortex and consequent increase in spanwise pressure variation and decrease in circulation. The 3-D computations with and without wall venting predict similar trends to experiment at low angles of attack, but either stall too earl or else overpredict lift levels near maximum lift by as much as 5%. Unstructured-grid computations demonstrate that mounting brackets lower die the levels near maximum lift conditions.

Rumsey, Christopher L.↗

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.↗