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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Plume-Free Stream Interaction Heating Effects During Orion Crew Module Reentry

During reentry of the Orion Crew Module (CM), vehicle attitude control will be performed by firing reaction control system (RCS) thrusters. Simulation of RCS plumes and their interaction with the oncoming flow has been difficult for the analysis community due to the large scarf angles of the RCS thrusters and the unsteady nature of the Orion capsule backshell environments. The model for the aerothermal database has thus relied on wind tunnel test data to capture the heating effects of thruster plume interactions with the freestream. These data are only valid for the continuum flow regime of the reentry trajectory. A Direct Simulation Monte Carlo (DSMC) analysis was performed to study the vehicle heating effects that result from the RCS thruster plume interaction with the oncoming freestream flow at high altitudes during Orion CM reentry. The study was performed with the DSMC Analysis Code (DAC). The inflow boundary conditions for the jets were obtained from Data Parallel Line Relaxation (DPLR) computational fluid dynamics (CFD) solutions. Simulations were performed for the roll, yaw, pitch-up and pitch-down jets at altitudes of 105 km, 125 km and 160 km as well as vacuum conditions. For comparison purposes (see Figure 1), the freestream conditions were based on previous DAC simulations performed without active RCS to populate the aerodynamic database for the Orion CM. Other inputs to the analysis included a constant Orbital reentry velocity of 7.5 km/s and angle of attack of 160 degrees. The results of the study showed that the interaction effects decrease quickly with increasing altitude. Also, jets with highly scarfed nozzles cause more severe heating compared to the nozzles with lower scarf angles. The difficulty of performing these simulations was based on the maximum number density and the ratio of number densities between the freestream and the plume for each simulation. The lowest altitude solutions required a substantial amount of computational resources (up to 1800 processors) to simulate approximately 2 billion molecules for the refined (adapted) solutions.

Marichalar, J.↗
NASA NTRS2006

Applying Parallel Adaptive Methods with GeoFEST/PYRAMID to Simulate Earth Surface Crustal Dynamics

This viewgraph presentation reviews the use Adaptive Mesh Refinement (AMR) in simulating the Crustal Dynamics of Earth's Surface. AMR simultaneously improves solution quality, time to solution, and computer memory requirements when compared to generating/running on a globally fine mesh. The use of AMR in simulating the dynamics of the Earth's Surface is spurred by future proposed NASA missions, such as InSAR for Earth surface deformation and other measurements. These missions will require support for large-scale adaptive numerical methods using AMR to model observations. AMR was chosen because it has been successful in computation fluid dynamics for predictive simulation of complex flows around complex structures.

GeoFEST↗
NASA NTRS2007

Comparison Between Simulated and Experimentally Measured Performance of a Four Port Wave Rotor

Performance and operability testing has been completed on a laboratory-scale, four-port wave rotor, of the type suitable for use as a topping cycle on a gas turbine engine. Many design aspects, and performance estimates for the wave rotor were determined using a time-accurate, one-dimensional, computational fluid dynamics-based simulation code developed specifically for wave rotors. The code follows a single rotor passage as it moves past the various ports, which in this reference frame become boundary conditions. This paper compares wave rotor performance predicted with the code to that measured during laboratory testing. Both on and off-design operating conditions were examined. Overall, the match between code and rig was found to be quite good. At operating points where there were disparities, the assumption of larger than expected internal leakage rates successfully realigned code predictions and laboratory measurements. Possible mechanisms for such leakage rates are discussed.

Paxson, Daniel E.↗
NASA NTRS2012

AdapChem

AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.

Oluwole, Oluwayemisi O.↗
NASA NTRS2016

Ice-Accretion Test Results for Three Large-Scale Swept-Wing Models in the NASA Icing Research Tunnel

