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,207 records · Page 67

Uncertainty Assessment of CFD Investigation of the Nonlinear Difference-Frequency Wave Loads on a Semisubmersible FOWT Platform

Current mid-fidelity modeling approaches for floating offshore wind turbines (FOWTs) have been found to underpredict the nonlinear, low-frequency wave excitation and the response of semisubmersible FOWTs. To examine the cause of this underprediction, the OC6 project is using computational fluid dynamics (CFD) tools to investigate the wave loads on the OC5-DeepCwind semisubmersible, with a focus on the nonlinear difference-frequency excitation. This paper focuses on assessing the uncertainty of the CFD predictions from simulations of the semisubmersible in a fixed condition under bichromatic wave loading and on establishing confidence in the results for use in improving mid-fidelity models. The uncertainty for the nonlinear wave excitation is found to be acceptable but larger than that for the wave-frequency excitation, with the spatial discretization error being the dominant contributor. Further, unwanted free waves at the difference frequency have been identified in the CFD solution. A wave-splitting and wave load-correction procedure are presented to remove the contamination from the free waves in the results. A preliminary comparison to second-order potential-flow theory shows that the CFD model predicted significantly higher difference-frequency wave excitations, especially in surge, suggesting that the CFD results can be used to better calibrate the mid-fidelity tools.

17 WIND ENERGY↗

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↗

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↗

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↗

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

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

Advanced turbulence models for turbomachinery

Development and assessment of the single-time-scale k-epsilon turbulence model with different near-wall treatments and the multi-scale turbulence model for rotating flows are presented. These turbulence models are coded as self contained module decks that can be interfaced with a number of CFD main flow solvers. For each model, a stand-alone module deck with its own formulation, discretization scheme, solver and boundary condition implementations is presented. These satellite decks will take as input (from a main flow solver) the velocity field, grid, boundary condition specifications and will deliver turbulent quantities as output. These modules were tested as separate entities and, although many logical and programming problems were overcome, only wider use and further testing can render the modules sufficiently 'fool proof'.

Hadid, Ali H.↗

An Incidence Loss Model for Wave Rotors with Axially Aligned Passages

A simple mathematical model is described to account for the losses incurred when the flow in the duct (port) of a wave rotor is not aligned with the passages. The model, specifically for wave rotors with axially aligned passages, describes a loss mechanism which is sensitive to incident flow angle and Mach number. Implementation of the model in a one-dimensional CFD based wave rotor simulation is presented. Comparisons with limited experimental results are consistent with the model. Sensitivity studies are presented which highlight the significance of the incidence loss relative to other loss mechanisms in the wave rotor.

Paxson, Daniel E.↗

Boundary-Layer Transition Prediction Through Loose Coupling of OVERFLOW and LASTRAC

Transition prediction based on linear stability theory is expected to more accurately reflect the causality of transition onset than phenomenological transition models based on RANS-like transport equations. To help achieve the CFD vision 2030 aim of building a CFD tool chain with automated prediction of boundary layer transition, a technique to loosely tie the NASA OVERFLOW CFD solver with the LASTRAC stability analysis tool is described. The coupled solver is then used to compute transition over selected over a flat plate in a freestream with sufficiently low levels of turbulence, NLF(1)-0416 airfoil, the 6:1 prolate spheroid at an angle of attack, and a NASA juncture flow model with symmetric wing configuration. The findings show that the loosely coupled approach can reliably predict the transition location accurately in scenarios that are dominated by a single transition mechanism involving Tollmien-Schlichting instabilities, crossflow instabilities, or separation bubble-induced transition, or include a mixture of selected mechanisms. The toolset presents here appears to be robust to the prescription of the initial transition location, and it can lead to a converged solution in four or five rounds of the mean flow calculation and stability analysis, with minimal input from the user.

boundary layer transition↗

New Flutter Analysis Technique for CFD-based Unsteady Aeroelasticity

