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,243 records · Page 69

Combined Uncertainty and A-Posteriori Error Bound Estimates for CFD Calculations: Theory and Implementation

Simulation codes often utilize finite-dimensional approximation resulting in numerical error. Some examples include, numerical methods utilizing grids and finite-dimensional basis functions, particle methods using a finite number of particles. These same simulation codes also often contain sources of uncertainty, for example, uncertain parameters and fields associated with the imposition of initial and boundary data,uncertain physical model parameters such as chemical reaction rates, mixture model parameters, material property parameters, etc.

CFD↗

Retooling CFD for hypersonic aircraft

The CFD facility requirements of hypersonic aircraft configuration design development are different from those thus far employed for reentry vehicle design, because (1) the airframe and the propulsion system must be fully integrated to achieve the desired performance; (2) the vehicle must be reusable, with minimum refurbishment requirements between flights; and (3) vehicle performance must be optimized for a wide range of Mach numbers. An evaluation is presently made of flow resolution within shock waves, transition and turbulence phenomenon tractability, chemical reaction modeling, and hypersonic boundary layer transition, with state-of-the-art CFD.

Dwoyer, Douglas L.↗

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗

Microbial Communities in Microgravity: Simulation in Lab and on the Computer

Microorganisms grow differently in spaceflight than they do on Earth. While much remains unexplained about how microgravity affects microbial growth, one dominant hypothesis is that the lack of density-driven convection in the liquid growth environment makes mixing diffusion-limited, and therefore slower. This is supported by evidence that individual microbial strains experience starvation and acid stress in microgravity. However, if it is true, then microgravity would also have measurable effects on microbes in mixed communities, because many interspecies interactions involve the exchange of soluble metabolites through the medium (cross-feeding). Specifically, cooperative cross-feeding communities would grow more slowly in microgravity, and cooperation would be less stable on evolutionary timescales. Here we describe our efforts to test this hypothesis by simulating microgravity in silico and in the lab, using a model system of Eschericia coli and Salmonella enterica that grow only when they can exchange methionine and acetate. We created CAMDLES (CFD-DEM Artificial Microgravity Developments for Living Ecosystem Simulation) as an extension of CFDEM®coupling software, to carry out computational modeling of biological flows, growth, and mass transfer in microgravity and also in laboratory artificial microgravity devices (rotating wall vessels, RWV). Using CAMDLES, we found distinct differences in growth rates between RWV and true microgravity, and we were able to identify several features, such as spatial distribution, biofilm formation, and product yield parameters, that influence the degree to which RWV growth recapitulates microgravity growth. In addition, we report on the development of a laboratory system for monitoring growth rates and species ratios of the community in RWVs, using fluorescent strains. Pairing CAMDLES with the laboratory model system allows us to generate quantitative predictions about the effects of spaceflight on organisms that will be essential to sustaining human space exploration in the long term.

microbiology↗

Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps

NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.

Deployable Entry Vehicle↗

Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps

NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.

Deployable Entry Vehicle↗

High speed civil transport: Sonic boom softening and aerodynamic optimization

An improvement in sonic boom extrapolation techniques has been the desire of aerospace designers for years. This is because the linear acoustic theory developed in the 60's is incapable of predicting the nonlinear phenomenon of shock wave propagation. On the other hand, CFD techniques are too computationally expensive to employ on sonic boom problems. Therefore, this research focused on the development of a fast and accurate sonic boom extrapolation method that solves the Euler equations for axisymmetric flow. This new technique has brought the sonic boom extrapolation techniques up to the standards of the 90's. Parallel computing is a fast growing subject in the field of computer science because of its promising speed. A new optimizer (IIOWA) for the parallel computing environment has been developed and tested for aerodynamic drag minimization. This is a promising method for CFD optimization making use of the computational resources of workstations, which unlike supercomputers can spend most of their time idle. Finally, the OAW concept is attractive because of its overall theoretical performance. In order to fully understand the concept, a wind-tunnel model was built and is currently being tested at NASA Ames Research Center. The CFD calculations performed under this cooperative agreement helped to identify the problem of the flow separation, and also aided the design by optimizing the wing deflection for roll trim.

