Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Fluid Dynamic (CFD)”

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 973 records · Page 54

Assessment of the Unstructured Grid Software TetrUSS for Drag Prediction of the DLR-F4 Configuration

An application of the NASA unstructured grid software system TetrUSS is presented for the prediction of aerodynamic drag on a transport configuration. The paper briefly describes the underlying methodology and summarizes the results obtained on the DLR-F4 transport configuration recently presented in the first AIAA computational fluid dynamics (CFD) Drag Prediction Workshop. TetrUSS is a suite of loosely coupled unstructured grid CFD codes developed at the NASA Langley Research Center. The meshing approach is based on the advancing-front and the advancing-layers procedures. The flow solver employs a cell-centered, finite volume scheme for solving the Reynolds Averaged Navier-Stokes equations on tetrahedral grids. For the present computations, flow in the viscous sublayer has been modeled with an analytical wall function. The emphasis of the paper is placed on the practicality of the methodology for accurately predicting aerodynamic drag data.

Pirzadeh, Shahyar Z.↗

Large-Eddy Simulation Code Developed for Propulsion Applications

A large-eddy simulation (LES) code was developed at the NASA Glenn Research Center to provide more accurate and detailed computational analyses of propulsion flow fields. The accuracy of current computational fluid dynamics (CFD) methods is limited primarily by their inability to properly account for the turbulent motion present in virtually all propulsion flows. Because the efficiency and performance of a propulsion system are highly dependent on the details of this turbulent motion, it is critical for CFD to accurately model it. The LES code promises to give new CFD simulations an advantage over older methods by directly computing the large turbulent eddies, to correctly predict their effect on a propulsion system. Turbulent motion is a random, unsteady process whose behavior is difficult to predict through computer simulations. Current methods are based on Reynolds-Averaged Navier- Stokes (RANS) analyses that rely on models to represent the effect of turbulence within a flow field. The quality of the results depends on the quality of the model and its applicability to the type of flow field being studied. LES promises to be more accurate because it drastically reduces the amount of modeling necessary. It is the logical step toward improving turbulent flow predictions. In LES, the large-scale dominant turbulent motion is computed directly, leaving only the less significant small turbulent scales to be modeled. As part of the prediction, the LES method generates detailed information on the turbulence itself, providing important information for other applications, such as aeroacoustics. The LES code developed at Glenn for propulsion flow fields is being used to both analyze propulsion system components and test improved LES algorithms (subgrid-scale models, filters, and numerical schemes). The code solves the compressible Favre-filtered Navier- Stokes equations using an explicit fourth-order accurate numerical scheme, it incorporates a compressible form of Smagorinsky s model for the subgrid-scale turbulence, and it uses generalized curvilinear coordinates to allow analysis of a wide range of geometries. The code runs in parallel on shared memory multiprocessor computers and is written in Fortran 90 with dynamic memory allocation. A sample result for a Mach-1.4 round jet is presented in the figure. Instantaneous Mach number contours in several cross-planes downstream of the nozzle exit are shown, illustrating how an LES captures the large unsteady three-dimensional turbulent structures present in the jet.

DeBonis, James R.↗

Computational Fluid Dynamics Simulation of Dual Bell Nozzle Film Cooling

