SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Role of Inlet Boundary Conditions on Fuel-Air Mixing at Supercritical Conditions

Advanced engine design and alternative fuels present the possibility of fuel injection at purely supercritical conditions in diesel engines and gas turbines. The complex interactions that govern this phenomenon still need significant research, particularly the boundary conditions for fuel injection are critical for accurate simulation. However, the flow inside the injector itself is often omitted to reduce the computational efforts, and thus, velocity, mass flux, or total pressure is specified at the injector exit (or domain inlet), often with simplified velocity profiles and turbulence levels. This simplified inlet boundary treatment has minimal effects on results for conventional fuel injection conditions, however, the validity of this approach at supercritical conditions has not been assessed. Comprehensive real-gas and binary fluid mixing models have been implemented for computational fluid dynamics (CFD) analysis of fuel-air mixing at supercritical conditions. The model is verified using prior CFD results from the literature. The model is used to investigate the effects of the shape of axial velocity and mass fraction profiles at the inlet boundary with the goal to improve the comparison of predictions to experimental data. Results show that the boundary conditions have a significant effect on the predictions, and none of the cases match precisely with experimental data. Furthermore, the study reveals that the physical location of the inlet boundary might be difficult to infer correctly from the experiments and highlights the need for high-quality, repeatable measurements at supercritical conditions to support the development of relevant high-fidelity models for fuel-air mixing.

02 PETROLEUM↗

NASA NTRS2007

PAB3D Simulations for the CAWAPI F-16XL

Numerical simulations of the flow around F-16XL are performed as a contribution to the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) using the PAB3D CFD code. Two turbulence models are used in the calculations: a standard k-! model, and the Shih-Zhu-Lumley (SZL) algebraic stress model. Seven flight conditions are simulated for the flow around the F-16XL where the free stream Mach number varies from 0.242 to 0.97. The range of angles of attack varies from 0° to 20°. Computational results, surface static pressure, boundary layer velocity profiles, and skin friction are presented and compared with flight data. Numerical results are generally in good agreement with flight data, considering that only one grid resolution is utilized for the different flight conditions simulated in this study. The ASM results are closer to the flight data than the k-ε model results. The ASM predicted a stronger primary vortex, however, the origin of the vortex and footprint is approximately the same as in the k-ε predictions.

Flight Data↗

NASA NTRSDate not supplied

Exploring data-driven modeling of boundary layer transition

Prediction of laminar-turbulent transition in boundary layer flows is an important component of predicting the aerodynamic performance of a number of aerospace configurations. According to the CFD Vision 2030 [1], transition modeling represents acriticalarea in CFD simulation capability that will remain a pacing item for the foreseeable future. The fact thattransition can take placevia either one of a myriad possible paths adds to the challenges inreliable transition predictions, despite a limited knowledge of the relevant input parameters. In the low disturbance environments typical of flight applications, transition is often initiated by small amplitude disturbances in the form of linear instability waves of the laminar boundary layer. These disturbances amplify linearly at first and eventually undergo a sequence of nonlinear interactions that result in transition to turbulence. Because the nonlinear phase is rather rapid, the amplification of boundary layer instabilities is governed by the linearstability theory over a majority of the distance leading up to the onset of transition. Semi-empirical transition correlations based on the linear stability theory have been successful in explaining the observed trends in transition location within a broad class of flows. However, the application of stability theory is highly non-robust and often requires a significant domain expertise. Recent work at the NASA Langley Research Center has beenaimed at bridging the gap between physics based transition analyses such as those based on linear stability theory and practical applications that require transition prediction by users that may not be well versed in transition physics. The applications of deep learning have been at the center of these efforts. This presentation will focus on the progress achieved thus far, highlighting the applications of neural networks to selectedtransition scenarios across a range of Mach numbers and flow configuration, as well as the lessons learnedand remaining challengeswithrespect to the selection of training data and neural networks architectures, hyperparameter tuning, and the physical insights distilled from the otherwise black-box models.

M. R. Malik↗

NASA NTRS2019

CFD Evaluation of Lean-Direct Injection Combustors for Commercial Supersonics Technology

An overview is given of an effort that focused on using CFD analysis to complement design and configuration definition of Lean-Direct Injection (LDI) combustion concepts for NASA's Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for second and third generation LDI configurations at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray modeling, ignition and kinetics modeling. Emissions (EINOx) characteristics were predicted for CST cruise conditions, and compared with emissions data from experimental measurements to evaluate the fidelity of the CFD modeling approach to predict emissions changes in response to changes in supersonic cycle conditions.

nitrogen oxides injectors↗

NASA NTRS2019

CFD Evaluation Of Lean-Direct Injection Combustors for Commercial Supersonics Technology

Ajmani, Kumud↗

NASA NTRS2006

Estimating Resolution Lengths of Hybrid Turbulence Models

