Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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,171 records · Page 65

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work uses an automated grid adaptation tool recently implemented in VULCAN-CFD, and explores this tool’s ability to solve highspeed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user-generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flowfield parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Furthermore, motivated by the grid convergence index (GCI) used for structured grid simulations, a grid convergence estimate (GCE) was developed and demonstrated for the grid adaptation. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex high-speed mixing problems with only a fraction of the grid cells required to obtain comparable results on a user-generated structured grid.

hypersonics↗

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work utilizes an automated grid adaptation tool recently implemented into VULCAN-CFD, and explores this tool’s ability to solve high-speed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flow field parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex highspeed mixing problems using only a fraction of the grid cells required to obtain comparable results using a user-generated structured grid.

grid adaptation↗

A numerical investigation of the finite element method in compressible primitive variable Navier-Stokes flow

The results of a comprehensive numerical investigation of the basic capabilities of the finite element method (FEM) for numerical solution of compressible flow problems governed by the two-dimensional and axis-symmetric Navier-Stokes equations in primitive variables are presented. The strong and weak points of the method as a tool for computational fluid dynamics are considered. The relation of the linear element finite element method to finite difference methods (FDM) is explored. The calculation of free shear layer and separated flows over aircraft boattail afterbodies with plume simulators indicate the strongest assets of the method are its capabilities for reliable and accurate calculation employing variable grids which readily approximate complex geometry and capably adapt to the presence of diverse regions of large solution gradients without the necessity of domain transformation.

Cook, C. H.↗

Simulations of the Unsteady Flow through the Fastrac Supersonic Turbine

Analysis of the unsteady aerodynamic environment in the Fastrac supersonic turbine is presented. Modal analysis of the turbine blades indicated possible resonance in crucial operating ranges of the turbopump. Unsteady computational fluid dynamics (CFD) analysis was conducted to support the aerodynamic and structural dynamic assessments of the turbine. Before beginning the analysis, two major problems with current unsteady analytical capabilities had to be addressed: modeling a straight centerline nozzle with the turbine blades and exit guide vanes (EGVs), and reducing run times significantly while maintaining physical accuracy. Modifications were made to the CFD code used in this study to allow the coupled nozzle/blade/EGV analysis and to incorporate Message Passing Interface (MPI) software. Because unsteadiness is a key issue for the Fastrac turbine [and future rocket engine turbines such as for the Reusable Launch Vehicle (RI.V)], calculations were performed for two nozzle-to-blade axial gaps. Calculations were also performed for the nozzle alone, and the results were imposed as an inlet boundary condition for a blade/EGV calculation for the large gap case. These results are compared to the nozzle/blade/EGV results.

Griffin, Lisa W.↗

Comparison of a 3-D CFD-DSMC Solution Methodology With a Wind Tunnel Experiment

A solution method for problems that contain both continuum and rarefied flow regions is presented. The methodology is applied to flow about the 3-D Mars Sample Return Orbiter (MSRO) that has a highly compressed forebody flow, a shear layer where the flow separates from a forebody lip, and a low density wake. Because blunt body flow fields contain such disparate regions, employing a single numerical technique to solve the entire 3-D flow field is often impractical, or the technique does not apply. Direct simulation Monte Carlo (DSMC) could be employed to solve the entire flow field; however, the technique requires inordinate computational resources for continuum and near-continuum regions, and is best suited for the wake region. Computational fluid dynamics (CFD) will solve the high-density forebody flow, but continuum assumptions do not apply in the rarefied wake region. The CFD-DSMC approach presented herein may be a suitable way to obtain a higher fidelity solution.

Glass, Christopher E.↗

On the Application of an Actuator Line Model for Rotorcraft Outwash Predictions