Icing simulation tools and computational fluid dynamics codes are reaching levels of maturity such that they are being proposed by manufacturers for use in certification of aircraft for flight in icing conditions with increasingly less reliance on natural-icing flight testing and icing-wind-tunnel testing. Sufficient high-quality data to evaluate the performance of these tools is not currently available. The objective of this work was to generate a database of ice-accretion geometry that can be used for development and validation of icing simulation tools as well as for aerodynamic testing. Three large-scale swept wing models were built and tested at the NASA Glenn Icing Research Tunnel (IRT). The models represented the Inboard (20 percent semispan), Midspan (64 percent semispan) and Outboard stations (83 percent semispan) of a wing based upon a 65 percent scale version of the Common Research Model (CRM). The IRT models utilized a hybrid design that maintained the full-scale leading-edge geometry with a truncated afterbody and flap. The models were instrumented with surface pressure taps in order to acquire sufficient aerodynamic data to verify the hybrid model design capability to simulate the full-scale wing section. A series of ice-accretion tests were conducted over a range of total temperatures from -23.8 to -1.4 C with all other conditions held constant. The results showed the changing ice-accretion morphology from rime ice at the colder temperatures to highly 3-D scallop ice in the range of -11.2 to -6.3 C. Warmer temperatures generated highly 3-D ice accretion with glaze ice characteristics. The results indicated that the general scallop ice morphology was similar for all three models. Icing results were documented for limited parametric variations in angle of attack, drop size and cloud liquid-water content (LWC). The effect of velocity on ice accretion was documented for the Midspan and Outboard models for a limited number of test cases. The data suggest that there are morphological characteristics of glaze and scallop ice accretion on these swept-wing models that are dependent upon the velocity. This work has resulted in a large database of ice-accretion geometry on large-scale, swept-wing models.

Wing Aerodynamics↗
NASA NTRS2016

Ice-Accretion Test Results for Three Large-Scale Swept-Wing Models in the NASA Icing Research Tunnel

Icing simulation tools and computational fluid dynamics codes are reaching levels of maturity such that they are being proposed by manufacturers for use in certification of aircraft for flight in icing conditions with increasingly less reliance on natural-icing flight testing and icing-wind-tunnel testing. Sufficient high-quality data to evaluate the performance of these tools is not currently available. The objective of this work was to generate a database of ice-accretion geometry that can be used for development and validation of icing simulation tools as well as for aerodynamic testing. Three large-scale swept wing models were built and tested at the NASA Glenn Icing Research Tunnel (IRT). The models represented the Inboard (20% semispan), Midspan (64% semispan) and Outboard stations (83% semispan) of a wing based upon a 65% scale version of the Common Research Model (CRM). The IRT models utilized a hybrid design that maintained the full-scale leading-edge geometry with a truncated afterbody and flap. The models were instrumented with surface pressure taps in order to acquire sufficient aerodynamic data to verify the hybrid model design capability to simulate the full-scale wing section. A series of ice-accretion tests were conducted over a range of total temperatures from -23.8 deg C to -1.4 deg C with all other conditions held constant. The results showed the changing ice-accretion morphology from rime ice at the colder temperatures to highly 3-D scallop ice in the range of -11.2 deg C to -6.3 deg C. Warmer temperatures generated highly 3-D ice accretion with glaze ice characteristics. The results indicated that the general scallop ice morphology was similar for all three models. Icing results were documented for limited parametric variations in angle of attack, drop size and cloud liquid-water content (LWC). The effect of velocity on ice accretion was documented for the Midspan and Outboard models for a limited number of test cases. The data suggest that there are morphological characteristics of glaze and scallop ice accretion on these swept-wing models that are dependent upon the velocity. This work has resulted in a large database of ice-accretion geometry on large-scale, swept-wing models.

Wind-Tunnel Tests↗
NASA NTRS2017

Ablation and Heating During Atmospheric Entry and Its Effect on Airburst Risk