A two-stage procedure has been devised for estimating the spatial resolution achievable in the simulation of a given flow on a given computational grid by a computational fluid dynamics (CFD) code that incorporates a hybrid model of turbulence. The hybrid models to which this procedure is especially relevant are those of the Reynolds-averaged Navier-Stokes (RANS) and the partial-averaged Navier-Stokes (PANS) approaches. This procedure represents the first step toward adding variable-resolution turbulence-modeling capabilities to CFD codes as part of a continuing effort to increase the accuracy and robustness of CFD simulations of unsteady flows. Some background information is prerequisite to a meaningful summary of the procedure. Among experts in CFD, it is well known that combination of the Reynolds-averaged Navier-Stokes (RANS) approach and eddy-viscosity turbulence models offers limited capability for simulating unsteady and complex flows. The RANS approach includes an assumption that most of the energy in a given flow is modeled through turbulence-transport equations and is resolved in a computational grid used to simulate the flow. RANS also overpredicts eddy viscosity, thereby yielding excessive damping of unsteady motion. The eddy viscosity attains an unphysically large value because of unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved flow field. One approach used to overcome this deficiency is to provide a mechanism for the RANS equations to resolve motion only on the largest scales and to use a hybrid model to represent effects at smaller scales. The RANS approach involves the use of a standard two-equation turbulence model in which the effect of turbulence is summarized by a viscosity that is a function of (1) the time-averaged kinetic- energy density (k) associated with the local fluctuating (turbulent) component of flow and (2) the time-averaged rate of dissipation of the turbulent-kinetic- energy density ( ). In PANS, which was developed to overcome the grid dependency associated with other hybrid turbulence models (including that of RANS), the standard two-equation turbulence model is replaced by another two-equation model in which one solves for the previously unresolved k and , which are now allowed to vary in space and time.

Abdol-Hamid, Khaled S.↗

NASA NTRS1999

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool for the analysis of flowfields near the injector face has been developed and used to analyze the Fastrac engine. The objective was to produce a computational design tool which was detailed enough to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows. To obtain a model which could be used to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe liquid and vapor sub- and super-critical flows, the model included thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. A homogeneous model was constructed such that the local state of the flow was determined directly, i.e. the quality of the flow was calculated. Such a model does not identify drops or their distribution, but it does allow the flow along the injector face and into the acoustic cavity to be predicted. It also allows the film coolant flow to be accurately described. The initial evaluation of the injector code was made by simulating cold flow from an unlike injector element and from a like-on-like overlapping fan (LOL) injector element. The predicted mass flux distributions of these injector elements compared well to cold flow test results. These are the same cold flow tests which serve as the data base for the JANNAF performance prediction codes. The flux distributions 1 inch downstream of the injector face are very similar; the differences were somewhat larger at further distances from the faceplate. Since the cold flow testing did not achieve good mass balances when integrations across the entire fan were made, the CFD simulation appears to be reasonable alternative to future cold flow testing. To simulate the Fastrac, an RP-1/LOX combustion model must be chosen. This submodel must be relatively simple to accomplish three-dimensional, multiphase flow simulations. Single RP-1 pyrolysis and partial oxidation steps were chosen and the combustion was completed with the wet CO mechanism. Soot was also formed with a single global reaction. To validate the combustion submodel, global data from gas generator tests and from subscale motor test were used to predict qualitatively correct mean molecular weights, temperature, and soot levels. Because such tests do not provide general kinetics rates, the methodology is not necessarily appropriate for other than rocket type flows conditions. Soot predictions were made so that radiation heating to the motor walls can be made. These initial studies of the Fastrac were for a small region close to the injector face and chamber wall which included a segment of the acoustic cavity. The region analyzed includes 11 individual orifice holes to represent the LOL elements and the H2 film coolant holes. Typical results of this simulation are shown in Figure 1. At this point the only available test data to verify the predictions are temperatures measured in the acoustic cavity. These temperatures are in reasonable agreement at about 2000R (1111 K). Future work is expected to include improving the computational efficiency or the CFD model and/or using more computer capacity than the single Pentium PC with which these simulations were made.

Farmer, Richard C.↗

DOE OSTI2022

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗

NASA NTRSDate not supplied

Comparing Planetary Helicopter Brownout Observations with Numerical Predictions for Ingenuity

Ingenuity is considered the Martian equivalent of the Wright Flyer. Its sustained flights on Mars proved that non-terrestrial powered flight is possible, even on planets with thin atmospheres. However, Ingenuity was small and carried no payload. Future planetary rotorcraft will be bigger and heavier to accomplish new scientific and exploratory goals. On dusty planets, the associated increase in rotor thrust poses a potential challenge well known to helicopter pilots who have flown in desert environments: dust clouds forming during take-off and landing due to rotor wash, referred to as brownout. This work evaluates the effectiveness of existing numerical simulation models in predicting the severity of this phenomenon for planetary rotorcraft. We compare the output of a CFD-based dust mobilization model with data recently gathered using image processing techniques of Ingenuity’s flights. The simulation results suggest that the sensitivity of the wall friction velocity to the rotorcraft height may be lower than previously assessed. We further evaluate the sensitivity of the total mass of mobilized dust to other model parameters such as the sandblasting efficiency and the saltation threshold velocity. Finally, we discuss the limitations of this comparison due to uncertainty in both the observation and the numerical model.

HEC↗

NASA NTRS2024