Cheung, Samson↗

Cryogenic Extension of NASA Species Polynomials Using Hydrogen and Oxygen at Stoichiometry

NASA has been conducting research in Rotating Detonation Rocket Engines (RDREs) for several years. A recent test at the Marshall Spaceflight Center has successfully shown operation of a RDRE that ran for 251 seconds which had combustion chamber cryogenic inlet conditions [1]. However, an overall pressure gain was not observed in the performance of the engine. Such an indication could mean an optimal chamber design is yet to be discovered. Such a design could significantly reduce parasitic pressure loses and improve overall engine efficiency beyond that of conventional rockets. To investigate different chamber configurations, new CFD capability is needed to accurately capture cryogenic thermophysical and transport property data. Such an effort is currently part of NASA’s Early Career Initiative (ECI) program. This paper focuses on the thermophysical property mixture model approach that makes use of NASA’s polynomial fits of species and another cryogenic model chosen from NIST’s Refprop program. Presented in this paper are results of a 1D denotation CFD simulation of stoichiometric Hydrogen and Oxygen mixture at a cryogenic upstream condition.

Combustion↗

Progress Towards the CFD Vision 2030

Outline: NASA Motivation; CFD Vision 2030 Study Summary - The Vision - Roadmap - Validation Experiments - CFD 2030 Activities Since Report Release; T3 Technical Challenges Supporting CFD Vision 2030 - Revolutionary Computational Aerosciences (RCA) [Aerodynamics - Aerothermodynamics - Aerostructures - Aeroacoustics - Propulsion Integration] - Multidisciplinary Design, Analysis, and Optimization (MDAO) - Combustion Modeling - Innovative Measurements (Proposed); Implementation Plan, Going Forward; Summary.

Vision 2030↗

Comparison of Computational Predictions of the Mach 0.80 Transonic Truss-Braced Wing Configuration with Experimental Data

The NASA Langley and Ames Research Centers have teamed together for comparisons of computational predictions of the Boeing Mach 0.80Transonic Truss-Braced Wing (TTBW) configuration with a high-speed experimental dataset. The Mach 0.80 TTBW vehicle is a high wing, high aspect ratio configuration, designed for a high lift-to-drag ratio. System studies have predicted significant fuel burn and emissions benefits with the TTBW technology moving toward meeting NASA Subsonic Transport Systems-Level-Metrics. A 4.5% scale Mach 0.80 design TTBW model was recently tested at the NASA Ames Research Center 11-by 11-Foot Transonic Wind Tunnel(11-Ft TWT),providing a valuable dataset to validate computational tools and investigate best practices as risk reduction efforts continue for the development of the advanced TTBW vehicle. The NASA Computational Fluid Dynamics (CFD) team has computed free-air flow solutions on the Mach 0.80 design flight configuration and two wind tunnel configuration variants using the USM3D and LAVA flow solvers. Accurate modeling of the configuration tested in the wind tunnel environment is critical to validating the CFD tools, thus the team has included the internal cavity region and sting in their modeling of the configuration, similar to that tested in the 11-Ft TWT. Overall, the CFD simulations compared well and show similar trends as the corrected experimental data for lift and drag polars. The CFD predicted lift curve is shifted in angle of attack from what was observed in the experiment. The shift in lift also was seen in the pitching moment comparison plots. CFD solutions were computed at constant CL test point values and showed overall very good agreement when comparing constant spanwise cuts of pressure coefficient data on the wing and strut with experimental data. CFD cavity corrections were also investigated using the 11-Ft TWT cavity correction method, similar to that used to correct the wind tunnel data. Results showed some improvement in pitching moment coefficient predictions, and an increase in drag, shifting the data to the right in drag polars, further from the experimental data at lower lift conditions, good agreement near the design CL, and slight improvement at the higher lift conditions.

Transonic Truss-Braced Wing↗

Development and Validation of a Slosh-Based Ullage Collapse Model