Marshall Space Flight Center (MSFC) and Armstrong Flight Research Center (AFRC) are working together to advance the technology readiness level (TRL) of the dual bell nozzle concept. Dual bell nozzles are a form of altitude compensating nozzle that consists of two connecting bell contours. At low altitude the nozzle flows fully in the first, relatively lower area ratio, nozzle. The nozzle flow separates from the wall at the inflection point which joins the two bell contours. This relatively low expansion results in higher nozzle efficiency during the low altitude portion of the launch. As ambient pressure decreases with increasing altitude, the nozzle flow will expand to fill the relatively large area ratio second nozzle. The larger area ratio of the second bell enables higher Isp during the high altitude and vacuum portions of the launch. Despite a long history of theoretical consideration and promise towards improving rocket performance, dual bell nozzles have yet to be developed for practical use and have seen only limited testing. One barrier to use of dual bell nozzles is the lack of control over the nozzle flow transition from the first bell to the second bell during operation. A method that this team is pursuing to enhance the controllability of the nozzle flow transition is manipulation of the film coolant that is injected near the inflection between the two bell contours. Computational fluid dynamics (CFD) analysis is being run to assess the degree of control over nozzle flow transition generated via manipulation of the film injection. A cold flow dual bell nozzle, without film coolant, was tested over a range of simulated altitudes in 2004 in MSFC's nozzle test facility. Both NASA centers have performed a series of simulations of that dual bell to validate their computational models. Those CFD results are compared to the experimental results within this paper. MSFC then proceeded to add film injection to the CFD grid of the dual bell nozzle. A series of nozzle pressure ratios and film coolant flow rates are investigated to determine the effect of the film injection on the nozzle flow transition behavior. The results of this CFD study of a dual bell with film injection are presented in this paper.

Braman, Kalen↗

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗

Comparisons of a Three-Dimensional, Full Navier Stokes Computer Model with High Mach Number Combuster Test Data

Comparisons between scramjet combustor data and a three-dimensional full Navier-Stokes calculation have been made to verify and substantiate computational fluid dynamics (CFD) codes and application procedures. High Mach number scramjet combustor development will rely heavily on CFD applications to provide wind tunnel-equivalent data of quality sufficient to design, build and fly hypersonic aircraft. Therefore. detailed comparisons between CFD results and test data are imperative. An experimental case is presented, for which combustor wall static pressures were measured and flow-fieid interferograms were obtained. A computer model was done of the experiment, and counterpart parameters are compared with experiment. The experiment involved a subscale combustor designed and fabricated for the National Aero-Space Plane Program, and tested in the Calspan Corporation 96" hypersonic shock tunnel. The combustor inlet ramp was inclined at a 20 angle to the shock tunnel nozzle axis, and resulting combustor entrance flow conditions simulated freestream M=10. The combustor body and cowl walls were instrumented with static pressure transducers, and the combustor lateral walls contained windows through which flowfield holographic interferograms were obtained. The CFD calculation involved a three-dimensional time-averaged full Navier-Stokes code applied to the axial flow segment containing fuel injection and combustion. The full Navier-Stokes approach allowed for mixed supersonic and subsonic flow, downstream-upstream communication in subsonic flow regions, and effects of adverse pressure gradients. The code included hydrogen-air chemistry in the combustor segment which begins near fuel injection and continues through combustor exhaust. Combustor ramp and inlet segments on the combustor lateral centerline were modelled as two dimensional. Comparisons to be shown include calculated versus measured wall static pressures as functions of axial flow coordinate, and calculated path-averaged density contours versus an holographic Interferogram.

Watkins, William B.↗

An investigation of turbulence models

The accuracy to which a turbulent boundary layer or wake can be predicted numerically depends on the validity of the turbulence closure model used. The modeling of turbulence physics is one of the most difficult problems in computational fluid dynamics (CFD). In fact, it is one of the pacing factors in the development of CFD. In general, there are three main approaches to the description of trubulence physics. First is turbulence modeling in which the Reynolds averaged Navier-Stokes equations are used and some closure approximation is made for the Reynolds stresses. A second approach to turbulence is large eddy simulation (LES) in which the computational mesh is taken to be fine enough that the large scale structure of the turbulence can be calculated directly. An empirical assumption must be made for the small scale sub-grid turbulence. The third approach is direct simulation. In this technique the Navier-Strokes equations are solved directly on a mesh which if fine enough to resolve the smallest length scale of the turbulence. The Reynolds averaged equations are not used and no closure assumption is required. These last two approaches require extensive computer resources and as such are not engineering tools. The purpose of the work was to investigate the various engineering turbulence models for accuracy and ease of programming. This involved comparison of the models with each other and with experimental data.