Measurements and Computational Analysis of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

The NASA Juncture Flow Experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This report presents and discusses the results of a recent experiment with the Juncture Flow Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of the test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model angle of incidence. Laser Doppler velocimetry (LDV) measurements were made at three model angles of incidence (0° : fully attached, 1° : incipient separation, and 5° : separated flow) and for each angle, mean-flow and Reynolds-stress data were obtained on the fuselage and at several locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and include model and tunnel wall static pressures, tunnel wall boundary-layer profiles, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. A companion set of Reynolds-averaged Navier-Stokes (RANS) CFD simulations were also performed and those results were used to support the discussion of the experimental test results.

Juncture Flow↗

NASA NTRS2010

Direct CFD Predictions of Low Frequency Sounds Generated by a Helicopter Main Rotor

The use of CFD to directly predict helicopter main rotor noise is shown to be quite promising as an alternative mean for low frequency source noise evaluation. Results using existing state-of-the-art grid structures and finite-difference schemes demonstrated that small perturbation pressures, associated with acoustics radiation, can be extracted with some degree of fidelity. Accuracy of the predictions are demonstrated via comparing to predictions from conventional acoustic analogy-based models, and with measurements obtained from wind tunnel and flight tests for the MD-902 helicopter at several operating conditions. Findings show that the direct CFD approach is quite successfully in yielding low frequency results due to thickness and steady loading noise mechanisms. Mid-to-high frequency contents, due to blade-vortex interactions, are not predicted due to CFD modeling and grid constraints.

Sim, Ben W.↗

NASA NTRSDate not supplied

Multiscale Mesh Adaptation for Transonic Aeroelastic Flutter Problems

This work applies multiscale mesh adaptation with refine to reduce spatial discretization error of aeroelastic computational fluid dynamics (CFD) simulations. Benchmark flutter models, such as the pitch and plunge NACA64A-010 airfoil and the benchmark supercritical wing, are studied with both a linearized frequency-domain solver and time-marching CFD coupled to a modal structural solver in FUN3D. The undeformed NASA Common Research Model (CRM), an aeroelastic jig shape variant of the CRM, is also studied with the linearized frequency-domain approach. For these cases, the adaptation process converges to comparable flutter predictions to hand-generated meshes but with smaller node counts. However the additional disciplines of the linearized frequency-domain analysis, the mesh deformation, and the unsteady finite-volume solver create robustness challenges that need to be addressed before it can be applied as a fully automated process for complex transonic aeroelastic problems. In particular, negative volumes are observed to be an issue for FUN3D’s linear elasticity mesh deformation solver when moving the adapted meshes.

Aeroelasticity↗

NASA NTRS2023

Measurements and Computations of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

The NASA Juncture Flow experiment is designed to acquire high-quality flowfield data deep in the corner of a wing-fuselage junction specifically for the purpose of computational fluid dynamics (CFD) validation and turbulence model improvement. This paper will present and discuss the results of a recent experiment with the juncture-flow model in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The main objective of that test was to expand the existing juncture-flow dataset with a symmetric wing case that displays fully attached, incipient separation, and separated flow in the corner of the wing-fuselage junction, depending on the model pitch angle. Laser Doppler velocimetry (LDV) measurements were made at three model pitch angles (0 deg: fully attached, 1 deg: incipient separation, and 5 deg: separated flow) and for each one, mean-flow and Reynolds-stress data was obtained on the fuselage and at several streamwise locations along the corner of the wing-fuselage junction. Supporting measurements were made during the test campaign and included model and tunnel wall static pressures, tunnel wall and ceiling boundary-layer rake data, oil-flow visualizations, and laser-based measurements of the as-built model geometry and model position in the test section. Comparisons between the experimental data and Reynolds averaged Navier-Stokes CFD results will also be presented and discussed.

Juncture Flow↗

NASA NTRSDate not supplied

Measurements and Computations of the Turbulent Corner Flow on the NASA Juncture-Flow Model with a Symmetric Wing

Juncture Flow↗

DOE OSTI2019

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗

NASA NTRS2013

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: 1) a k-epsilon turbulence model computation on a 6.8 million point grid using wall functions, 2) a k-epsilon turbulence model computation on a 14 million point grid integrating to the wall, and 3) a k-omega turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-epsilon and k-omega computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-epsilon computations which extended from 70% to 100% span at the exit rating plane, whereas the k-omega computation had reversed flow from 95% to 100% span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-epsilon case resulted in an under-prediction in adiabatic efficiency of 8.3 points, whereas the k-omega case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

NASA NTRS2013

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: (1) a k-ε turbulence model computation on a 6.8 million point grid using wall functions, (2) a k-ε turbulence model computation on a 14 million point grid integrating to the wall, and (3) a k-ω turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-ε and k-ω computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-ε computations which extended from 70 to 100 percent span at the exit rating plane, whereas the k-ω computation had reversed flow from 95 to 100 percent span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-ε case resulted in an underprediction in adiabatic efficiency of 8.3 points, whereas the k-ω case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

NASA NTRS1995

Turbulence modeling for non-equilibrium flow

The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Durbin, P. A.↗