As the next generation of Vertical Take-Off and Landing (VTOL) vehicles develops and the VTOL industry makes significant progress in utilizing them for public transportation, it is crucial to have validated and cost-effective computational models accessible to both industry professionals and academics in the rotorcraft community. This study investigates the implementation of a new actuator line model within NASA’s OVERFLOW computational fluid dynamics (CFD) solver, with a focus on predicting rotorcraft outwash. This reduced-order rotor model was developed to provide reasonably accurate outwash predictions while reducing computational costs compared to traditional blade-resolved CFD simulations. To validate the actuator line model for outwash predictions, the model was compared against existing experimental data across three validation cases: single rotor Out-of-Ground Effect (OGE), single rotor In-Ground Effect (IGE), and tandem rotor IGE. All three cases were based on the rotor geometry and configuration of the CH-47D. To evaluate the feasibility and advantages of using an actuator line model for outwash predictions the results were compared to a series of high-fidelity blade-resolved CFD simulations as well as comprehensive analyses simulations based on CHARM predictions. The results of this study demonstrate the feasibility of leveraging this new actuator line model for reasonably accurate outwash predictions. The outwash results obtained from the actuator line model closely matched those from high-fidelity blade-resolved simulations while significantly reducing the overall computational costs of the simulations. Results from this study further demonstrated the feasibility of using CHARM for efficient, and reasonably accurate time-averaged outwash predictions for multi-rotor configurations.

CFD↗

Sensitivity of Age-of-Air Calculations to the Choice of Advection Scheme

The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been interpreted in terms of the models' large-scale circulation field. This study shows, however, that in addition to the simulated large-scale circulation, three-dimensional age calculations can also be affected by the choice of advection scheme employed in solving the tracer continuity equation, Specifically, using the 3.0deg latitude X 3.6deg longitude and 40 vertical level version of the Geophysical Fluid Dynamics Laboratory SKYHI GCM and six online transport schemes ranging from Eulerian through semi-Lagrangian to fully Lagrangian, it will be demonstrated that the oldest ages are obtained using the nondiffusive centered-difference schemes while the youngest ages are computed with a semi-Lagrangian transport (SLT) scheme. The centered- difference schemes are capable of producing ages older than 10 years in the mesosphere, thus eliminating the "young bias" found in previous age-of-air calculations. At this stage, only limited intuitive explanations can be advanced for this sensitivity of age-of-air calculations to the choice of advection scheme, In particular, age distributions computed online with the National Center for Atmospheric Research Community Climate Model (MACCM3) using different varieties of the SLT scheme are substantially older than the SKYHI SLT distribution. The different varieties, including a noninterpolating-in-the-vertical version (which is essentially centered-difference in the vertical), also produce a narrower range of age distributions than the suite of advection schemes employed in the SKYHI model. While additional MACCM3 experiments with a wider range of schemes would be necessary to provide more definitive insights, the older and less variable MACCM3 age distributions can plausibly be interpreted as being due to the semi-implicit semi-Lagrangian dynamics employed in the MACCM3. This type of dynamical core (employed with a 60-min time step) is likely to reduce SLT's interpolation errors that are compounded by the short-term variability characteristic of the explicit centered-difference dynamics employed in the SKYHI model (time step of 3 min). In the extreme case of a very slowly varying circulation, the choice of advection scheme has no effect on two-dimensional (latitude-height) age-of-air calculations, owing to the smooth nature of the transport circulation in 2D models. These results suggest that nondiffusive schemes may be the preferred choice for multiyear simulations of tracers not overly sensitive to the requirement of monotonicity (this category includes many greenhouse gases). At the same time, age-of-air calculations offer a simple quantitative diagnostic of a scheme's long-term diffusive properties and may help in the evaluation of dynamical cores in multiyear integrations. On the other hand, the sensitivity of the computed ages to the model numerics calls for caution in using age of air as a diagnostic of a GCM's large-scale circulation field.

Eluszkiewicz, Janusz↗

Flow Characterization of the NASA Langley Unitary Plan WindTunnel, Test Section 2: Computational Results

