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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Europa Clipper Preliminary Design Review Propellant Slosh Analysis

Propellant slosh was analyzed for both the oxidizer and the fuel for the Europa Clipper propulsion system. Slosh was examined for various fill fractions for cases where acceleration was on the order of magnitude of 10(exp -2) m/sq. s using the computational fluid dynamics software package STAR-CCM+ and at various fill fractions for cases where acceleration was on the order of magnitude of 10(exp -5) m/sq. s using Surface Evolver. Equivalent mechanical model parameters were derived from the CFD data using MATLAB for both the higher and the lower acceleration slosh cases. These parameters were plotted and can be used to interpolate mechanical model parameters at fill fractions not analyzed by CFD or Surface Evolver.

CFD↗
NASA NTRS2018

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗
NASA NTRS1998

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗
NASA NTRS2023

Turbulence Model Validation Through Joint Experimental /Computational Studies of Separated Flow Over A Three-Dimensional Tapered Bump: Part I - Experimental Investigation

The results of a three-year experimental campaign aimed at comprehensively documenting the separated flow over a three-dimensional bump are presented with the purpose of generating a benchmark experimental database useful in validating computational fluid dynamics (CFD) flow simulations. The bump model geometry was designed to provide well-defined and repeatable smooth-body flow separation conditions that were suitable for both experiments and simulations. The bump had a Gaussian streamwise profile with a constant height equal to 8.5% of its width over the central 60% of the test section width. The remaining 40% were outboard spanwise portions that gradually taper to zero using an error function profile to minimize side-wall boundary layer effects. The model was immersed in a turbulent boundary layer that was developed on a suspended flat plate in the Notre Dame Mach 0.6 Wind Tunnel. In order to document the effect of the incoming boundary layer thickness on the flow separation, the bump model could be located at two streamwise positions. The mean velocity and turbulence intensity of the wind tunnel freestream flow field and approaching turbulent boundary layer were fully documented. The measurements of the flow separation region included surface visualization, wall shear stress using oil-film interferometry, mean and dynamic surface pressure, and planar and stereoscopic particle image velocimetry. The experiments were conducted over a range of Mach numbers from 0.05 to 0.2 corresponding to a range of Reynolds numbers based on the test section spanwise dimension (0.914m) of 1.0 × 10 6 ≤ Re L = U ∞ L/ν ≤ 4.0 × 10 6 . The bulk of the results are presented for the higher Mach number conditions of 0.1 and 0.2 with Re L = 2.0 × 10 6 and 4.0 × 10 6 , respectively. Extensive uncertainty analysis of the data was performed. The data is archived in the NASA Langley Turbulence Modeling Resource website at https://turbmodels.larc.nasa.gov/Other_exp_Data/speedbump_sep_exp.html. In addition to the experiments, a computational effort was made in parallel by the CFD group at Boeing Research & Technology highlighting the usefulness of the data set, which is outlined in the accompanying CFD report.

Flow separation↗
NASA NTRSDate not supplied

Low-cost Quantification of Fluid Flow Parameter Sensitivity using Reduced-order Modeling

Uncertainty quantification of computational fluid dynamics (CFD) simulations is a complicated procedure which still relies in many cases on engineering judgment and factors of safety. This is in part because the computational cost of measuring the simulation's sensitivity to all meaningful parameters (e.g., body surface roughness) and hyperparameters (e.g., subiteration convergence criterion) is intractable for even a single simulation. Reduced-order modeling dramatically lowers this computational cost of simulating fluid flows, but usually only where similar data is already available. In this work, fluid reduced-order models are utilized to quantify a flow's sensitivity to certain physical parameters for the purposes of improved uncertainty quantification. Characteristic observability and sensitivity are both explored. The resulting sensitivity quantifications enable more informed CFD frameworks and more rigorous uncertainty bounds on the resulting data.

uncertainty quantification↗
NASA NTRSDate not supplied

Low-cost Quantification of Fluid Flow Parameter Sensitivity using Reduced-order Modeling