Doria, Michael L.↗

Applied transonics at Grumman

A review of several applications of Computational Fluid Dynamics (CFD) to various aspects of aerodynamic design recently carried out at Grumman is presented. The emphasis is placed on project-oriented applications where the ease of use of the methods and short start-to-completion times are required. Applications cover transonic wing design/optimization, wing mounted stores load prediction, transonic buffet alleviation, fuselage loads estimation, and compact offset diffuser design for advanced aircraft configurations. Computational methods employed include extended transonic small disturbance (automatic grid embedding) formulation for analysis/design/optimization and a thin layer Navier-Stokes formulation for both external and internal flow analyses.

Davis, W. H.↗

SINDA/FLUINT and Thermal Desktop Multi-Node Settled and Unsettled Propellant Tank Modeling of Zero Boil Off Test

Cryogenic propellant storage tank self-pressurization involves complex physical phenomena which are usually analytically modelled via complex multidimensional CFD (Computational Fluid Dynamics) codes. Unfortunately these codes, even when modelling axisymmetric domains, may take weeks or longer to obtain transient pressure and temperature information for relatively short periods of time (several seconds to several hours). Propellant tank storage end-to-end mission simulations can last a duration of days to weeks to months. Multi-node modelling of propellant tanks is a viable alternative to traditional CFD modelling and presents the advantage of greatly reduced run times on the order of hours and days compared to the weeks or longer for CFD codes. A multi-node model represents the fluid within the storage tank, as well as the storage tank itself, as a fluid-thermal network. This type of setup is not necessarily geometrically based. This can be accomplished using a commercial generalized fluid-thermal network code, such as SINDA/FLUINT (SF). The advantage of using a fluid-thermal network code like SF lies in its extensive ability to model the external environment of the storage tank through the graphical user interface, Thermal Desktop (TD). The total heat load into the tank may be a function of heaters and a complex radiative environment as well. Thermal Desktop may be used to address the detailed radiative environment of the tank as well as building a geometrically accurate depiction of the storage tank itself.

Thermal Desktop↗

Nearfield Summary and Analysis of the Third AIAA Sonic Boom Prediction Workshop C608 Low Boom Demonstrator

A summary and statistical analysis of the nearfield Computational Fluid Dynamics (CFD)submissions for the Third AIAA Sonic Boom Prediction Workshop is provided with a focus onthe C608 Low Boom Flight Test Demonstrator. The C608 is more complex in terms of geometryand propulsion boundary conditions than previous workshop cases and is more representativeof vehicles with lower ground loudness and the potential for lower annoyance. The nearfieldsignatures submitted by the participants are propagated to the ground to compute statisticsof loudness measures over the vehicle sonic boom carpet. Principle component analysis isused to extract the primary variation modes. Context from previous sonic boom workshopsindicates that this workshop has the lowest variation even though the case is more challengingbecause it is quieter and more complex. The results documented in this summary indicate thatthe international state-of-the-art for nearfield CFD has a variation that is small enough formeaningful low-boom design and analysis. Low-boom configuration analysis methods with lowvariation are important tools to develop certification processes for addressing the prohibitionon overland supersonic commercial flight.

sonic boom CFD workshop c608↗

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Microgravity↗

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Crystal Growth↗

Machine Learning-Based Regression Models for Ironmaking Blast Furnace Automation

Computational fluid dynamics (CFD)-based simulation has been the traditional way to model complex industrial systems and processes. One very large and complex industrial system that has benefited from CFD-based simulations is the steel blast furnace system. The problem with the CFD-based simulation approach is that it tends to be very slow for generating data. The CFD-only approach may not be fast enough for use in real-time decisionmaking. To address this issue, in this work, the authors propose the use of machine learning techniques to train and test models based on data generated via CFD simulation. Regression models based on neural networks are compared with tree-boosting models. In particular, several areas (tuyere, raceway, and shaft) of the blast furnace are modeled using these approaches. The results of the model training and testing are presented and discussed. The obtained R 2 metrics are, in general, very high. The results appear promising and may help to improve the efficiency of operator and process engineer decisionmaking when running a blast furnace.