Flow in the empty Unitary Planform Wind Tunnel at the NASA Langley Research Center issimulated with computational fluid dynamics methods. The objectives are to assess CFD’s truepredictive capability and to generate flow-field maps upstream of the tunnel’s test section foruse in associated vehicle tests. Multiple CFD solvers, grid adaption methods, and turbulencemodels are used by five teams doing the simulations. The simulation domain is as large aspractical: from the start of the settling chamber, just downstream of the last set of turningvanes, to well downstream of the test section. Simulations are done at three Mach numbers andseveral Reynolds numbers that cover most of the tunnel’s operating range. The sensitivity ofthe CFD solutions to simulation process details, including turbulence model, spatial resolutionand settling-chamber inflow velocity profile, is characterized. Significant sensitivities in somedetails of the predicted test section flow are described. CFD predicts that details of flow inthe settling chamber can affect the flow in the test section. Turning of the flow through thetunnel contraction and past nozzle-block hardware generates streamwise vortices, some ofwhich survive into the core of the test section, principally at low Mach number. Preliminarycomparisons to experimental measurements are given. The major flow characteristics andthe dynamical vortices predicted by CFD exist in the experiment. Several details in CFD andexperiment are observed to differ, including the effects of the vortices and features that appearto be Mach waves in the test section. Sensitivity to grid resolution and turbulence modelingis noted in the vortices; no potential cause for the differences in Mach waves has yet beenidentified in the CFD sensitivity studies.

wind tunnel↗

Modeling Laminar-to-Turbulent Transition in the Panel Test Facility Arcjet

A computational fluid dynamics (CFD) laminar-to-turbulence transition model was developed for the NASA Ames Research Center’s 20 MW Panel Test Facility (PTF). Surface pressure, heat flux to a water-cooled plate, and surface temperature on a tile plate coated with reaction-cured glass were measured across several conditions in the facility and compared with laminar and fully-turbulent CFD simulations. The potential for bypass transition in the PTF nozzle was assessed via application of the Langtry-Menter four-equation transitional shear-stress transport (SST) model. Results from the transition model were inconsistent with measurements. Flow interaction with a boundary conditioning plate feature inside the nozzle was also investigated as a potential source of laminar-to-turbulent transition using two turbulence models with specified transition locations. The Baldwin-Lomax turbulence model was configured to simulate a transition at the upstream edge of the boundary conditioning plate and produced results consistent with the surface pressure measurements but not the cold-wall heat flux. Finally, the SST turbulence model was calibrated to transition at the upstream edge of the boundary conditioning plate and produced results consistent with both the surface pressure and cold-wall heat flux measurements. The SST-based model demonstrated reasonable agreement with surface temperature measurements on the reaction-cured glass tile, albeit with some discrepancies.

Computational Fluid Dynamics↗

Unsteady Three-Dimensional Simulation of a Shear Coaxial GO2/GH2 Rocket Injector with RANS and Hybrid-RAN-LES/DES Using Flamelet Models

Historically, the analysis and design of liquid rocket engines (LREs) has relied on full‐scale testing and one‐dimensional empirical tools. The testing is extremely expensive and the one‐dimensional tools are not designed to capture the highly complex, and multi‐dimensional features that are inherent to LREs. Recent advances in computational fluid dynamics (CFD) tools have made it possible to predict liquid rocket engine performance, stability, to assess the effect of complex flow features, and to evaluate injector‐driven thermal environments, to mitigate the cost of testing. Extensive efforts to verify and validate these CFD tools have been conducted, to provide confidence for using them during the design cycle. Previous validation efforts have documented comparisons of predicted heat flux thermal environments with test data for a single element gaseous oxygen (GO2) and gaseous hydrogen (GH2) injector. The most notable validation effort was a comprehensive validation effort conducted by Tucker et al. [1], in which a number of different groups modeled a GO2/GH2 single element configuration by Pal et al [2]. The tools used for this validation comparison employed a range of algorithms, from both steady and unsteady Reynolds Averaged Navier‐Stokes (U/RANS) calculations, large‐eddy simulations (LES), detached eddy simulations (DES), and various combinations. A more recent effort by Thakur et al. [3] focused on using a state‐of‐the‐art CFD simulation tool, Loci/STREAM, on a two‐dimensional grid. Loci/STREAM was chosen because it has a unique, very efficient flamelet parameterization of combustion reactions that are too computationally expensive to simulate with conventional finite‐rate chemistry calculations. The current effort focuses on further advancement of validation efforts, again using the Loci/STREAM tool with the flamelet parameterization, but this time with a three‐dimensional grid. Comparisons to the Pal et al. heat flux data will be made for both RANS and Hybrid RANSLES/ Detached Eddy simulations (DES). Computation costs will be reported, along with comparison of accuracy and cost to much less expensive two‐dimensional RANS simulations of the same geometry.