Ullage collapse modeling is an essential tool for the success of spacecraft powered by cryogenic propellants. The degree of ullage collapse can have a major impact on the design of cryogenic propellant tanks as well as the mission’s concept of operations. Additional commodities (e.g., helium gas for pressure control) and their quantities are dependent on the amount of heat transfer and phase change of propellants expected during flight. Currently, high-fidelity computational fluid dynamics (CFD) with heat transfer and phase change is the most accurate method for predicting ullage collapse. CFD simulations of this kind, however, can become very costly for problems requiring large domains, long simulation times, or significant liquid motion. Consequently, there is a need for a reduced order model to provide fast responses which can be used in engineering design decisions. The proposed numerical model is a nodal approach, augmented by outputs from CFD analyses that are generated rapidly and reliably. These outputs include transient values for liquid-vapor interface area, interface velocity, and liquid position using a volume of fluid (VoF) method. There are several benefits of combining a nodal approach with slosh predictions from CFD simulations: leveraging of the Loci/STREAM-VoF solver that has been validated in high-g and low-g environments, rapid modeling of VoF with isothermal propellant simplifications and relatively coarse meshes, and improved model predictions from the nodal tool. The nodal approach, informed by the slosh predictions from CFD, assumes uniform conditions in the liquid and gas domains. Modifications to account for temperature stratification in the liquid and gas may be implemented as a future improvement. Thermodynamic changes within the tank are solved at user-defined timesteps. First, impacts due to mass entering or leaving the control volume are captured. Second, heat transfer impacts are accounted for by tracking the volume of the gas which has been contacted by liquid in that timestep and rebalancing the internal energy. Finally, phase change is approximated using the area and velocity of the interface. The tank conditions are re-computed by conserving mass and energy. Calibration of the model to several test flights and experiments are required to baseline the predictive power and uncertainty. The final product will enable designers and engineers to create quick-turnaround predictions for ullage collapse and inform design feasibility.

Computational Fluid Dynamics↗

Rotor Airloads Prediction Using Unstructured Meshes and Loose CFD/CSD Coupling

The FUN3D unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids has been modified to allow prediction of trimmed rotorcraft airloads. The trim of the rotorcraft and the aeroelastic deformation of the rotor blades are accounted for via loose coupling with the CAMRAD II rotorcraft computational structural dynamics code. The set of codes is used to analyze the HART-II Baseline, Minimum Noise and Minimum Vibration test conditions. The loose coupling approach is found to be stable and convergent for the cases considered. Comparison of the resulting airloads and structural deformations with experimentally measured data is presented. The effect of grid resolution and temporal accuracy is examined. Rotorcraft airloads prediction presents a very substantial challenge for Computational Fluid Dynamics (CFD). Not only must the unsteady nature of the flow be accurately modeled, but since most rotorcraft blades are not structurally stiff, an accurate simulation must account for the blade structural dynamics. In addition, trim of the rotorcraft to desired thrust and moment targets depends on both aerodynamic loads and structural deformation, and vice versa. Further, interaction of the fuselage with the rotor flow field can be important, so that relative motion between the blades and the fuselage must be accommodated. Thus a complete simulation requires coupled aerodynamics, structures and trim, with the ability to model geometrically complex configurations. NASA has recently initiated a Subsonic Rotary Wing (SRW) Project under the overall Fundamental Aeronautics Program. Within the context of SRW are efforts aimed at furthering the state of the art of high-fidelity rotorcraft flow simulations, using both structured and unstructured meshes. Structured-mesh solvers have an advantage in computation speed, but even though remarkably complex configurations may be accommodated using the overset grid approach, generation of complex structured-mesh systems can require months to set up. As a result, many rotorcraft simulations using structured-grid CFD neglect the fuselage. On the other hand, unstructured-mesh solvers are easily able to handle complex geometries, but suffer from slower execution speed. However, advances in both computer hardware and CFD algorithms have made previously state-of-the-art computations routine for unstructured-mesh solvers, so that rotorcraft simulations using unstructured grids are now viable. The aim of the present work is to develop a first principles rotorcraft simulation tool based on an unstructured CFD solver.

