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,009 records · Page 56

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↗

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.

Ajmani, Kumud↗

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

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

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

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↗

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↗

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

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↗

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↗

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↗

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↗

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↗

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↗

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

Modeling the impact of the fuel injection strategy on the combustion and performance characteristics of a heavy-duty GCI engine

Gasoline compression ignition (GCI) is a promising strategy to achieve high thermal efficiency and low emissions with limited modifications to the conventional diesel engine hardware. It is a partially premixed concept, which derives its superiority from higher volatility and longer ignition delay of gasoline-like fuels combined with higher compression ratio typical of diesel engines. Here, the present study investigates the combustion process in a GCI engine operating with different injection strategies using computational fluid dynamics (CFD). Simulations are carried out on a single cylinder of a multi cylinder heavy-duty compression ignition engine, which operates at a compression ratio of 17:1 and an engine speed of 1038 rev/min. Two different injection strategies viz., late injection (LI), and early pilot injection (EP) are investigated to understand their impact on combustion and performance of the engine. Renormalized group (RNG) k-ε model is used to describe in-cylinder turbulence and KH RT model is used to simulate the fuel spray breakup. The developed CFD methodology is validated against relevant experimental data under a wide range of operating conditions for each injection strategy. The developed CFD methodology was found to capture the engine combustion behavior quite well. Based on the validated CFD model, the differences in the progress of combustion event for the two injection strategies is highlighted. It was found that a larger pilot fuel mass fraction results in a steeper rise in the initial heat release rate which in turn influences the transition to mixing controlled combustion. In line with the experimental data, the study showed that the late pilot injection strategy with three injection pulses, results in higher performance compared to the other conditions.

33 ADVANCED PROPULSION SYSTEMS↗

Validation of Two-Phase Flow Modeling Techniques in Loci-Stream Using Axial Jet Pressure Control and Autogenous Pressurization Experiments

Improvements to computational fluid dynamics modeling of cryogenic fluid management operations are presented with experimental validations. Modeling improvements includes utilizing improved delayed detaches eddy simulation to more effectively model turbulence compared to the Reynolds-averaged Navier Stokes modeling, modeling temperature dependent tank solid properties to capture large variation in thermodynamic properties at cryogenic conditions, modeling real gas equation of state for hydrogen instead of ideal gas behavior, and an on-going effort to improve interface turbulence transport for the volume of fluid method. Two flight-scale ground tests, five K-site tank axial jet pressure control experiments and an Engineering Design Unit tank autogenous pressurization experiment are validated and show the effectiveness of the improved modeling techniques. The successful validation establishes further trust in these models for flight project support and provides a pathway for maturation of in-development interface turbulence damping model for the volume of fluid method for more generalized applications.

CFD↗

Computational Methods for HSCT-Inlet Controls/CFD Interdisciplinary Research

A program aimed at facilitating the use of computational fluid dynamics (CFD) simulations by the controls discipline is presented. The objective is to reduce the development time and cost for propulsion system controls by using CFD simulations to obtain high-fidelity system models for control design and as numerical test beds for control system testing and validation. An interdisciplinary team has been formed to develop analytical and computational tools in three discipline areas: controls, CFD, and computational technology. The controls effort has focused on specifying requirements for an interface between the controls specialist and CFD simulations and a new method for extracting linear, reduced-order control models from CFD simulations. Existing CFD codes are being modified to permit time accurate execution and provide realistic boundary conditions for controls studies. Parallel processing and distributed computing techniques, along with existing system integration software, are being used to reduce CFD execution times and to support the development of an integrated analysis/design system. This paper describes: the initial application for the technology being developed, the high speed civil transport (HSCT) inlet control problem; activities being pursued in each discipline area; and a prototype analysis/design system in place for interactive operation and visualization of a time-accurate HSCT-inlet simulation.

Cole, Gary L.↗