Westra, Doug G.↗

Mixing and NOx Emission Calculations of Confined Reacting Jet Flows in Cylindrical and Annular Ducts

Rapid mixing of cold lateral jets with hot cross-stream flows in confined configurations is of practical interest in gas turbine combustors as it strongly affects combustor exit temperature quality, and gaseous emissions in for example rich-lean combustion. It is therefore important to further improve our fundamental understanding of the important processes of dilution jet mixing especially when the injected jet mass flow rate exceeds that of the cross-stream. The results reported in this report describe some of the main flow characteristics which develop in the mixing process in a cylindrical duct. A three-dimensional computational fluid dynamics (CFD) code has been used to predict the mixing flow field characteristics and NOx emission in a quench section of a rich-burn/quick-mix/lean-burn (RQL) combustor. Sixty configurations have been analyzed in both circular and annular geometries in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying the number of orifices per row and orifice shape. Other parameters such as J (momentum-flux ratio), MR (mass flowrate ratio), DR (density ratio), and mixer sector orifice ACd (effective orifice area) were maintained constant throughout the entire study. The results indicate that the mixing flow field can be correlated with the NOx production if they are referenced with the stoichiometric equivalence ratio value and not the equilibrium value. The mixing flowfields in both circular and annular mixers are different. The penetration of equal jets in both annular and circular geometries is vastly different which significantly affects the performance of the mixing section. In the computational results with the circular mixer, most of the NOx formation occurred behind the orifice starting at the orifice wake region. General trends have been observed in the NOx production as the number of orifices is changed and this appears to be common for all hole configurations and mixer types (circular or annular). The performance of any orifice shape (in producing minimum NOx) appears to be acceptable if the number of orifices can be freely varied in order to attain the optimum jet penetration.

Oechsle, Victor L.↗

Ocean modelling on the CYBER 205 at GFDL

At the Geophysical Fluid Dynamics Laboratory, research is carried out for the purpose of understanding various aspects of climate, such as its variability, predictability, stability and sensitivity. The atmosphere and oceans are modelled mathematically and their phenomenology studied by computer simulation methods. The present state-of-the-art in the computer simulation of large scale oceans on the CYBER 205 is discussed. While atmospheric modelling differs in some aspects, the basic approach used is similar. The equations of the ocean model are presented along with a short description of the numerical techniques used to find their solution. Computational considerations and a typical solution are presented in section 4.

Cox, M.↗

A survey of simulation and diagnostic techniques for hypersonic nonequilibrium flows

The possible means of simulating nonequilibrium reacting flows in hypersonic environments, and the required diagnostic techniques, are surveyed in two categories: bulk flow behavior and determination of chemical rate parameters. Flow visualization of shock shapes for validation of computational-fluid dynamic calculations is proposed. The facilities and the operating conditions necessary to produce the required nonequilibrium conditions, the suitable optical techniques, and their sensitivity requirements, are surveyed. Shock-tubes, shock-tunnels, and ballistic ranges in a wide range of sizes and strengths are found to be useful for this purpose, but severe sensitivity requirements are indicated for the optical instruments, which can be met only by using highly-collimated laser sources. Likewise, for the determination of chemical parameters, this paper summarizes the quantities that need to be determined, required facilities and their operating conditions, and the suitable diagnostic techniques and their performance requirements. Shock tubes of various strengths are found to be useful for this purpose. Vacuum ultraviolet absorption and fluorescence spectroscopy and coherent anti-Stokes Raman spectroscopy are found to be the techniques best suited for the measurements of the chemical data.