Uncertainty quantification of computational fluid dynamics (CFD) simulations is a complicated procedure which still relies in many cases on engineering judgment and factors of safety. This is in part because the computational cost of measuring the simulation's sensitivity to all meaningful parameters (e.g., body surface roughness) and hyperparameters (e.g., subiteration convergence criterion) is intractable for even a single simulation. Reduced-order modeling dramatically lowers this computational cost of simulating fluid flows, but usually only where similar data is already available. In this work, fluid reduced-order models are utilized to quantify a flow's sensitivity to certain physical parameters for the purposes of improved uncertainty quantification. Characteristic observability and sensitivity are both explored. The resulting sensitivity quantifications enable more informed CFD frameworks and more rigorous uncertainty bounds on the resulting data.

uncertainty quantification↗
NASA NTRS1997

Enhanced Control Effector Designs for Airbreathing Transatmospheric Vehicles

A study was conducted to evaluate the potential effectiveness of a moveable cowl-trailing-edge design for airbreathing hypersonic single-stage-to-orbit (SSTO) configurations, which can be extended or deflected from the nominal position in order to provide additional pitching moment capability. This additional pitching moment capability may reduce the necessary deflection angle from conventional control surfaces and the associated trim drag penalty. Calculations for a generic SSTO configuration with baseline and modified cowl trailing edge geometries were performed at Mach 6 and 10 initially with an engineering analysis code in order to examine several design parametrics. In order to more accurately model geometries and flow physics, 2-D viscous computational fluid dynamics (CFD) predictions were obtained. FInally, a limited set of 3D CFD predictions were obtained at Mach 6 in order to show the effects of modeling 3D flow fields as well as the full 3D vehicle geometry. Comparisons of aftbody surface pressures and force and moment predictions show differences between initial predictions and 2-D CFD solutions due to geometry modeling and calculation method differences. The 3-D CFD predictions confirm the trends observed in the 2-D solutions and provide additional information on 3D effects. These analyses shows that cowl-trailing-edge extensions were effective in providing additional (nose-down) pitching moment increments as well as increased thrust compared to the baseline geometry. These effects reduce the control surface deflection angle required to trim the vehicle and the associated trim drag. Cowl-trailingedge deflections were not as effective, generating more nose-up pitching moment and decreased thrust compared to non-deflected cases.

Cockrell, Charles E.↗
NASA NTRS2004

Evaluation of CFD to Determine Two-Dimensional Airfoil Characteristics for Rotorcraft Applications

The efficient prediction of helicopter rotor performance, vibratory loads, and aeroelastic properties still relies heavily on the use of comprehensive analysis codes by the rotorcraft industry. These comprehensive codes utilize look-up tables to provide two-dimensional aerodynamic characteristics. Typically these tables are comprised of a combination of wind tunnel data, empirical data and numerical analyses. The potential to rely more heavily on numerical computations based on Computational Fluid Dynamics (CFD) simulations has become more of a reality with the advent of faster computers and more sophisticated physical models. The ability of five different CFD codes applied independently to predict the lift, drag and pitching moments of rotor airfoils is examined for the SC1095 airfoil, which is utilized in the UH-60A main rotor. Extensive comparisons with the results of ten wind tunnel tests are performed. These CFD computations are found to be as good as experimental data in predicting many of the aerodynamic performance characteristics. Four turbulence models were examined (Baldwin-Lomax, Spalart-Allmaras, Menter SST, and k-omega).

Smith, Marilyn J.↗
NASA NTRS1996

LES on unstructured deforming meshes: Towards reciprocating IC engines

