Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,261 records · Page 70

Benchmarking Simulations of the Compton Spectrometer and Imager with Calibrations

The Compton Spectrometer and Imager (COSI) is a balloon-borne 𝛾-ray (0.2-5 MeV) telescope designed to study astrophysical sources. COSI employs a compact Compton telescope design utilizing 12 high-purity germanium double-sided strip detectors and is inherently sensitive to polarization. In 2016, COSI was launched from Wanaka, New Zealand and completed a successful 46-day flight on NASA’s new Super Pressure Balloon. In order to perform imaging, spectral, and polarization analysis of the sources observed during the 2016 flight, we compute the detector response from well-benchmarked simulations. As required for accurate simulations of the instrument, we have built a comprehensive mass model of the instrument and developed a detailed detector effects engine which applies the intrinsic detector performance to Monte Carlo simulations. The simulated detector effects include energy, position, and timing resolution, thresholds, dead strips, charge sharing, charge loss, crosstalk, dead time, and detector trigger conditions. After including these effects, the simulations closely resemble the measurements, the standard analysis pipeline used for measurements can also be applied to the simulations, and the responses computed from the simulations are accurate. We have computed the systematic error that we must apply to measured fluxes at certain energies, which is 6.3% on a rage. Here we describe the detector effects engine and the benchmarking tests performed with calibrations.

Sleator, Clio C.↗

How Should Life Support Be Modeled and Simulated?

Why do most space life support research groups build and investigate large models for systems simulation? The need for them seems accepted, but are we asking the right questions and solving the real problems? The modeling results leave many questions unanswered. How then should space life support be modeled and simulated? Life support system research and development uses modeling and simulation to study dynamic behavior as part of systems engineering and analysis. It is used to size material flows and buffers and plan contingent operations. A DoD sponsored study used the systems engineering approach to define a set of best practices for modeling and simulation. These best practices describe a systems engineering process of developing and validating requirements, defining and analyzing the model concept, and designing and testing the model. Other general principles for modeling and simulation are presented. Some specific additional advice includes performing a static analysis before developing a dynamic simulation, applying the mass and energy conservation laws, modeling on the appropriate system level, using simplified subsystem representations, designing the model to solve a specific problem, and testing the model on several different problems. Modeling and simulation is necessary in life support design but many problems are outside its scope.

Jones, Harry W.↗

Simulating the Benesov Bolide Flowfield and Spectrum at Altitudes of 47 and 57 km

This paper develops a computational fluid dynamics capability for simulating the radiative emission from a meteor shock-layer and wake to an external observer. The developed capability includes the impact of radiation and ablation on the meteor flowfield, where coupled radiation accounts for the impact of radiation on the flowfield energy equations and coupled ablation accounts for the injection of ablation species from the meteoroid surface into the ow. This capability includes updated flowfield chemistry and nonequilibrium radiation models, which are compiled from the literature. To provide a level of validation for this capability, the Benesov spectral measurements are considered. Although the meteoroid was likely fragmented at the 47 and 57 km altitudes considered, and the measurement uncertainty is roughly ±50%, these measurements represent the best available spectral measurements for a relatively large (meter-class) bolide. To determine the equivalent meteoroid diameter to simulate, the continuum component of the spectrum in the 570 to 610 nm range is considered. This range is dominated by air emission from the high pressure and temperature shock-layer in front of the meteoroid, which is simulated with a relatively small uncertainty (because of the abundance of relevant shock tube measurements) and sensitive to the equivalent meteoroid diameter. This allows an equivalent diameter between 0.62 to 0.9 m to be determined, based on comparisons between measurements and simulations in the 570 to 610 nm range. The range of equivalent diameters is due to the various meteoroid geometries considered, where increasing the bluntness decreases the required diameter. Applying these diameters to simulations at 47 and 57 km results in simulated spectra that compare within 30% of the measured values, considering the wavelength integrated values between 400 and 650 nm. This agreement is well within the estimated measurement uncertainty of ±50%. All major spectral features are captured by the simulations. This agreement provides a level of validation for the developed model that has not been previously available.

Luminous efficiency↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. The quadrupole method has been shown as a viable technique for rapidedly simulating the thermal repsonse of layered systems. Often, the effort has focused on a one-dimentional models, in particular for improved analysis of thermal data. For composites, the in-plane heat diffusion often significantly impacts the thermal response of defects of interest, therefore three-dimensional simulations are desirable. This paper discusses the extension of the quadrupole methodology to perform simulations of thermographic responses in three-dimensional configuration. This enables the more realistic simulations of the thermal response of delaminiations in composites. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are known obtained from other sources.