This paper presents a flutter analysis technique for the transonic flight regime. The technique uses an iterative approach to determine the critical dynamic pressure for a given mach number. Unlike other CFD-based flutter analysis methods, each iteration solves for the critical dynamic pressure and uses this value in subsequent iterations until the value converges. This process reduces the iterations required to determine the critical dynamic pressure. To improve the accuracy of the analysis, the technique employs a known structural model, leaving only the aerodynamic model as the unknown. The aerodynamic model is estimated using unsteady aeroelastic CFD analysis combined with a parameter estimation routine. The technique executes as follows. The known structural model is represented as a finite element model. Modal analysis determines the frequencies and mode shapes for the structural model. At a given mach number and dynamic pressure, the unsteady CFD analysis is performed. The output time history of the surface pressure is converted to a nodal aerodynamic force vector. The forces are then normalized by the given dynamic pressure. A multi-input multi-output parameter estimation software, ERA, estimates the aerodynamic model through the use of time histories of nodal aerodynamic forces and structural deformations. The critical dynamic pressure is then calculated using the known structural model and the estimated aerodynamic model. This output is used as the dynamic pressure in subsequent iterations until the critical dynamic pressure is determined. This technique is demonstrated on the Aerostructures Test Wing-2 model at NASA's Dryden Flight Research Center.

Pak, Chan-gi↗

Verification, Validation, and Solution Quality in Computational Physics: CFD Methods Applied to Ice Sheet Physics

Procedures and methods for veri.cation of coding algebra and for validations of models and calculations used in the aerospace computational fluid dynamics (CFD) community would be ef.cacious if used by the glacier dynamics modeling community. This paper presents some of those methods, and how they might be applied to uncertainty management supporting code veri.cation and model validation for glacier dynamics. The similarities and differences between their use in CFD analysis and the proposed application of these methods to glacier modeling are discussed. After establishing sources of uncertainty and methods for code veri.cation, the paper looks at a representative sampling of veri.cation and validation efforts that are underway in the glacier modeling community, and establishes a context for these within an overall solution quality assessment. Finally, a vision of a new information architecture and interactive scienti.c interface is introduced and advocated.

Thompson, David E.↗

Computational Results for the KTH-NASA Wind-Tunnel Model Used for Acquisition of Transonic Nonlinear Aeroelastic Data

A status report is provided on the collaboration between the Royal Institute of Technology (KTH) in Sweden and the NASA Langley Research Center regarding the aeroelastic analyses of a full-span fighter configuration wind-tunnel model. This wind-tunnel model was tested in the Transonic Dynamics Tunnel (TDT) in the summer of 2016. Large amounts of data were acquired including steady/unsteady pressures, accelerations, strains, and measured dynamic deformations. The aeroelastic analyses presented include linear aeroelastic analyses, CFD steady analyses, and analyses using CFD-based reduced-order models (ROMs).

Silva, Walter A.↗

Effect of Turbulence Modeling on an Excited Jet

The flow dynamics in a high-speed jet are dominated by unsteady turbulent flow structures in the plume. Jet excitation seeks to control these flow structures through the natural instabilities present in the initial shear layer of the jet. Understanding and optimizing the excitation input, for jet noise reduction or plume mixing enhancement, requires many trials that may be done experimentally or computationally at a significant cost savings. Numerical simulations, which model various parts of the unsteady dynamics to reduce the computational expense of the simulation, must adequately capture the unsteady flow dynamics in the excited jet for the results are to be used. Four CFD methods are considered for use in an excited jet problem, including two turbulence models with an Unsteady Reynolds Averaged Navier-Stokes (URANS) solver, one Large Eddy Simulation (LES) solver, and one URANS/LES hybrid method. Each method is used to simulate a simplified excited jet and the results are evaluated based on the flow data, computation time, and numerical stability. The knowledge gained about the effect of turbulence modeling and CFD methods from these basic simulations will guide and assist future three-dimensional (3-D) simulations that will be used to understand and optimize a realistic excited jet for a particular application.

Brown, Clifford A.↗

A method for the modelling of porous and solid wind tunnel walls in computational fluid dynamics codes

Porous wall wind tunnels have been used for several decades and have proven effective in reducing wall interference effects in both low speed and transonic testing. They allow for testing through Mach 1, reduce blockage effects and reduce shock wave reflections in the test section. Their usefulness in developing computational fluid dynamics (CFD) codes has been limited, however, by the difficulties associated with modelling the effect of a porous wall in CFD codes. Previous approaches to modelling porous wall effects have depended either upon a simplified linear boundary condition, which has proven inadequate, or upon detailed measurements of the normal velocity near the wall, which require extensive wind tunnel time. The current work was initiated in an effort to find a simple, accurate method of modelling a porous wall boundary condition in CFD codes. The development of such a method would allow data from porous wall wind tunnels to be used more readily in validating CFD codes. This would be beneficial when transonic validations are desired, or when large models are used to achieve high Reynolds numbers in testing. A computational and experimental study was undertaken to investigate a new method of modelling solid and porous wall boundary conditions in CFD codes. The method utilized experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution was then imposed as a pressure boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the effect of porosity changes on the wall was investigated. Also, the effectiveness of this method when only sparse experimental measurements were available has been investigated. The current work demonstrated this approach for low speed flows and compared the results with experimental data obtained from a heavily instrumented variable porosity test section. The approach developed was simple, computationally inexpensive, and did not require extensive or intrusive measurements of the boundary conditions during the wind tunnel test. It may be applied to both solid and porous wall wind tunnel tests.