A variable explicit/implicit characteristics-based advection scheme that is second-order accurate in space and time has been developed recently for unstructured deforming meshes (O'Rourke & Sahota 1996a). To explore the suitability of this methodology for Large-Eddy Simulation (LES), three subgrid-scale turbulence models have been implemented in the CHAD CFD code (O'Rourke & Sahota 1996b): a constant-coefficient Smagorinsky model, a dynamic Smagorinsky model for flows having one or more directions of statistical homogeneity, and a Lagrangian dynamic Smagorinsky model for flows having no spatial or temporal homogeneity (Meneveau et al. 1996). Computations have been made for three canonical flows, progressing towards the intended application of in-cylinder flow in a reciprocating engine. Grid sizes were selected to be comparable to the coarsest meshes used in earlier spectral LES studies. Quantitative results are reported for decaying homogeneous isotropic turbulence, and for a planar channel flow. Computations are compared to experimental measurements, to Direct-Numerical Simulation (DNS) data, and to Rapid-Distortion Theory (RDT) where appropriate. Generally satisfactory evolution of first and second moments is found on these coarse meshes; deviations are attributed to insufficient mesh resolution. Issues include mesh resolution and computational requirements for a specified level of accuracy, analytic characterization of the filtering implied by the numerical method, wall treatment, and inflow boundary conditions. To resolve these issues, finer-mesh simulations and computations of a simplified axisymmetric reciprocating piston-cylinder assembly are in progress.

Haworth, D. C.↗
NASA NTRS1990

CFD propels NASP propulsion progress

The most complex aerothermodynamics encountered in the National Aerospace Plane (NASP) propulsion system are associated with the fuel-mixing and combustion-reaction flows of its combustor section; adequate CFD tools must be developed to model shock-wave systems, turbulent hydrogen/air mixing, flow separation, and combustion. Improvements to existing CFD codes have involved extension from two dimensions to three, as well as the addition of finite-rate hydrogen-air chemistry. A novel CFD code for the treatment of reacting flows throughout the NASP, designated GASP, uses the most advanced upwind-differencing technology.

Povinelli, Louis A.↗
NASA NTRS2003

Internal Flow Thermal/Fluid Modeling of STS-107 Port Wing in Support of the Columbia Accident Investigation Board

As part of the aero-thermodynamics team supporting the Columbia Accident Investigation Board (CAB), the Marshall Space Flight Center was asked to perform engineering analyses of internal flows in the port wing. The aero-thermodynamics team was split into internal flow and external flow teams with the support being divided between shorter timeframe engineering methods and more complex computational fluid dynamics. In order to gain a rough order of magnitude type of knowledge of the internal flow in the port wing for various breach locations and sizes (as theorized by the CAB to have caused the Columbia re-entry failure), a bulk venting model was required to input boundary flow rates and pressures to the computational fluid dynamics (CFD) analyses. This paper summarizes the modeling that was done by MSFC in Thermal Desktop. A venting model of the entire Orbiter was constructed in FloCAD based on Rockwell International s flight substantiation analyses and the STS-107 reentry trajectory. Chemical equilibrium air thermodynamic properties were generated for SINDA/FLUINT s fluid property routines from a code provided by Langley Research Center. In parallel, a simplified thermal mathematical model of the port wing, including the Thermal Protection System (TPS), was based on more detailed Shuttle re-entry modeling previously done by the Dryden Flight Research Center. Once the venting model was coupled with the thermal model of the wing structure with chemical equilibrium air properties, various breach scenarios were assessed in support of the aero-thermodynamics team. The construction of the coupled model and results are presented herein.

Sharp, John R.↗
NASA NTRS2021

The Causes of Propeller Pitching Moment and the Conditions for its Significance

Recent development of vertical takeoff and landing (VTOL) aircraft has renewed interest in the study of propellers. One metric in particular, the propeller pitching moment, has been observed to be important to VTOL aircraft stability and control in the past. Propellers at angles of attack could not be accurately modeled in generations past due to a lack of computational power, but even with advances in computer technology, modern designers seem to possess insufficient knowledge in this area. In this dissertation, we study the physics behind propeller pitching moment in the context of an isolated propeller and a propeller upstream of a wing. An unsteady 3D vortex lattice method is developed specifically to model propellers at angles of attack and is validated by comparing to high-fidelity CFD analyses. We then use the model to isolate velocity influences to show that the propeller pitching moment is largely caused by two effects: a skewed wake and the presence of wing circulation. Generated maps of propeller pitching moment over a range of operational parameters corresponding to VTOL transition show that the low flight speeds and high angles of attack encountered during transition lead to significant magnitudes of propeller pitching moment that would be difficult to trim using passive methods. Also, derivation of a generalizable metric of significance shows that the peak contribution of propeller pitching moment to aircraft stability is comparable to a longitudinal displacement of the center of gravity by several percent of the wing chord. Finally, we give a concluding discussion on the impact of propeller pitching moment on VTOL aircraft design.