Large meteoroids and asteroids entering the atmosphere endure tremendous heating from the shock heated air, and thereby lose a significant fraction of their mass during atmospheric entry a process known as ablation. The predicted evolution of the asteroids mass as it passes through the atmosphere can affect both the predicted energy deposition profile relevant to an airburst event, or the residual mass that strikes the ground in the case of an impact event. This presentation is divided roughly into two parts. In the first part, an overview of traditional models for heat transfer and ablation that are historically used in the meteor physics community is presented, and the validity in the asteroid entry regime discussed. Sensitivity analyses performed using the recently developed Fragment-Cloud Model (FCM) will be presented which show illustrate the range of sizes and entry parameters for which the predicted asteroid threat is most sensitive to the models for ablation and heat transfer. The second part of the presentation shall focus on recent work done under NASAs Asteroid Threat Assessment Project (ATAP) to develop new models for heat transfer and ablation using high-fidelity numerical simulation in concert with state-of-the-art experiments. Coupled computational fluid dynamics (CFD)radiation transport simulations preformed using the state-of-the-art entry modeling tools at NASA show that, for large meteoroids and asteroids, there can significant attenuation of the heat transfer to the surface (95 in some cases) by the products of ablation. In addition to the heat transfer, new models for the material response and ablation of asteroidal material have been developed [cite]. In the current work, we present finding from recent novel experiments performed in the arc jet facility at NASA Ames, which allows us to, in part, simulate the extreme environment experienced by the asteroid during entry. Briefly, the experimental set-up was comprised of a 1.5 conical article of machined H5 chondrite, which was exposed to a high-enthalpy flow resulting in approximately 4 kWcm2 of heating to the surface. A still frame capture from high-speed video taken during this experiment can be seen in Figure 1. In this figure, we can observe some of the major mechanisms for meteoroid ablation, such as melt flow, spallation (mechanical removal of material), and vaporization. Major findings from this, and other experiments will be discussed, as well progress on utilizing the data from the experiments to inform and develop improved models for ablation.

Asteroids↗
NASA NTRS2019

Generation of a Reduced-Order Model of an Unmanned Combat Air Vehicle Using Indicial Response Functions

Within the NATO STO AVT-251 Task Group, a generic Unmanned Combat Air Vehicle (UCAV) planform was redesigned based on requirements derived from parts of the flight envelope of the defined mission. Because of the lambda-shape planform and associated flow phenomena including shocks and vortices, the aerodynamic design process relies heavily on high-fidelity predictions from computational fluid dynamics (CFD). These simulations enable accurate flight dynamics predictions, allowing for the identification of any potential performance issues early in the design process. Potentially, a complete flight dynamics model can be derived from CFD calculations provided that the underlying reduced order model (ROM) captures the underlying physics contained within the CFD results. This work investigates the creation of a nonlinear ROM using indicial response functions. The response functions are obtained using a grid motion approach that separates the effects of angle of attack and pitch rate. This approach is demonstrated using three different CFD codes from various organizations: ENSOLV at Royal Netherlands Aerospace Centre, USM3D at NASA Langley Research Center, and Kestrel at the United States Air Force Academy. ROM predictions were generated for a manoeuver resembling a low-speed pull-up. The predictions are found to be sensitive to the quality of the steady-state solutions from which step-response calculations were started, as well as the convergence of the step responses themselves. Nevertheless, the indicial response method is shown to provide aerodynamic predictions with acceptable accuracy including transient effects, and as such is a promising method for computationally efficient flight dynamics predictions.

Michel P.C. van Rooij↗
NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗
NASA NTRS2023

Turbulence Model Validation Through Joint Experimental /Computational Studies of Separated Flow Over A Three-Dimensional Tapered Bump: Part I - Experimental Investigation

The results of a three-year experimental campaign aimed at comprehensively documenting the separated flow over a three-dimensional bump are presented with the purpose of generating a benchmark experimental database useful in validating computational fluid dynamics (CFD) flow simulations. The bump model geometry was designed to provide well-defined and repeatable smooth-body flow separation conditions that were suitable for both experiments and simulations. The bump had a Gaussian streamwise profile with a constant height equal to 8.5% of its width over the central 60% of the test section width. The remaining 40% were outboard spanwise portions that gradually taper to zero using an error function profile to minimize side-wall boundary layer effects. The model was immersed in a turbulent boundary layer that was developed on a suspended flat plate in the Notre Dame Mach 0.6 Wind Tunnel. In order to document the effect of the incoming boundary layer thickness on the flow separation, the bump model could be located at two streamwise positions. The mean velocity and turbulence intensity of the wind tunnel freestream flow field and approaching turbulent boundary layer were fully documented. The measurements of the flow separation region included surface visualization, wall shear stress using oil-film interferometry, mean and dynamic surface pressure, and planar and stereoscopic particle image velocimetry. The experiments were conducted over a range of Mach numbers from 0.05 to 0.2 corresponding to a range of Reynolds numbers based on the test section spanwise dimension (0.914m) of 1.0 × 10 6 ≤ Re L = U ∞ L/ν ≤ 4.0 × 10 6 . The bulk of the results are presented for the higher Mach number conditions of 0.1 and 0.2 with Re L = 2.0 × 10 6 and 4.0 × 10 6 , respectively. Extensive uncertainty analysis of the data was performed. The data is archived in the NASA Langley Turbulence Modeling Resource website at https://turbmodels.larc.nasa.gov/Other_exp_Data/speedbump_sep_exp.html. In addition to the experiments, a computational effort was made in parallel by the CFD group at Boeing Research & Technology highlighting the usefulness of the data set, which is outlined in the accompanying CFD report.