Biedron, Robert T.↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗

Three-Dimensional Effects in Multi-Element High Lift Computations

In an effort to discover the causes for disagreement between previous two-dimensional (2-D) computations and nominally 2-D experiment for flow over the three-element McDonnell Douglas 30P-30N airfoil configuration at high lift, a combined experimental/CFD investigation is described. The experiment explores several different side-wall boundary layer control venting patterns, documents venting mass flow rates, and looks at corner surface flow patterns. The experimental angle of attack at maximum lift is found to be sensitive to the side-wall venting pattern: a particular pattern increases the angle of attack at maximum lift by at least 2 deg. A significant amount of spanwise pressure variation is present at angles of attack near maximum lift. A CFD study using three-dimensional (3-D) structured-grid computations, which includes the modeling of side-wall venting, is employed to investigate 3-D effects on the flow. Side-wall suction strength is found to affect the angle at which maximum lift is predicted. Maximum lift in the CFD is shown to be limited by the growth of an off-body corner flow vortex and consequent increase in spanwise pressure variation and decrease in circulation. The 3-D computations with and without wall venting predict similar trends to experiment at low angles of attack, but either stall too early or else overpredict lift levels near maximum lift by as much as 5%. Unstructured-grid computations demonstrate that mounting brackets lower the lift levels near maximum lift conditions.

Rumsey, Christopher L.↗

Turbine Internal and Film Cooling Modeling For 3D Navier-Stokes Codes

The aim of this research project is to make use of NASA Glenn on-site computational facilities in order to develop, validate and apply aerodynamic, heat transfer, and turbine cooling models for use in advanced 3D Navier-Stokes Computational Fluid Dynamics (CFD) codes such as the Glenn-" code. Specific areas of effort include: Application of the Glenn-HT code to specific configurations made available under Turbine Based Combined Cycle (TBCC), and Ultra Efficient Engine Technology (UEET) projects. Validating the use of a multi-block code for the time accurate computation of the detailed flow and heat transfer of cooled turbine airfoils. The goal of the current research is to improve the predictive ability of the Glenn-HT code. This will enable one to design more efficient turbine components for both aviation and power generation. The models will be tested against specific configurations provided by NASA Glenn.

DeWitt, Kenneth↗

A Computational Investigation of Gear Windage

A CFD method has been developed for application to gear windage aerodynamics. The goals of this research are to develop and validate numerical and modeling approaches for these systems, to develop physical understanding of the aerodynamics of gear windage loss, including the physics of loss mitigation strategies, and to propose and evaluate new approaches for minimizing loss. Absolute and relative frame CFD simulation, overset gridding, multiphase flow analysis, and sub-layer resolved turbulence modeling were brought to bear in achieving these goals. Several spur gear geometries were studied for which experimental data are available. Various shrouding configurations and free-spinning (no shroud) cases were studied. Comparisons are made with experimental data from the open literature, and data recently obtained in the NASA Glenn Research Center Gear Windage Test Facility. The results show good agreement with experiment. Interrogation of the validative and exploratory CFD results have led, for the first time, to a detailed understanding of the physical mechanisms of gear windage loss, and have led to newly proposed mitigation strategies whose effectiveness is computationally explored.

Hill, Matthew J.↗

Overview and Summary of the Second AIAA High Lift Prediction Workshop

The second AIAA CFD High-Lift Prediction Workshop was held in San Diego, California, in June 2013. The goals of the workshop continued in the tradition of the first high-lift workshop: to assess the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in landing/takeoff (high-lift) configurations. This workshop analyzed the flow over the DLR-F11 model in landing configuration at two different Reynolds numbers. Twenty-six participants submitted a total of 48 data sets of CFD results. A variety of grid systems (both structured and unstructured) were used. Trends due to grid density and Reynolds number were analyzed, and effects of support brackets were also included. This paper analyzes the combined results from all workshop participants. Comparisons with experimental data are made. A statistical summary of the CFD results is also included.

Rumsey, Christopher L.↗