Xiaofan Fei↗
NASA NTRS2004

The Acoustic Analogy: A Powerful Tool in Aeroacoustics with Emphasis on Jet Noise Prediction

The acoustic analogy introduced by Lighthill to study jet noise is now over 50 years old. In the present paper, Lighthill s Acoustic Analogy is revisited together with a brief evaluation of the state-of-the-art of the subject and an exploration of the possibility of further improvements in jet noise prediction from analytical methods, computational fluid dynamics (CFD) predictions, and measurement techniques. Experimental Particle Image Velocimetry (PIV) data is used both to evaluate turbulent statistics from Reynolds-averaged Navier-Stokes (RANS) CFD and to propose correlation models for the Lighthill stress tensor. The NASA Langley Jet3D code is used to study the effect of these models on jet noise prediction. From the analytical investigation, a retarded time correction is shown that improves, by approximately 8 dB, the over-prediction of aft-arc jet noise by Jet3D. In experimental investigation, the PIV data agree well with the CFD mean flow predictions, with room for improvement in Reynolds stress predictions. Initial modifications, suggested by the PIV data, to the form of the Jet3D correlation model showed no noticeable improvements in jet noise prediction.

Farassat, F.↗
NASA NTRSDate not supplied

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55-deg Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of testsin a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 55-deg sphere-cone models placed in a free jet downstream of the 13-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream.Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon gages and five pressure tabs.The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. These simulations take into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗
NASA NTRSDate not supplied

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55-deg Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 55-deg sphere-cone models placed in a free jet downstream of the 13-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream. Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon gages and five pressure tabs. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. These simulations take into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗
NASA NTRSDate not supplied

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55-deg Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 55-deg sphere-cone models placed in a free jet downstream of the 13-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream. Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon gages and five pressure tabs. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. These simulations take into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗
NASA NTRS2011

Two-Point Turbulence Closure Applied to Variable Resolution Modeling

Variable resolution methods have become frontline CFD tools, but in order to take full advantage of this promising new technology, more formal theoretical development is desirable. Two general classes of variable resolution methods can be identified: hybrid or zonal methods in which RANS and LES models are solved in different flow regions, and bridging or seamless models which interpolate smoothly between RANS and LES. This paper considers the formulation of bridging methods using methods of two-point closure theory. The fundamental problem is to derive a subgrid two-equation model. We compare and reconcile two different approaches to this goal: the Partially Integrated Transport Model, and the Partially Averaged Navier-Stokes method.

Girimaji, Sharath S.↗
NASA NTRS2002

Effect of Favorable Pressure Gradients on Turbine Blade Pressure Surface Heat Transfer

Recent measurements on a turbine rotor showed significant relaminarization effects. These effects were evident on the pressure surface heat transfer measurements. The character of the heat transfer varied with Reynolds number. Data were obtained for exit Reynolds numbers between 500,000 and 880,000. Tests were done with a high level of inlet turbulence, 7.5%. At lower Reynolds numbers the heat transfer was similar to that for laminar flow, but at a level higher than for laminar flow. At higher Reynolds numbers the heat transfer was similar to turbulent flow, when the acceleration parameter, K, was sufficiently small. The proposed paper discusses the experimental results, and also discusses approaches to calculating the surface heat transfer for the blade surface. Calculations were done using a three-dimensional Navier-Stokes CFD analysis. The results of these tests, when compared with previous blade tests in the same facility, illustrate modeling difficulties that were encountered in CFD predictions. The two blades were in many ways similar. However, the degree of agreement between the same analysis and the experimental data was significantly different. These differences are highlighted to illustrate where improvements in modeling approaches are needed for transitional flows.

Boyle, Robert J.↗