Thermography↗

A Distributed Simulation-to-Flight Framework to Support Investigating Trust/Trustworthiness in Multi-Agent Systems

As autonomous systems continue to grow both in use and complexity, the necessity for robust and extensible simulation-to-flight frameworks is paramount for establishing an effective architecture for autonomous systems. Hardware test flights are time-consuming and cost prohibitive during early system design and development. Simulation environments can be useful tools to accelerate algorithm development and testing. However, transitions from simulation to flight (sim-to-flight) can be challenging, unless systems are designed with this transition in mind and with the necessary capabilities built into the architecture and framework. One of the objectives of Autonomy Teaming and TRAjectories for Complex Trusted Operational Reliability (ATTRACTOR) was to design and develop a distributed mixed-reality simulation environment to begin establishing a basis for certification of autonomous systems via research into trust and trustworthiness. ATTRACTOR’s objective was to construct computational concepts of trustworthiness and justifiable trust in multi-agent autonomous teams, to inform future certification of safety-critical and time-critical autonomous systems in aviation. In this paper, we present an autonomous systems architecture and development framework paired with a persistent distributed modeling and simulation (ModSim) environment for test and evaluation of autonomous systems. They were designed under ATTRACTOR in order to measure and establish trustworthiness and trust in single-and multi-agent human-machine systems whether these machines are fixed-wing general aviation, rotary-wing Unmanned Aerial Vehicles (UAVs), ground rovers, or even spacecraft. The Autonomous Entity Operational Network (AEON) framework enables autonomous system development with an easily extensible collection of libraries and plug-n-play nodes facilitated by the Data Distribution Service (DDS) communication protocol standard. The Baseline Environment for Autonomous Modeling (BEAM) simulation environment is a distributed mixed-reality Unity™-based environment built around the same DDS communication paradigm allowing for easy integration with AEON-based autonomous applications, enabling sim-to-flight with minimal configuration changes. Using AEON and BEAM, source code that runs in simulation ports directly to hardware and has successfully flown in the National Airspace System (NAS) at NASA LaRC many times over the lifetime of ATTRACTOR.

Benjamin N Kelley↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Autonomous Lunar PNT Simulator (ALPS)

In this paper we present a simulation tool that has been used to develop an architecture for a low-cost Position, Navigation and Time (PNT) system to provide PNT services on and around the Moon using non-dedicated low-cost orbital and ground assets. The simulation tool has been developed to be flexible and is capable of modeling and analyzing the many different capabilities and configurations that the non-dedicated assets could support. The tool models the creation of an ad hoc swarm, the localization of this swarm and the subsequent provision of PNT services from this swarm. We present results from several studies of select configurations chosen to reflect existing and future real-world needs and capabilities. Over the next few decades there is expected to be a substantial increase in Lunar missions supporting and inspired by NASA’s Artemis Program. It is expected that a large fraction of these missions will be low cost, utilizing rideshares and CubeSats, just as has seen in Earth orbit over the past decade. Many of these missions will need navigation capabilities but may be unable to support the large power, mass and weight that a weak GNSS or DSN based navigation solution would entail. As Lunar PNT service demand is not likely to be needed over the entire lunar surface 24/7, the creation of a dedicated Lunar GNSS constellation cannot be justified. Hence, our proposed architecture envisions utilizing existing Lunar science and exploration assets to create ad hoc and on demand Lunar PNT swarms capable of providing PNT services to low cost missions. The simulation reflects the two distinct parts of the architecture: the creation of a Lunar PNT swarm using non-dedicated existing assets, and the quantitative modeling of the quality of the PNT services that this swarm provides. To support the former, the simulation supports various swarm localization techniques, including centralized and distributed EKFs both of which support pluggable dynamics models. In modeling the PNT service performance the simulation adopts standard techniques from the GNSS community, including providing degree of precision (DOP) estimates for theoretical end users. Crucially, all asset capabilities, including clock accuracy, independent location self-knowledge and timing measurement precision can be set independently for each asset, reflecting the key concept of utilizing non-dedicated assets. The simulation is predominantly implemented in MATLAB, with GMAT being used for the propagation of orbital assets. The paper will present results from the simulation reflecting the tool’s flexibility and focusing on scenarios that match real-world proposed missions, including scenarios designed to provide PNT support to lunar surface missions similar to NASA Ames' forthcoming VIPER mission. The performance of centralized and decentralized swarm configuration and localization techniques will be compared. Finally, performance of the PNT service provided by an autonomous Lunar PNT swarm will be compared to existing radiometric and weak GNSS methods.

