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,531 records · Page 85

CFD Simulations of Landing and Takeoff CRM High-Lift Configurations

Numerical simulations have been performed for a conventional high-lift version of the Common Research Model (CRM) model corresponding to landing and take- off configurations. Computed values of lift and drag for the landing configuration are compared with the experimental data acquired in the 14- by 22-Foot Subsonic Tunnel (14×22) at the NASA Langley Research Center (LaRC). Simulations replicated the experimentally observed improvements in the lift characteristics in the presence of a nacelle chine. Flow visualization images indicate that vortices generated by the nacelle chine reduce flow separation regions on the upper surface of the wing at higher angles of attack. Based on such observations, the take-off configuration considered is also equipped with a nacelle chine, and testing of this configuration is planned in the 14×22 tunnel in near future. Preliminary solutions are also presented to explore the feasibility of using a localized flap gap blowing (LFGB) active flow control concept for improving the aerodynamic performance at take-off conditions.

aerodynamics↗

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

Surface definition and grid generation about an Assured Crew Return Vehicle (ACRV) for Space Station Freedom

The surface definition and grid generation about an Assured Crew Return Vehicle (ACRV) for the Space Station Freedom are described. The purpose of the surface definition and grid generation is to provide the necessary geometry information for CFD calculations about the vehicle. There are two salient features in this description. The first is that the numerical model representing the ACRV configuration is obtained from the measurement of an existing wind tunnel model. The method for smoothing the measured data and obtaining the numerical model is described. The second feature is the description of the algebraic grid generation method and software to compute volume grids about the ACRV. The methods and software allow rapid computation of volume grids for a wide range of flow conditions.

Smith, R. E.↗

Validation of CFD for aircraft design: A Boeing view

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

Vas, Irwin E.↗

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.

Venkateswaran,S.↗

Hybrid MPI+OpenMP Programming of an Overset CFD Solver and Performance Investigations

This report describes a two level parallelization of a Computational Fluid Dynamic (CFD) solver with multi-zone overset structured grids. The approach is based on a hybrid MPI+OpenMP programming model suitable for shared memory and clusters of shared memory machines. The performance investigations of the hybrid application on an SGI Origin2000 (O2K) machine is reported using medium and large scale test problems.

Djomehri, M. Jahed↗

Wind Tunnel Model Design for the Study of Plume Effects on Sonic Boom for Isolated Exhaust Nozzles

A low cost test capability was developed at the NASA Glenn Research Center 1- by 1-Foot Supersonic Wind Tunnel (SWT), with a goal to reduce the disturbance caused by supersonic aircraft flight over populated areas. This work focused on the shock wave structure caused by the exhaust nozzle plume. Analysis and design was performed on a new rig to test exhaust nozzle plume effects on sonic boom signature. Test capability included a baseline nozzle test article and a wind tunnel model consisting of a strut, a nose cone and an upper plenum. Analysis was performed on the external and internal aerodynamic configuration, including the shock reflections from the wind tunnel walls caused by the presence of the model nosecone. This wind tunnel model was designed to operate from Mach 1.4 to Mach 3.0 with nozzle pressure ratios from 6 to 12 and altitudes from 30,000 ft (4.36 psia) to 50,000 ft (1.68 psia). The model design was based on a 1 in. outer diameter, was 9 in. in overall length, and was mounted in the wind tunnel on a 3/8 in. wide support strut. For test conditions at 50,000 ft the strut was built to supply 90 psia of pressure, and to achieve 20 psia at the nozzle inlet with a maximum nozzle pressure of 52 psia. Instrumentation was developed to measure nozzle pressure ratio, and an external static pressure probe was designed to survey near field static pressure profiles at one nozzle diameter above the rig centerline. Model layout placed test nozzles between two transparent sidewalls in the 1x1 SWT for Schlieren photography and comparison to CFD analysis.

Castner, Raymond S.↗

Preliminary Simulations of the Ullage Dynamics in Microgravity during the Jet Mixing Portion of Tank Pressure Control Experiments

The results of CFD simulations of microgravity tank pressure control experiments performed on the Space Shuttle are presented. A 13.7 liter acrylic model tank was used in these experiments. The tank was filled to an 83 percent fill fraction with Freon refrigerant to simulate cryogenic propellants stored in space. In the experiments, a single liquid jet near the bottom of the tank was used for mixing the tank. Simulations at a range of jet Weber numbers were performed. Qualitative comparisons of the liquid and gas interface dynamics observed and recorded in the experiments and those computed are shown and discussed. The simulations were able to correctly capture jet penetration of the ullage, qualitatively reproduce ullage shapes and dynamics, as well as the final equilibrium position of the ullage.

