Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD simulations”

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 487 records · Page 27

A Green’s Function Sensor Fusion Approach for Evaluating Spacecraft Entry Heating From on-Board Thermal Instrumentation

During atmospheric entry, distributed thermal measurements are critical to enable the evaluation of heat loads on spacecraft thermal protection systems (TPS). In recent space exploration missions, Schmidt-Boelter-type heat flux gauges have been integrated into the TPS alongside conventional temperature measurement instrumentation to measure total (convective and radiative) and radiative heat transfer rates1,2. While direct heat flux sensors (HFS) are able to provide valuable information detailing the thermal loads experienced by spacecraft, the interpretation of these measurements in unsteady, convective environments requires a correction factor to account for local heating augmentations at the cold wall HFS surface2,3. Current efforts to estimate cold wall correction factors, and thus recover the hot wall TPS heat flux, rely on time-marching computational fluid dynamics (CFD) simulations2. Simulation-based methods are susceptible to large uncertainties, however, as they require estimations of vehicle trajectory, gas kinetics, wall catalysis models, and other flight conditions as input parameters2,4. Furthermore, CFD simulations are computationally expensive and cannot efficiently survey all possible entry scenarios, exacerbating the uncertainty of reconstructed hot wall heat flux values. These drawbacks motivate the development of alternative hot wall heat flux reconstruction methods that are not reliant on CFD-based correction factors.

Kenneth McAfee↗

Numerical Simulation of Low-Density Shock-Wave Interactions

Computational Fluid Dynamics (CFD) numerical simulations of low-density shock-wave interactions for an incident shock impinging on a cylinder have been performed. Flow-field density gradient and surface pressure and heating define the type of interference pattern and corresponding perturbations. The maximum pressure and heat transfer level and location for various interaction types (i.e., shock-wave incidence with respect to the cylinder) are presented. A time-accurate solution of the Type IV interference is employed to demonstrate the establishment and the steadiness of the low-density flow interaction.

Glass, Christopher E.↗

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was generated in a previous low-boom multidisciplinary optimization (MDO) study that used computational fluid dynamics (CFD) off-body pressure to compute the undertrack sonic boom ground signature and calibrated low-fidelity aerodynamics analyses to compute the mission performance metrics. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, can carry 40 passengers for a low-boom overland mission with cruise Mach 1.7 for the airport pairs over the continental US and has the potential to achieve an undertrack sonic boom ground noise level below 70 perceived level of decibels at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) for the mission analysis and modeled as a flow-through nacelle for CFD-based sonic boom analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned in this paper after replacing the flow-through nacelles by CFD engines for sonic boom analysis using aeropropulsive CFD simulation. A process is developed for approximation of the NPSS engine at SOC by a CFD engine for aeropropulsive CFD simulation. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the outer mold line of the configuration with the CFD engines is redesigned to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept at SOC. Finally, the CFD simulation data for the redesigned concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of supersonic transports that satisfy both the low-boom and mission performance requirements.

Low-boom supersonic transport↗

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between liquid and vapor phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures, temperatures, and pressurant requirements agree with experimental data, with pressurant mass predictions within 16\% of observed values. This study highlights the importance of selecting appropriate turbulence models to accurately simulate complex flow and heat transfer phenomena during tank pressurization and expulsion. By enhancing the accuracy and reliability of CFD models for liquid hydrogen under cryogenic conditions, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗

Computational Fluid Dynamics Simulation of the Swordfish Oxygen/Methane Rotating Detonation Rocket Engine

Rotating Detonation Rocket Engines (RDREs) provide several theoretical performance benefits over traditional deflagration combustion rocket engines. RDREs can potentially generate higher specific impulse (ISP) for a given propellant combination enabling significant consumable mass savings for a spaceflight vehicle. Conversely, RDREs are a relatively new technology with a shallower analytical modeling base compared to traditional combustion engine architectures. In this work, we utilize time-accurate, Computational Fluid Dynamics (CFD) modeling to simulate the complex fluid environment inside an RDRE combustion chamber and propellant manifolds. The selected engine of study is the Swordfish RDRE which was developed and tested and at the National Aeronautics and Space Administration (NASA) Marshall Space Flight Center (MSFC). CFD simulations were conducted using the as-tested Swordfish geometry and conditions to facilitate comparison between the analytical predictions and experimental results.

Brian R. Richardson↗

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between vapor and liquid phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved autogenous pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures and temperatures agree well with experimental data, whereas the pressurant mass predictions are within 16% of observed values. This study emphasizes the importance of selecting appropriate turbulence models to accurately simulate the complex flow and heat transfer phenomena during tank autogenous pressurization and expulsion. By improving the accuracy and reliability of CFD models for these processes, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was designed for a cruise Mach of 1.7 and 40 passengers. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, was generated using computational fluid dynamics (CFD) based sonic boom analysis at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) and modeled as a flow-through nacelle for CFD analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned using an aeropropulsive CFD simulation. A process for approximation of the NPSS engine at SOC by a CFD engine is developed. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the configuration with the CFD engines is optimized to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept. Finally, the CFD simulation data for the optimized concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of low-boom supersonic transports.