Kelley Elizabeth Hashemi↗

Multibody Based Digital Astronaut Dynamics Simulation

BACKGROUND: This study provides the Software, Robotics, & Simulation Division at the NASA Johnson Space Center with a verification tool for multibody dynamics simulation requiring human motion. The motivation stems from current studies of several Vibration Isolation & Stabilization(VIS)system designs that attenuate the moments and forces which would be transmitted to a spacecraft during an exercise. A multibody dynamics model for a proposed VIS was available previously[1], therefore modeling of the VIS was not needed for this work. The interest here is in creating the multibody dynamics model of an astronaut in motion which may be utilized independently or while attached to a mechanism. An existing simulation [2] that utilizes OpenSim [3,4] and an in-house multibody dynamics package (MBDyn) [5] is used in order to verify the astronaut model. The main advantage this model will have over the existing simulation is that everything will be processed in one tool. METHODS AND RESULTS: Creating the simulation required; estimation of Body Segment Inertial Parameters (BSIP),a multibody model of the human-VIS system, joint acceleration profiles, and input files for MBDyn, which is used for this analysis. The scaling factors provided by Dumas et al. [6] are utilized in estimating the BSIP. Anthropometric data are used for estimating these parameters, the Anthropometric Survey of US Army Personnel (ANSUR II) [7] was the source. The astronaut model consists of 15 bodies, 14 joints and 32 degrees of freedom, with the dynamics topology generated using the center of mass locations and anthropometric data. MBDyn has an option for prescribed joint motion (PJM), which requires joint acceleration data as input. The joint angle data is first obtained from a motion capture system and then processed through code that has been created to generate approximate joint acceleration profiles. The topology tree for the astronaut model begins at the right foot up to the pelvis where there is one branch for going down the left leg and another for the torso. The torso branch leads to branches for the arms and a leaf body for the head/neck segment. For attachment to the VIS, the heel of the right foot is connected to the VIS platform through a fixed joint, resembling a foot restraint. The left foot does not attach to the platform in order to prevent a system with a closed loop. Topology and symmetry of the astronaut model were verified through kinematic analysis. Further verification of the forces and moments transmitted to the VIS were verified against the existing simulation. There was a satisfactory level of agreement when testing a simple motion, for example, rocking back and forth. Full exercise motions are to be tested soon. The main outcome has been a novel application of MBDyn for biomechanics modeling that is now available for dynamic simulations involving human motion. The estimation of BSIP was another useful result of this study, requiring only 15 inputs for generating mass properties of a theoretical astronaut model. Expansion on this work is possible by going through an alternative route in obtaining the joint motion data. Instead of high-tech and often expensive motion capture systems, an individual may watch videos with high focus and at a slow motion for each individual segment in order to determine the initial and final time and angle for that specific degree of freedom. Synthetic trajectories may also be created if there is no video reference available.

F N Matari↗

Investigations of Pole Erosion Mechanisms in the 12.5 kW HERMeS Hall Thruster Using Numerical Simulations and Ion Velocity Measurements

When the Hall Effect Rocket with Magnetic Shielding (HERMeS) was tested at the NASA GRC’s VF-5 the erosion rate on the pole covers was found to be at least a factor of two higher at 300 V and 20.8 A than at other operating conditions in the 400 V to 600 V range. Simulations using our hydrodynamics code Hall2De do not predict increased erosion rates at 300 V but accurately compute erosion rates similar to the measurements at the other operating conditions. We investigate the source of the discrepancy between measurements and simulations at 300 V using a combination of numerical simulations and experimental measurements of the ion velocity fields in the acceleration region of the thruster. By examining previous simulations that predicted the measured erosion rates for other thrusters and operating conditions, we determine that sputtering of the pole surfaces by high energy ions is the most likely mechanism behind the erosion rates at 300 V. High energy ions sputter the pole surfaces when the acceleration region of the thruster is downstream of the pole surface plane. When the latter occurs, the curvature of the plasma potential contours at the edges of the channel accelerates a small fraction of the high energy ions radially. We find that neither our Hall2De simulations nor the experimental measurements produce ions of sufficiently high energy and flux to the inner pole cover to explain the measured erosion there. Only at the outer edge of the cover is where we find the ions needed to yield simulation results that are comparable to the erosion measurements. We argue that one possible mechanism for the higher erosion is high energy ions that graze the channel corner of the pole and become trapped in a sheath that develops downstream. Within the sheath, we find that the electric field is large enough to turn the ions towards the surface of the pole. A simplified calculation shows that the erosion rates produced by sputtering of trapped ions are similar to the measurements. We also propose that the presence of local oscillations and high-energy cathode ions may be sources of the enhanced erosion We conclude this article by proposing a combination of analysis, simulation, and experimental measurements that can be used to address the validity of these hypotheses.