jet mixing↗

Preliminary Simulations of the Ullage Dynamics in Microgravity During the Jet Mixing Portion of Tank Pressure Control Experiments

The results of CFD simulations of microgravity tank pressure control experiments performed on the Space Shuttle are presented. A 13.7 liter acrylic model tank was used in these experiments. The tank was filled to an 83 percent fill fraction with Freon refrigerant to simulate cryogenic propellants stored in space. In the experiments, a single liquid jet near the bottom of the tank was used for mixing the tank. Simulations at a range of jet Weber numbers were performed. Qualitative comparisons of the liquid and gas interface dynamics observed and recorded in the experiments and those computed are shown and discussed. The simulations were able to correctly capture jet penetration of the ullage, qualitatively reproduce ullage shapes and dynamics, as well as the final equilibrium position of the ullage.

simulation↗

Liquefaction Study of Gaseous Oxygen Inside Mars Ascent Vehicle Propellant Tank

The in-situ production of propellants for Mars missions will utilize carbon dioxide (CO2) in the Mars atmosphere to produce oxygen. The oxygen then needs to be cooled, liquefied, and stored to be available for Mars ascent propulsion, which could be up to 2 years after liquefaction starts. Recent investigations have demonstrated the feasibility of both achieving zero boiloff and controlling the pressure of oxygen within a tank using high-efficiency reverse turbo-Brayton-cycle cryocoolers. A tube-on-tank configuration is being studied in this work. The cooling fluid circulating in the cryocooler system is routed through a network of cooling tubes on the oxygen tank. The oxygen gas produced from the in-situ production process is introduced into the chilled tank. A series of analysis of this configuration has been performed to investigate the liquefaction rate inside the tank, the thermal gradient near the top of the tank where the oxygen gas feeding tubing is located. The analyses include 2D axisymmetric CFD analysis using ANSYS Fluent, 1D thermal analysis using Matlab, and 3D thermal analysis using MSC Patran/pthermal. These three models correlate and validate each other.

Liquefaction↗

Cryogenic Fluid Management Using CFD - Development and Applications

The Fluid Dynamics Branch has positioned itself to support a wide range of customers in need of Cryogenic Fluid Management (CFM) analysis. A computational fluid dynamics (CFD) tool used for all manner of internal and external propulsion applications has been extended and refined to better model cryogenic propellant storage and tanking operations. Through the CFM Portfolio project, several validation activities were initiated. Validation of propellant tank self-pressurization, autogenous pressurization, slosh-induced ullage collapse, and jet-induced mixing all aid in defining model accuracy. The on-going validation effort has enabled confident application of the tool to in-line design and evaluation of CFM hardware and operations. Recent project support included design of no-vent fill operations using various techniques for the Human Landing System (HLS). Propellant transfer dynamics were predicted to ensure complete propellant delivery. Many more activities in support of HLS and other projects and programs were supported. The branch has engaged the CFM community to share recent findings and capabilities through conferences. To reach the wider community, a consolidation of fundamental findings from work on in-space pulsed settling maneuvers was distilled in an AIAA paper. Development and demonstration of CFM modeling capabilities continues in order to meet the needs of the agency and its industry partners in the endeavor to sustainably reach the Moon and beyond.

Jacob Brodnick↗

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown↗

Computation of Effective Mechanical Properties and Mechanical Erosion Modeling of TPS Materials

The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.

Mechanical Properties↗

Comprehensive Forced Response Analysis of J2X Turbine Bladed-Discs with 36- Degree Variation in CFD Loading

At NASA/MSFC, Structural Dynamics personnel continue to perform advanced analysis for the turbomachinery in the J2X Rocket Engine, which is under consideration for the new Space Launch System. One of the most challenging analyses in the program is predicting turbine blade structural capability. Resonance was predicted by modal analysis, so comprehensive forced response analyses using high fidelity cyclic symmetric finite element models were initiated as required. Analysis methodologies up to this point have assumed the flow field could be fully described by a sector, so the loading on every blade would be identical as it travelled through it. However, in the J2X the CFD flow field varied over the 360 deg of a revolution because of the flow speeds and tortuous axial path. MSFC therefore developed a complex procedure using Nastran Dmap's and Matlab scripts to apply this circumferentially varying loading onto the cyclically symmetric structural models to produce accurate dynamic stresses for every blade on the disk. This procedure is coupled with static, spin, and thermal loading to produce high cycle fatigue safety factors resulting in much more accurate analytical assessments of the blades.