97 MATHEMATICS AND COMPUTING↗

Computational aerodynamics - The next generation

Continued advances in the various elements that comprise the field of computational fluid dynamics (CFD) are promoting a radically different approach to the aerodynamic design and analysis of aerospace vehicles and systems. The elements of CFD generally include numerical algorithm development, transition and turbulence modeling, surface modeling, and grid generation, scientific visualization and validation methodologies. This paper discusses the research progress and prospects for the future in each of these elements within NASA's CFD and Experimental Validation Program. The applicability of computational methods for the purposes of understanding complex flow phenomena, exploring aerodynamic concepts, and providing vehicle-design input is also addressed.

Hessenius, Kristin A.↗

Experimental and Computational Analysis of Unidirectional Flow Through Stirling Engine Heater Head

A high efficiency Stirling Radioisotope Generator (SRG) is being developed for possible use in long-duration space science missions. NASA s advanced technology goals for next generation Stirling convertors include increasing the Carnot efficiency and percent of Carnot efficiency. To help achieve these goals, a multi-dimensional Computational Fluid Dynamics (CFD) code is being developed to numerically model unsteady fluid flow and heat transfer phenomena of the oscillating working gas inside Stirling convertors. In the absence of transient pressure drop data for the zero mean oscillating multi-dimensional flows present in the Technology Demonstration Convertors on test at NASA Glenn Research Center, unidirectional flow pressure drop test data is used to compare against 2D and 3D computational solutions. This study focuses on tracking pressure drop and mass flow rate data for unidirectional flow though a Stirling heater head using a commercial CFD code (CFD-ACE). The commercial CFD code uses a porous-media model which is dependent on permeability and the inertial coefficient present in the linear and nonlinear terms of the Darcy-Forchheimer equation. Permeability and inertial coefficient were calculated from unidirectional flow test data. CFD simulations of the unidirectional flow test were validated using the porous-media model input parameters which increased simulation accuracy by 14 percent on average.

Wilson, Scott D.↗

Temporally continuous thermofluidic–thermomechanical modeling framework for metal additive manufacturing

Additive manufacturing (AM) is known to generate large magnitudes of residual stresses (RS) within builds due to steep and localized thermal gradients. In the current state of commercial AM technology, manufacturers generally perform heat treatments in effort to reduce the generated RS and its detrimental effects on part distortion and in-service failure. Computational models that effectively simulate the deposition process can provide valuable insights to improve RS distributions. Accordingly, it is common to employ Computational fluid dynamics (CFD) models or finite element (FE) models. While CFD can predict geometric and thermal-fluid behavior, it cannot predict the structural response (e.g., stress–strain) behavior. On the other hand, an FE model can predict mechanical behavior, but it lacks the ability to predict geometric and fluid behavior. Thus, an effectively integrated thermofluidic–thermomechanical modeling framework that exploits the benefits of both techniques while avoiding their respective limitations can offer valuable predictive capability for AM processes. In contrast to previously published efforts, the work herein describes a one-way coupled CFD-FEA framework that abandons major simplifying assumptions, such as geometric steady-state conditions, the absence of material plasticity, and the lack of detailed RS evolution/accumulation during deposition, as well as insufficient validation of results. Here, the presented framework is demonstrated for a directed energy deposition (DED) process, and experiments are performed to validate the predicted geometry and RS profile. Both single- and double-layer stainless steel 316L builds are considered. Geometric data is acquired via 3D optical surface scans and X-ray micro-computed tomography, and residual stress is measured using neutron diffraction (ND). Comparisons between the simulations and measurements reveal that the described CFD-FEA framework is effective in capturing the coupled thermomechanical and thermofluidic behaviors of the DED process. The methodology presented is extensible to other metal AM processes, including power bed fusion and wire-feed-based AM.

42 ENGINEERING↗