Flow separation↗
NASA NTRS2013

A Non-Cut Cell Immersed Boundary Method for Use in Icing Simulations

This paper describes a computational fluid dynamic method used for modelling changes in aircraft geometry due to icing. While an aircraft undergoes icing, the accumulated ice results in a geometric alteration of the aerodynamic surfaces. In computational simulations for icing, it is necessary that the corresponding geometric change is taken into consideration. The method used, herein, for the representation of the geometric change due to icing is a non-cut cell Immersed Boundary Method (IBM). Computational cells that are in a body fitted grid of a clean aerodynamic geometry that are inside a predicted ice formation are identified. An IBM is then used to change these cells from being active computational cells to having properties of viscous solid bodies. This method has been implemented in the NASA developed node centered, finite volume computational fluid dynamics code, FUN3D. The presented capability is tested for two-dimensional airfoils including a clean airfoil, an iced airfoil, and an airfoil in harmonic pitching motion about its quarter chord. For these simulations velocity contours, pressure distributions, coefficients of lift, coefficients of drag, and coefficients of pitching moment about the airfoil's quarter chord are computed and used for comparison against experimental results, a higher order panel method code with viscous effects, XFOIL, and the results from FUN3D's original solution process. The results of the IBM simulations show that the accuracy of the IBM compares satisfactorily with the experimental results, XFOIL results, and the results from FUN3D's original solution process.

Sarofeen, Christian M.↗
NASA NTRS1986

Computational analysis of the SSME fuel preburner flow

A computational fluid dynamics model which simulates the steady state operation of the SSME fuel preburner is developed. Specifically, the model will be used to quantify the flow factors which cause local hot spots in the fuel preburner in order to recommend experiments whereby the control of undesirable flow features can be demonstrated. The results of a two year effort to model the preburner are presented. In this effort, investigating the fuel preburner flowfield, the appropriate transport equations were numerically solved for both an axisymmetric and a three-dimensional configuration. Continuum's VAST (Variational Solution of the Transport equations) code, in conjunction with the CM-1000 Engineering Analysis Workstation and the NASA/Ames CYBER 205, was used to perform the required calculations. It is concluded that the preburner operational anomalies are not due to steady state phenomena and must, therefore, be related to transient operational procedures.

Wang, T. S.↗
NASA NTRS1990

Development of comprehensive numerical schemes for predicting evaporating gas-droplets flow processes of a liquid-fueled combustor

An existing Computational Fluid Dynamics code for simulating complex turbulent flows inside a liquid rocket combustion chamber was validated and further developed. The Advanced Rocket Injector/Combustor Code (ARICC) is simplified and validated against benchmark flow situations for laminar and turbulent flows. The numerical method used in ARICC Code is re-examined for incompressible flow calculations. For turbulent flows, both the subgrid and the two equation k-epsilon turbulence models are studied. Cases tested include idealized Burger's equation in complex geometries and boundaries, a laminar pipe flow, a high Reynolds number turbulent flow, and a confined coaxial jet with recirculations. The accuracy of the algorithm is examined by comparing the numerical results with the analytical solutions as well as experimented data with different grid sizes.

Chen, C. P.↗
NASA NTRS1991

Prediction of jet mean flow structure in support of HSCT noise suppression concepts