Elrod, David↗

Numerical Simulations For the F-16XL Aircraft Configuration

Numerical simulations of flow around the F-16XL are presented as a contribution to the Cranked Arrow Wing Aerodynamic Project International II (CAWAPI-II). The NASA Tetrahedral Unstructured Software System (TetrUSS) is used to perform numerical simulations. This CFD suite, developed and maintained by NASA Langley Research Center, includes an unstructured grid generation program called VGRID, a postprocessor named POSTGRID, and the flow solver USM3D. The CRISP CFD package is utilized to provide error estimates and grid adaption for verification of USM3D results. A subsonic high angle-of-attack case flight condition (FC) 25 is computed and analyzed. Three turbulence models are used in the calculations: the one-equation Spalart-Allmaras (SA), the two-equation shear stress transport (SST) and the kε turbulence models. Computational results, and surface static pressure profiles are presented and compared with flight data. Solution verification is performed using formal grid refinement studies, the solution of Error Transport Equations, and adaptive mesh refinement. The current study shows that the USM3D solver coupled with CRISP CFD can be used in an engineering environment in predicting vortex-flow physics on a complex configuration at flight Reynolds numbers.

Elmiligui, Alaa A.↗

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗

CFD and Experimental Data Comparisons for Conventional and AFC-Enabled CRM High-Lift Configurations

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of a representative conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface of the flap shoulder, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-Configuration to produce the aerodynamic performance (lift) comparable to that of the CRM-HL configuration over a large angle of attack range encompassing maximum lift conditions. For comparison purposes, computations are also per-formed for the CRM-HL configuration. In the current paper, PowerFLOWR©, a CFD code based on the Lattice Boltzmann method (LBM), is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with a combination of two different types of AFC actuators. The transonic version of the PowerFLOWR© code that has been validated for high-speed flows is used to simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. This study is focused on the AFC systems and actuator arrangements that had emerged based on the parametric studies conducted previously at the nominal landing condition. Comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations over a broad angle of attack range. The numerical solutions predict the expected trends in aerodynamic forces with angle of attack variation.

Veer N Vatsa↗

Predicting Rocket or Jet Noise in Real Time

A semi-empirical theoretical model and a C++ computer program that implements the model have been developed for use in predicting the noise generated by a rocket or jet engine. The computer program, entitled the Realtime Rocket and Jet Engine Noise Analysis and Prediction Software, is one of two main subsystems of the Acoustic Prediction/Measurement Tool, which comprises software, acoustic instrumentation, and electronic hardware combined to afford integrated capabilities for real-time prediction and measurement of noise emitted by rocket and jet engines. [The other main subsystem, consisting largely of acoustic instrumentation and electronic hardware, is described in Wireless Acoustic Measurement System, which appears elsewhere in this section.] The theoretical model was derived from the fundamental laws of fluid mechanics, as first was done by M. J. Lighthill in his now famous theory of aerodynamically generated sound. The far-field approximation of the Lighthill theory is incorporated into this model. Many other contributions from various researchers have also been introduced into the model. The model accounts for two noise components: shear noise and self noise. The final result of the model is expressed in terms of a volume integral of the acoustic intensities attributable to these two components, subject to various directivity coefficients. The computer program was written to solve the volume integral. The inputs required by the program are two data files from a computational fluid dynamics (CFD) simulation of the flow of interest: the computational-grid file and the solution file. The CFD solution should be one that has been obtained for conditions that closely approximate those of an experimental test that is yet to be performed. In the current state of development of the model and software, it is recommended that the observation points lie along a radius at an angle >60 from the jet axis. The software provides, and is driven via, a graphical user interface, which facilitates its use. Optionally, the program accepts additional input in the form of data on the measured sound pressure level as a function of frequency at a given far-field location, preferably at an angle of 90 from the jet axis. The user is prompted to use default empirical constants or to choose constants based the measurement data. The user can view the results and compare them with other computational or experimental data. Once satisfied with the results, the user can save a graph of the results in a file that can be imported into documents.

Frendi, Kader↗