Mikellides, Ioannis G.↗

POLARRIS: A POLArimetric Radar Retrieval and Instrument Simulator

This paper introduces a synthetic polarimetric radar simulator and retrieval package, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), for evaluating cloud‐resolving models (CRMs). POLARRIS is composed of forward (POLARRIS‐f) and inverse (retrieval and diagnostic) components (iPOLARRIS) to generate not only polarimetric radar observables (Zh, Zdr, Kdp, ρhv) but also radar‐consistent geophysical parameters such as hydrometeor identification, vertical velocity, and rainfall rates retrieved from CRM data. To demonstrate its application and uncertainties, POLARRIS is applied to simulations of a mesoscale convective system over the Southern Great Plains on 23 May 2011, using the Weather Research and Forecasting model with both spectral bin microphysics (SBM) and the Goddard single‐moment bulk 4ICE microphysics. Statistical composites reveal a significant dependence of simulated polarimetric observables (Zdr, Kdp) on the assumptions of the particle axis ratio (oblateness) and orientation angle distributions. The simulated polarimetric variables differ considerably between the SBM and 4ICE microphysics in part due to the differences in their ice particle size distributions as revealed by comparisons with aircraft measurements. Regardless of these uncertainties, simulated hydrometeor identification distributions overestimate graupel and hail fractions, especially from the simulation with SBM. To minimize uncertainties in forward model, the particle shape and orientation angle distributions of frozen particles should be predicted in a microphysics scheme in addition to the size distributions and particle densities.

Toshi Matsui↗

Comparison between Hayabusa 2 Spectral Measurements and Simulations

This paper compares recent Hayabusa 2 spectral measurements with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the coupled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator, based on the limited available published information. The comparison between the simulations and measurements focuses on two atomic nitrogen lines and two atomic oxygen lines, as well as the CN Violet band system. For the atomic lines, the measurements and simulations agree within 25% over most of the trajectory. This excellent agreement is unprecedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield/radiation modeling. To provide a link between these observed radiation measurements and the radiative heating to a vehicle surface, the recently developed flowfield property binning approach is used to identify flowfield properties that provide the dominant emission contribution to the measured spectrum. These identified flowfield properties are then shown to match those for the surface radiative heating to the currently developed Mars Sample Return (MSR) Earth Entry System (EES) at a specific trajectory point and surface location. This matching of flowfield properties indicates that this EES radiative heating and Hayabusa 2 observed radiation simulations are equivalent radiation problems. Therefore, the excellent agreement between Hayabusa 2 measurements and simulations may be leveraged to inform the radiation heating margin for the MSR EES.

Christopher O Johnston↗

Validation of Small Satellite Dynamics Simulation Modules using ASTERIA Flight Data