The paper describes the application of techniques based on computational fluid dynamics to the simulation of jet flowfields. A solution code for the Reynolds-averaged Navier-Stokes equations is supplemented by conventional two-equation turbulence models based on the Boussinesq approximation. The axisymmetric SCIPVIS code is enhanced with the PARCH and CRAFT codes to examine plug-jet flowfields and imperfectly expanded axisymmetric free round jets. The sensitivity of shock/boundary layer interactions is observed in simulations of the plug case, and the adaptive gridding in the disk region and turbulence levels generated at the triple point are identified as areas in the Mach case that require improvement. Jet-wave structure in the region beyond the first several shock cells can be predicted, and turbulence modeling can be undertaken with respect to improving compressibility, length scale, vorticity, and energy budget. The mean flow structure of imperfectly expanded jets can be studied to develop related noise suppression concepts for the High-Speed Civilian Transport (HSCT).

Sinha, N.↗
NASA NTRS1993

A multiarchitecture parallel-processing development environment

A description is given of the hardware and software of a multiprocessor test bed - the second generation Hypercluster system. The Hypercluster architecture consists of a standard hypercube distributed-memory topology, with multiprocessor shared-memory nodes. By using standard, off-the-shelf hardware, the system can be upgraded to use rapidly improving computer technology. The Hypercluster's multiarchitecture nature makes it suitable for researching parallel algorithms in computational field simulation applications (e.g., computational fluid dynamics). The dedicated test-bed environment of the Hypercluster and its custom-built software allows experiments with various parallel-processing concepts such as message passing algorithms, debugging tools, and computational 'steering'. Such research would be difficult, if not impossible, to achieve on shared, commercial systems.

Townsend, Scott↗
NASA NTRS2000

A 3-D Coupled CFD-DSMC Solution Method With Application to the Mars Sample Return Orbiter

A method to obtain coupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC), 3-D flow field solutions for highly blunt bodies at low incidence is presented and applied to one concept of the Mars Sample Return Orbiter vehicle as a demonstration of the technique. CFD is used to solve the high-density blunt forebody flow defining an inflow boundary condition for a DSMC solution of the afterbody wake flow. By combining the two techniques in flow regions where most applicable, the entire mixed flow field is modeled in an appropriate manner.

Glass, Christopher E.↗
NASA NTRS2000

Numerical Study of Rarefied Hypersonic Flow Interacting with a Continuum Jet

An uncoupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC) technique is developed and applied to provide solutions for continuum jets interacting with rarefied external flows. The technique is based on a correlation of the appropriate Bird breakdown parameter for a transitional-rarefied condition that defines a surface within which the continuum solution is unaffected by the external flow-jet interaction. The method is applied to two problems to assess and demonstrate its validity; one of a jet interaction in the transitional-rarefied flow regime and the other in the moderately rarefied regime. Results show that the appropriate Bird breakdown surface for uncoupling the continuum and non-continuum solutions is a function of a non-dimensional parameter relating the momentum flux and collisionality between the two interacting flows. The correlation is exploited for the simulation of a jet interaction modeled for an experimental condition in the transitional-rarefied flow regime and the validity of the correlation is demonstrated. The uncoupled technique is also applied to an aerobraking flight condition for the Mars Global Surveyor spacecraft with attitude control system jet interaction. Aerodynamic yawing moment coefficients for cases without and with jet interaction at various angles-of-attack were predicted, and results from the present method compare well with values published previously. The flow field and surface properties are analyzed in some detail to describe the mechanism by which the jet interaction affects the aerodynamics.

Glass, Christopher E.↗
NASA NTRS2003

Designing a Hybrid Laminar-Flow Control Experiment: The CFD-Experiment Connection

The NASA/Boeing hybrid laminar flow control (HLFC) experiment, designed during 1993-1994 and conducted in the NASA LaRC 8-foot Transonic Pressure Tunnel in 1995, utilized computational fluid dynamics and numerical simulation of complex fluid mechanics to an unprecedented extent for the design of the test article and measurement equipment. CFD was used in: the design of the test wing, which was carried from definition of desired disturbance growth characteristics, through to the final airfoil shape that would produce those growth characteristics; the design of the suction-surface perforation pattern that produced enhanced crossflow-disturbance growth: and in the design of the hot-wire traverse system that produced minimal influence on measured disturbance growth. These and other aspects of the design of the test are discussed, after the historical and technical context of the experiment is described.

Streett, C. L.↗