Beutner, Thomas John↗

A Modification and Analysis of Lagrangian Trajectory Modeling and Granular Dynamics of Lunar Dust Particles

A previously developed mathematical model is amended to more accurately incorporate the effects of lift and drag on single dust particles in order to predict their behavior in the wake of high velocity gas flow. The model utilizes output from a CFD or DSMC simulation of exhaust from a rocket nozzle hot gas jet. An extension of the Saffman equation for lift based on the research of McLaughlin (1991) and Mei (1992) is used, while an equation for the Magnus force modeled after the work of Oesterle (1994) and Tsuji et al (1985) is applied. A relationship for drag utilizing a particle shape factor (phi = 0.8) is taken from the work of Haider and Levenspiel (1989) for application to non-spherical particle dynamics. The drag equation is further adjusted to account for rarefaction and compressibility effects in rarefied and high Mach number flows according to the work of Davies (1945) and Loth (2007) respectively. Simulations using a more accurate model with the correction factor (Epsilon = 0.8 in a 20% particle concentration gas flow) given by Richardson and Zaki (1954) and Rowe (1961) show that particles have lower ejection angles than those that were previously calculated. This is more prevalent in smaller particles, which are shown through velocity and trajectory comparison to be more influenced by the flow of the surrounding gas. It is shown that particles are more affected by minor changes to drag forces than larger adjustments to lift forces, demanding a closer analysis of the shape and behavior of lunar dust particles and the composition of the surrounding gas flow.

Long, Jason M.↗

Three-dimensional numerical predictions of the flow behind a rearward-facing-step in a supersonic combustor

A CFD code was developed to compute the mixing and combustion of hydrogen fuel in the turbulent flow fields of scramjets. The code is based on the complete Reynold's time-averaged three-dimensional Navier-Stokes equations, and chemical species equations. The code includes a global finite rate hydrogen-air chemical reaction model and the Baldwin-Lomax (1978) turbulence model. The validity of this CFD code was demonstrated by comparing computed results from the code with existing data for three configurations, each containing an important feature of the scramjet combustor flow. The code was then applied to the reacting flow field in a complex combustor model in which hydrogen fuel was injected transversely into a supersonic hot airstream through several fuel injectors.

Uenishi, K.↗

Evaluation of grid generation technologies from an applied perspective

An analysis of the grid generation process from the point of view of an applied CFD engineer is given. Issues addressed include geometric modeling, structured grid generation, unstructured grid generation, hybrid grid generation and use of virtual parts libraries in large parametric analysis projects. The analysis is geared towards comparing the effective turn around time for specific grid generation and CFD projects. The conclusion was made that a single grid generation methodology is not universally suited for all CFD applications due to both limitations in grid generation and flow solver technology. A new geometric modeling and grid generation tool, CFD-GEOM, is introduced to effectively integrate the geometric modeling process to the various grid generation methodologies including structured, unstructured, and hybrid procedures. The full integration of the geometric modeling and grid generation allows implementation of extremely efficient updating procedures, a necessary requirement for large parametric analysis projects. The concept of using virtual parts libraries in conjunction with hybrid grids for large parametric analysis projects is also introduced to improve the efficiency of the applied CFD engineer.

Hufford, Gary S.↗

PDF approach for compressible turbulent reacting flows

The objective of the present work is to develop a probability density function (pdf) turbulence model for compressible reacting flows for use with a CFD flow solver. The probability density function of the species mass fraction and enthalpy are obtained by solving a pdf evolution equation using a Monte Carlo scheme. The pdf solution procedure is coupled with a compressible CFD flow solver which provides the velocity and pressure fields. A modeled pdf equation for compressible flows, capable of capturing shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed, and an averaging procedure is developed to provide smooth Monte-Carlo solutions to ensure convergence. Two supersonic diffusion flames are studied using the proposed pdf model and the results are compared with experimental data; marked improvements over CFD solutions without pdf are observed. Preliminary applications of pdf to 3D flows are also reported.

Hsu, A. T.↗