SEARCH · Engineering Papers
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.
This paper presents detailed description of a novel CFD procedure and comparison of its solution results to that obtained by other available CFD codes as well as actual flight and wind tunnel test data pertaining to the GIII aircraft, currently undergoing flight testing at AFRC.
Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD researchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where simplex elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identifies the reason behind the difficulties in use of such high-aspect ratio simplex elements is formulated using two different approaches and presented here. Drawing insights from the analysis, a potential solution to avoid that pitfall is also provided as part of this work. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, how the gradient evaluation procedures of the CESE framework can be effectively used to produce accurate and stable results on such high-aspect ratio simplex meshes is also showcased.
Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.
This report presents a comprehensive computational fluid dynamics (CFD) study of the biomass pyrolysis vapor phase upgrading and the catalyst regeneration in different reactors under various operational conditions. This study used the open-source software Multiphase Flow with Interphase eXchanges (MFiX), developed at National Energy Technology Laboratory (NETL) for simulating hydrodynamics, heat transfer, and chemical reactions in multiphase systems. This simulation solved the conservation of mass, momentum, energy, species for gas and solid phases with chemical reactions.
Dynamic stability analysis of the Multi-Purpose Crew Vehicle (MPCV) is performed us-ing the US3D flow solver and a plugin library to solve rigid-body dynamics, Free-Flight CFD(FF-CFD). Previous efforts investigated the free-flight behavior of higher speed, open-backaeroshell shapes. Three primary experimental sources are used to evaluate the predictivecapability of the FF-CFD solver in the low supersonic range (Mach≈1). First, the ballisticrange results obtained at the HFFAF facility at NASA Ames is used to verify the six degree-of-freedom (6-DoF) dynamic capability of the FF-CFD solver. Next, FF-CFD simulationsare preformed using restricted motion and the resultant trajectories are post-processed toobtain pitch damping coefficient as a function of angle-of-attack. The comparison of thepitch damping results obtained with one degree-of-freedom (1-DoF) FF-CFD compare wellwith experimental fits derived from ballistic range data. Finally, the atmospheric flight ca-pability of FF-CFD is compared to data from the Ascent Abort-2 (AA-2) flight experiment.Two simulations were performed using forced and forced-free flight to investigate surfacepressure predictive capability and free-flight aerodynamic performance through a varyingatmosphere at low speed (Mach 0.6-0.2). Surface pressure predicted with FF-CFD agreeswell with experimental trends, with slight over-prediction near the end of the trajectory.The total angle-of-attack for the free-flight portion agrees well with experimental data.
This project primarily sought to address some of the computational cost concerns of detailed simulations by developing improved methods of handling species transport, and chemical kinetics evaluation in a commercial 3D Computational Fluid Dynamics (CFD) environment. Secondary goals were to apply these techniques to fuels and conditions of interest to better understand the key species and reactions required to adequately capture cycle to cycle coupling. Improved modeling will lead to better understanding of these combustion modes and their dependence on fuel composition, which can then enable more clean and efficient engines, ultimately benefiting the consumer as well as the general public. Two approaches were used to address the computational cost. A “Dynamic Species Reduction” (DSR) procedure was developed to remove chemical species from the simulation during periods when chemical reactions were not expected to be important, particularly during gas exchange when temperatures are low and little fuel remains. This modelling procedure automatically detects relevant species to retain in the simulation domain based on their local concentration, removes species below a specified concentration threshold, and then adjusts the remaining species mass to conserve not only the number of H, C, and O atoms in each cell, but also the relative proportions between species and the heating value of the mixture in every cell. The second procedure was a “Product Directed Remapping” (PDR) updates the algorithm used to group individual computational cells for chemical kinetics evaluation to account for non-uniform temperature distributions and the CO to CO 2 ratio in the cells. The model techniques developed in this work successfully demonstrated computational performance improvements for a range of conditions relevant for Highly Dilute SI and HCCI engine operation. Runtime reductions of 10% were observed for small mechanisms, with further reductions of up to 36% observed for larger mechanisms. Improvements were primarily related to reducing the number of chemical species tracked during the gas exchange process using the Dynamic Species Reduction method. With smaller benefits observed from changes to the kinetics binning and evaluation strategy in the post combustion region using the PDR method. The results of this work show the potential for improved computational runtime for complicated simulations. They can and should be extended to additional conditions and new bio-derived and renewable fuels as they are developed and new kinetic mechanisms become available.
The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.
A series of ballistic range tests were performed on a scaled model of the Supersonic Flight Demonstration Test (SFDT) intended to test the Supersonic Inflatable Aerodynamic Decelerator (SIAD) geometry. The purpose of these experiments were to provide aerodynamic coefficients of the vehicle to aid in mission and vehicle design. The experimental data spans the moderate Mach number range, $3.8-2.0$, with a total angle of attack ($alpha_T$) range, $10o-20o$. These conditions are intended to span the Mach-$alpha$ space for the majority of the SFDT experiment. In an effort to validate the predictive capabilities of Computational Fluid Dynamics (CFD) for free-flight aerodynamic behavior, numerical simulations of the ballistic range experiment are performed using the unstructured finite volume Navier-Stokes solver, US3D. Comparisons to raw vehicle attitude, and post-processed aerodynamic coefficients are made between simulated results and experimental data. The resulting comparisons for both raw model attitude and derived aerodynamic coefficients show good agreement with experimental results. Additionally, near body pressure field values for each trajectory simulated are investigated. Extracted surface and wake pressure data gives further insights into dynamic flow coupling leading to a potential mechanism for dynamic instability.
Explore the source record for details and available documents.
Flow visualization, Particle Image Velocimetry (PIV) and dissolved oxygen charging experiments are carried out in a small-scale cylindrical air-water bubble column at different superficial gas velocities. Void fraction profile, total gas holdup, liquid phase velocity field and volumetric mass transfer coefficient are obtained. Details of each experiment are given, with a new procedure for obtaining a 2D void fraction profile from flow visualization by image processing and time averaging which could provide useful information for CFD model validation. The experimental results are analyzed and compared with simulations where good agreement is found for both local and global flow patterns. The averaged governing equations and constitutive relations used in the simulation are discussed briefly, with emphasis on the proper comparison between simulations and experiments based on the derivation of the two-fluid model. Finally, by conducting three sets of experiments related to phase distribution, velocity field and mass transfer in the same apparatus, a validated CFD model for multiphase mass transfer is established and can be used in the future design of a multi-phase reactor.
For long-duration in-space storage of cryogenic propellants, an axial jet mixer is one concept for controlling tank pressure and reducing thermal stratification. Extensive ground-test data from the 1960s to the present exist for tank diameters of 10 ft or less. The design of axial jet mixers for tanks on the order of 30 ft diameter, such as those planned for the Ares V Earth Departure Stage (EDS) LH2 tank, will require scaling of available experimental data from much smaller tanks, as well designing for microgravity effects. This study will assess the ability for Computational Fluid Dynamics (CFD) to handle a change of scale of this magnitude by performing simulations of existing ground-based axial jet mixing experiments at two tank sizes differing by a factor of ten. Simulations of several axial jet configurations for an Ares V scale EDS LH2 tank during low Earth orbit (LEO) coast are evaluated and selected results are also presented. Data from jet mixing experiments performed in the 1960s by General Dynamics with water at two tank sizes (1 and 10 ft diameter) are used to evaluate CFD accuracy. Jet nozzle diameters ranged from 0.032 to 0.25 in. for the 1 ft diameter tank experiments and from 0.625 to 0.875 in. for the 10 ft diameter tank experiments. Thermally stratified layers were created in both tanks prior to turning on the jet mixer. Jet mixer efficiency was determined by monitoring the temperatures on thermocouple rakes in the tanks to time when the stratified layer was mixed out. Dye was frequently injected into the stratified tank and its penetration recorded. There were no velocities or turbulence quantities available in the experimental data. A commercially available, time accurate, multi-dimensional CFD code with free surface tracking (FLOW-3D from Flow Science, Inc.) is used for the simulations presented. Comparisons are made between computed temperatures at various axial locations in the tank at different times and those observed experimentally. The affect of various modeling parameters on the agreement obtained are assessed.
Several wall-modeled large eddy simulation (WMLES) methods are tested by simulating an unsteady Mach 2.0 compression corner geometry in Mississippi State’s Loci/CHEM solver. This study is conducted to evaluate the usage and requirements of these WMLES methods for applications regarding fluctuating pressure environments on launch vehicles with computational fluid dynamics (CFD). Two hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulation (LES) methods, Dynamic Hybrid RANS-LES (DHRL) and Improved Delayed Detached Eddy Simulation (IDDES), and one wall-stress-model, the Algebraic Wall Model for Wall-Modeled LES (AWMLES), are tested on varying grid and timestep refinement levels. These grid and timestep sizes are chosen to test the minimum requirements for successfully running these WMLES methods. The simulations are evaluated based on turbulent boundary layer properties in the developed boundary layer as well as unsteady quantities relating to fluctuating pressure environments in the region of the compression corner. The DHRL method shows good agreement with the comparison wind tunnel data and shows good grid and timestep convergence. The results from the IDDES and WMLES simulations show good agreement for several quantities with some discrepancies regarding others. The results presented in this paper will be used to inform further studies in predicting unsteady environments on higher-complexity geometries.
Explore the source record for details and available documents.
Navier-Stokes computational fluid dynamics calculations are presented for isolated, half-span, and full-span V-22 tiltrotor hover configurations. These computational results extend the validity of CFD hover methodology beyond conventional rotorcraft applications to tiltrotor configurations. Computed steady-state, isolated rotor performance agrees well with experimental measurements, showing little sensitivity to grid resolution. However, blade-vortex interaction flowfield details are sensitive to numerical dissipation and are more difficult to model accurately. Time-dependent, dynamic, half- and full-span installed configurations show sensitivities in performance to the tiltrotor fountain flow. As such, the full-span configuration exhibits higher rotor performance and lower airframe download than the half-span configuration. Half-span rotor installation trends match available half-span data, and airframe downloads are reasonably well predicted. Overall, the CFD solutions provide a wealth of flowfield details that can be used to analyze and improve tiltrotor aerodynamic performance.
Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.
Explore the source record for details and available documents.
It is notoriously difficult to build an accurate universal model for biomass pyrolysis due to its sensitivity to a wide number of critical material attributes such as chemical species and physical sizes. In this work, a biomass pyrolysis kinetics with 32 heterogeneous reactions and 59 species was implemented in an open-source multiphase computational fluid dynamics (CFD) software MFiX and validated against two different experimental pyrolysis data sets that provided detailed data describing chemical component yields. The reaction scheme was then used to build a surrogate model and assess the sensitivity of pyrolysis yields to feedstock compositions. The sensitivity analysis determined that the yield of bio-char showed a strong positive sensitivity to the carbon-rich lignin and tannin pseudo-species in the reaction scheme while the bio-oil and bio-gas were correlated to oxygen-rich lignin pseudo-species. The reaction scheme was then integrated into a coarse-grained discrete element model to simulate fast pyrolysis in a bubbling fluidized bed over a range of feedstock particle sizes. The reactor simulations showed further sensitivity to particle size and hydrodynamics. Notably, particles under 0.5 mm have small heat transfer limitations but left the reactor before completely converting and thus reduced the bio-oil yield. Results from this study can be used to guide future development of highly accurate models for fast pyrolysis reactors with a variety of feedstock properties and operating conditions.
Data from the recent UH-60A Airloads Test in the National Full-Scale Aerodynamics Complex 40- by 80- Foot Wind Tunnel at NASA Ames Research Center are presented and compared to predictions computed by a loosely coupled Computational Fluid Dynamics (CFD)/Comprehensive analysis. Primary calculations model the rotor in free-air, but initial calculations are presented including a model of the tunnel test section. The conditions studied include a speed sweep at constant lift up to an advance ratio of 0.4 and a thrust sweep at constant speed into deep stall. Predictions show reasonable agreement with measurement for integrated performance indicators such as power and propulsive but occasionally deviate significantly. Detailed analysis of sectional airloads reveals good correlation in overall trends for normal force and pitching moment but pitching moment mean often differs. Chord force is frequently plagued by mean shifts and an overprediction of drag on the advancing side. Locations of significant aerodynamic phenomena are predicted accurately although the magnitude of individual events is often missed.