Low-boom supersonic transport↗

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercial CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and or these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were run to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire left wing without the wheel. well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was nom1 to the spar. It included all o f the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas H.↗

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercia2 CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and for these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were m to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire lee wing without the wheel well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was normal to the spar. It included all of the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was Found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗

CFD2030 Grand Challenge: CFD-in-the-Loop Monte Carlo Simulation for Space Vehicle Design

Space vehicle design and certification differs widely from aircraft design relying more on probabilistic approaches than deterministic. Monte Carlo simulation plays an important role in the probabilistic design of space vehicles to ensure robust and reliable operation. Today, Monte Carlo flight simulation requires 1000’s of trajectory simulations that use databases to provide aerosciences models. These databases can be extremely expensive and time consuming to develop. Replacing these databases with unsteady computational fluid dynamics directly in the simulation loop has potential to significantly reduce the time required to analyze space vehicle concepts, improve simulation accuracy, and reduce the cost of space vehicle development. The CFD Vision 2030 Study outlined gaps and roadblocks to meeting the vision described in the study. The geometric, physical, and computational challenges associated with CFD-in-the-loop Monte Carlo simulation for space vehicle design are substantial and serve as an excellent grand challenge to advance the CFD 2030 vision.

CFD2030↗

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗

Advancing xEMU Lunar Dust Mitigation Devices

Fine, electrically charged, glass like dust particles caused significant damage to the Apollo EMU9 during lunar EVAs, identified as one of the greatest challenges to future exploration. Passive Lunar Dust Mitigation Devices (LDMD) were developed, within a SBIR Phase II, to prohibit this dust from interrupting venting space suit component operation. A Computational Fluid-Dynamics and Discrete Element Method (CFD-DEM) Simulation Tool was developed at the University of Colorado, Boulder to predict venting gas flow ability to self-clean adhered dust particles from LDMD surfaces. Lunar dust properties (i.e., adhesion and cohesion strengths) required to complete Simulation Tool analysis are relatively unknown due to considerable differences between the Earth and the Moon (i.e., gravity, humidity) and due to an absence of dust particles in their native state. Analytically determining gas flow velocity, density and direction within the fluid Boundary Layer, microns from LDMD surfaces presented a second challenge. Dusty Plasma Chamber testing is being performed at Auburn University, Auburn to observe electrostatically charged dust behavior as it adheres to LDMD prototypes and specific geometry and is then blown away by metered gas flow. Test articles were developed to offer insight into the impact of different flow geometries, surface roughness and dust removal within the gas Boundary-Layer. Observed dust behavior is currently being developed to support the CFD-DEM analysis. Many Simulation Tool analytical cases have been processed to support the intention of completing sensitivity studies to assess how different dust adherence values and Boundary Layer fluid properties impact LDMD self-cleaning effectivity.

Thomas J Stapleton↗

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very 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 used 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. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, 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 will be able to 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. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

Computational Modeling in Plasma Processing for 300 mm Wafers

Migration toward 300 mm wafer size has been initiated recently due to process economics and to meet future demands for integrated circuits. A major issue facing the semiconductor community at this juncture is development of suitable processing equipment, for example, plasma processing reactors that can accomodate 300 mm wafers. In this Invited Talk, scaling of reactors will be discussed with the aid of computational fluid dynamics results. We have undertaken reactor simulations using CFD with reactor geometry, pressure, and precursor flow rates as parameters in a systematic investigation. These simulations provide guidelines for scaling up in reactor design.

Meyyappan, Meyya↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Multiblock Navier-Stokes solutions about the F/A-18 wing-LEX-fuselage configuration

Three-dimensional thin-layer Navier-Stokes computations are presented for the F/A-18 configuration. The modeled configuration includes an accurate surface representation of the fuselage, leading-edge extension (LEX), and wing, both with and without leading-edge flap deflection. A multiblock structured grid strategy is employed to decompose the computational flowfield domain around the subject configuration. Steady-state solutions are obtained from an algorithm that solves the compressible Navier-Stokes equations with an upwind-biased, flux-difference splitting approach. The results presented are based on a fully turbulent flow assumption, simulating the high Reynolds number flow conditions that correspond to a recent FI A-18 flight experiment. Good agreements between the computations and the flight test results are obtained for both surface flow patterns as well as surface pressure distributions. Furthermore, a correlation between the computed LEX vortex-core and the flight test results, observed by way of smoke visualization, is also presented.

Aircraft configurations↗