ASTERIA (Arcsecond Space Telescope Enabling Research in Astrophysics) was a CubeSat space telescope that operated in low-Earth orbit, having been deployed from the International Space Station in 2017. The spacecraft has achieved sub-arcsecond pointing stability and millikelvin thermal stability over 20-minute observations. A key enabling feature of the ASTERIA mission is the low level of pointing error that has been demonstrated under ASTERIA’s fine pointing control mode. Prior to launch, the ASTERIA mission performed analysis and simulations to estimate the in-flight pointing performance. This analysis used reaction wheel models provided by JPL’s Small Satellite Dynamics Testbed (SSDT), which had performed Kistler table testing to characterize the jitter caused by the reaction wheels of the Blue Canyon Technologies XACT attitude control unit within ASTERIA. These models were subsequently incorporated into the SSDT simulation. The main source of jitter, the high-frequency attitude disturbance over the camera exposure time, on ASTERIA is the set of rotating reaction wheels that are used for maintaining fine pointing during observations. The reaction wheels impart disturbance forces and torques continuously that cause unwanted motion during imaging. Increased jitter consequently leads to blurring of the image. Therefore, to benefit all SmallSat missions devoted to photometric or spectroscopic astrophysics applications, an on-orbit flight data acquisition and jitter testing campaign was performed in an attempt to help validate the SSDT’s simulation models. This paper describes the process used to validate several simulation models and attempt to characterize the ASTERIA jitter by analyzing the size of the resulting spot size for the target stars as a function of wheel speed, having commanded four wheel speeds for each of three stars of known brightness. After summarizing the on-orbit observation and jitter level measurement process, this paper compares the results obtained from the flight operational environment to the corresponding set of simulation jitter levels. Though the flight data obtained during this experiment was not sufficient to validate the simulation’s jitter model, other models have been validated with flight data, including the magnetic field, orbit propagation, sun position, and long-duration orbit decay models, which may be used by all other projects that use the SSDT’s simulation to increase system capability knowledge.

Mohan, Swati↗

A Simulation Architecture for Air Traffic Over Urban Environments Supporting Autonomy Research in Advanced Air Mobility

NASA is conducting investigations into Advanced Air Mobility (AAM) concepts, aircraft, and operations. One of the most challenging scenarios for AAM will be enabling safe routine access into densely populated urban centers. To address challenges in the urban environment, a moderately high-fidelity simulation capability is needed to investigate AAM flight operations over a regional area for the development and evaluation of autonomy technologies. This paper describes a system to simulate flight operations around regions such as the San Francisco-Oakland Bay area at a moderately-high scale (10's-100's of flights) that incorporates detailed vehicle models and control necessary to support research in airborne autonomy. The flight vehicle utilizes NASA AAM concept vehicle dynamics integrated with a custom flight management system and flight control system to accurately simulate all phases of flight. Glass cockpit displays have been developed for monitoring aircraft operation over a detailed simulated urban environment. Simulation models have been integrated to simulate air and ground-based sensors, such as radar and LIDAR. The commercial X-Plane software package is used as a rendering engine to mimic vision-based sensors (such as onboard and ground-based cameras) at various times of day and in various weather conditions over a relatively detailed graphical model of the city. The paper presents a detailed illustration of the simulation and software architecture used for traffic over this urban region. This system is enabling the evaluation of NASA research concepts in autonomy for urban AAM operations on the path toward aircraft flight test evaluation.

Keerthana Kannan↗

Qualitative Comparison of Pathfinder Three-Dimensional Transient Gas Granular Simulation Results to Imagery from a Subscale, Reduced Pressure Plume Surface Interaction Ground Test

NASA’s Game Changing Development Program, funded by NASA’s Space Technology Mission Directorate, the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~600 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the PSI Project also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the first Pathfinder three-dimensional, transient simulation of the test configuration used to acquire erosion and cratering test data as part of the Physics Focused Ground Test conducted by the PSI Project in FY2021. This first pathfinder simulation was conducted using the Loci/GGFS (Gas Granular Flow Solver), a computational fluid dynamics tool developed in part with funding provided by the PSI Project. Crater growth and ejecta flow paths from the pathfinder simulation is compared to raw high-speed imagery PFGT results of the crater growth. There is significant similarity of test and simulation to warrant expending resources on a more rigorous quantitative Validation Assessment.

J S West↗

Qualitative Comparison of Pathfinder Three-Dimensional Transient Gas Granular Simulation Results to Imagery from a Subscale, Reduced Pressure Plume Surface Interaction Ground Test

NASA’s Game Changing Development Program, funded by NASA’s Space Technology Mission Directorate, the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~600 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the PSI Project also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the first Pathfinder three-dimensional, transient simulation of the test configuration used to acquire erosion and cratering test data as part of the Physics Focused Ground Test conducted by the PSI Project in FY2021. This first pathfinder simulation was conducted using the Loci/GGFS (Gas Granular Flow Solver), a computational fluid dynamics tool developed in part with funding provided by the PSI Project. Crater growth and ejecta flow paths from the pathfinder simulation is compared to raw high-speed imagery PFGT results of the crater growth. There is significant similarity of test and simulation to warrant expending resources on a more rigorous quantitative Validation Assessment.

Jeff West↗

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre-and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre-and post-buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the span wise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