Sharma, Surendra P.↗

Development of the KIVA-2 CFD code for rocket propulsion applications

The KIVA-2 code, originally developed to solve computational fluid dynamics problems in internal combustion engines, has been developed to solve rocket propulsion type flows. The objective of the work was to develop a code such that both liquid and solid particle motion could be simulated for arbitrary geometry and high speed as well as low speed reacting flows. Modification to the original code include: incorporating independently specific supersonic and subsonic inflows and outflows; symmetric as well as periodic boundary conditions; and the capability to use generalized single or multi-specie thermodynamic data and transport coefficients allowing the user to specify arbitrary wall temperature/heat flux distributions. This code has been shown to successfully solve rocket propulsion flows as well as flows with entrained particles for several different rocket nozzles.

Shannon, Robert V., Jr.↗

SSME HPOTP impeller backcavity CFD analysis

The ball bearings behind the Space Shuttle Main Engine (SSME) HPOTP preburner pump have a history of premature wear requiring their replacement. Extensive tests have been conducted in an attempt to identify the operating factors that contribute to the wear. It has been conjectured that the coolant inflow velocity swirl pattern can aid bearing operation by matching ball orbit speed and thus affect bearing life. However, control of the velocity distribution up to now could only be achieved by trial and error following hardware testing. Observation of hardware from recent flight and development operation led to the hypothesis that certain assemblies with more extensive grinding patterns on the backwall of the impeller for rotor balancing correlated with improved bearing wear. To analytically evaluate the effect of cavity configuration on the flowfield, 3-D computational fluid dynamics (CFD) analyses of various geometries was successfully executed using REACT3D. Height of the anti-vortex ribs on the stationary wall was varied, as was the configuration of the rotating wall, from smooth to simulations of various grindout patterns. The results obtained indicate the effects of the various geometries and provide valuable guidelines for cavity modification to optimize bearing cooling.

Hsu, W. W.↗

Validation of CFD for aircraft design: A Boeing view

Information is given in viewgraph form on computational fluid dynamics (CFD) for aircraft design. Topics covered include CFD validation for advanced systems, cavity flow, transonic flow, separated flow, boundary layer interaction, hypersonic flow, heat transfer, zonal modeling, the mathematical foundation for Navier-Stokes simulation, hypersonic inlets, and the role of wind tunnel tests.

Vas, Irwin E.↗

Parallel computational fluid dynamics '91; Conference Proceedings, Stuttgart, Germany, Jun. 10-12, 1991

A conference was held on parallel computational fluid dynamics and produced related papers. Topics discussed in these papers include: parallel implicit and explicit solvers for compressible flow, parallel computational techniques for Euler and Navier-Stokes equations, grid generation techniques for parallel computers, and aerodynamic simulation om massively parallel systems.

Reinsch, K. G.↗

Overview of ICE Project: Integration of Computational Fluid Dynamics and Experiments

Researchers at the NASA Glenn Research Center have developed a prototype integrated environment for interactively exploring, analyzing, and validating information from computational fluid dynamics (CFD) computations and experiments. The Integrated CFD and Experiments (ICE) project is a first attempt at providing a researcher with a common user interface for control, manipulation, analysis, and data storage for both experiments and simulation. ICE can be used as a live, on-tine system that displays and archives data as they are gathered; as a postprocessing system for dataset manipulation and analysis; and as a control interface or "steering mechanism" for simulation codes while visualizing the results. Although the full capabilities of ICE have not been completely demonstrated, this report documents the current system. Various applications of ICE are discussed: a low-speed compressor, a supersonic inlet, real-time data visualization, and a parallel-processing simulation code interface. A detailed data model for the compressor application is included in the appendix.

Stegeman, James D.↗