Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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 991 records · Page 55

Surface Measurements of Solar Spectral Radiative Flux in the Cloud-Free Atmosphere

Recent studies (Charlock, et al.; Kato, et. al) have indicated a potential discrepancy between measured solar irradiance in the cloud-free atmosphere and model derived downwelling solar irradiance. These conclusions were based primarily on broadband integrated solar flux. Extinction (both absorption and scattering) phenomena, however, typically have spectral characteristics that would be present in moderate resolution (e.g., 10 nm) spectra, indicating the need for such measurements to thoroughly investigate the cause of any discrepancies. The 1996 Department of Energy Atmospheric Radiation Measurement Program (ARM) Intensive Observation Period (IOP), held simultaneously with the NASA Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) Program, provided an opportunity for two simultaneous but independent measurements of moderate resolution solar spectral downwelling irradiance at the surface. The instruments were the NASA Ames Solar Spectral Flux Radiometer and the Analytical Spectral Devices, Inc., FieldSpecT-FR. Spectral and band integrated quantities from both sets of measurements will be presented, along with estimates of the downwelling solar irradiance from band model and line by line calculations, in an effort to determine the compatibility between measured and calculated solar irradiance in the cloud-free atmosphere.

Pilewskie, Peter↗

The effectiveness of D 2 pellet injection in reducing intra-ELM and inter-ELM tungsten divertor erosion rates in DIII-D during the Metal Rings Campaign

Abstract Edge localized modes (ELMs) in H-mode plasmas erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes D 2 pellet injection on the DIII-D fusion experiment used as an ELM mitigation technique applied during the 2016 tungsten Metal Rings Campaign to reduce W erosion during ELMs. The 400.9 nm photon wavelength line emission intensity of tungsten atoms (WI) filterscope channels and Langmuir probes were used to infer the gross erosion rate of tungsten-coated tiles installed in the divertor of DIII-D. D 2 mass injection rates ranging from 34 to 41 arbitrary units (A.U.) and no D 2 injection resulted in a similar total W erosion rate during ELMs (intra-ELM). On average, results show a 29% increase in the total gross W erosion rate with intermediate mass injection rates (∼13–23 A.U.) compared to the no pellets and the highest injection rate cases. On average, the fast D 2 mass injection rate cases had 15% less erosion in the inter-ELM phase than the case with no pellets. Generally, higher D 2 mass injection rates increased the ELM frequency, and the highest injection rates reduced the average erosion per ELM and fractional carbon impurities at the top of the pedestal by nearly 40% when compared to the no-pellet case. As expected, a higher D 2 pellet injection rate led to a higher plasma density and lower plasma temperature in the divertor. Additionally, an increasing divertor inter-ELM plasma electron density directly correlated to more frequent pellet injection and a decrease in both the average gross intra-ELM W erosion and the total gross intra-ELM W erosion rate. Simulations of intra-ELM erosion using the ‘free-streaming plus recycling model’ (FSRM) underestimate W erosion during pellet injection by about 30% on average. The discrepancies between the experimental measurements and the FSRM intra-ELM W erosion predictions are postulated to be due to C/W material mixing. A simple analytic mixed-material model is presented and results in better agreement with the experimental data. These results highlight the importance of incorporating the effects of a mixed-material layer in the analysis of PFC erosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

vaned diffuser↗

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

centrifugal compressor↗

Benchmarking of three DWM-based wake models at below-rated wind speeds

Wind turbine wake models are essential tools for predicting power losses and structural loads in wind farms. Among these, the dynamic wake meandering (DWM) model, included as a recommended approach in the International Electrotechnical Commission design standard, is a widely used engineering-fidelity method that balances accuracy and computational cost. This study compares the performance of three DWM-based wake model implementations (from the Technical University of Denmark, the National Renewable Energy Laboratory, and the Institute for Energy Technology) under below-rated wind speed conditions. Model predictions of wake flow, power output, and structural loads for a four-turbine row are evaluated across different ambient turbulence levels and wind-direction misalignments and compared against high-fidelity large-eddy simulation results. All three models captured the overall wake evolution and mean turbine performance with reasonable accuracy; their predicted time-averaged thrust and power were typically within 5 %–10 % of the large-eddy simulation benchmark. However, notable differences emerged in wake structure and unsteady load predictions, with discrepancies increasing for turbines further downstream. These differences highlight the importance of modelling choices such as wake summation and turbulence treatment, which strongly influence power-deficit and fatigue-load predictions. Comparison with large-eddy simulations reveals each approach's strengths and weaknesses, indicating where improvements are needed. Overall, the findings point to specific refinements for DWM models to improve their fidelity, ultimately enabling more robust wake predictions for wind farm design and operation.

17 WIND ENERGY↗

Comparing Global Atmospheric CO2 Flux and Transport Models with Remote Sensing (and Other) Observations

We report recent progress derived from comparison of global CO2 flux and transport models with new remote sensing and other sources of CO2 data including those from satellite. The overall objective of this activity is to improve the process models that represent our understanding of the workings of the atmospheric carbon cycle. Model estimates of CO2 surface flux and atmospheric transport processes are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, to provide the basic framework for carbon data assimilation, and ultimately for future projections of carbon-climate interactions. Models can also be used to test consistency within and between CO2 data sets under varying geophysical states. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 2000 through 2009. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at 1x1 degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-3), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to remote sensing observations from TCCON, GOSAT, and AIRS as well as relevant in situ observations. Examples of the influence of key process representations are shown from both forward and inverse model comparisons. We find that the model can resolve much of the synoptic, seasonal, and interannual variability in the observations, although reasons for persistent discrepancies in northern hemisphere vegetation uptake are examined. At this time, we do not find any serious shortcomings in the model transport representation, but this is still the subject of close scrutiny. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in CO2 fluxes and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗

Using MCNP6 to Estimate Fission Neutron Properties of a Reflected Plutonium Sphere

The purpose of this project was to determine the fission multiplicity distribution, p(v), for the Beryllium Reflected Plutonium (BeRP) ball and to determine whether or not it changed appreciably for various High Density Polyethylene (HDPE) reflected configurations. The motivation for this project was to determine whether or not the average number of neutrons emitted per fission, v, changed significantly enough to reduce the discrepancy between MCNP6 and Robba, Dowdy, Atwater (RDA) point kinetic model estimates of multiplication. The energy spectrum of neutrons that induced fissions in the BeRP ball, NIF (E), was also computed in order to determine the average energy of neutrons inducing fissions, NIF . p(v) was computed using the FMULT card, NIF (E) and NIF were computed using an F4 tally with an FM tally modifier (F4/FM) card, and the multiplication factor, k eff , was computed using the KCODE card. Although NIF (E) changed significantly between bare and HDPE reflected configurations of the BeRP ball, the change in p(v), and thus the change in v, was insignificant. This is likely due to a difference between the way that NIF is computed using the FMULT and F4/FM cards. The F4/FM card indicated that NIF (E) was essentially Watt-fission distributed for a bare configuration and highly thermalized for all HDPE reflected configurations, while the FMULT card returned an average energy between 1 and 2 MeV for all configurations, which would indicate that the spectrum is Watt-fission distributed, regardless of the amount of HDPE reflector. The spectrum computed with the F4/FM cards is more physically meaningful and so the discrepancy between it and the FMULT card result is being investigated. It is hoped that resolving the discrepancy between the FMULT and F4/FM card estimates of NIF(E) will provide better v estimates that will lead to RDA multiplication estimates that are in better agreement with MCNP6 simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of Turbulence in Wind-Tunnel Measurements

This paper gives some quantitative measurements of wind tunnel turbulence and its effect on the air resistance of spheres and airship models, measurements made possible by the hot wire anemometer and associated apparatus in its original form was described in Technical Report no. 320 and some modifications are presented in an appendix to the present paper. One important result of the investigation is a curve by means of which measurements of the air resistance of spheres can be interpreted to give the turbulence quantitatively. Another is the definite proof that the discrepancies in the results on the N. P. L. Standard airship models are due mainly to differences in the turbulences of the wind tunnels in which the tests were made. An attempt is made to interpret the observed results in terms of the boundary layer theory and for this purpose a brief account is given of the physical bases of this theory and of conceptions that have been obtained by analogy with the laws of flow in pipes.

Dryden, H L↗

Wind Turbine Aeroelastic Stability in OpenFAST

Wind turbines are growing in size and increasingly suffer from aeroelastic instabilities. Unfortunately, numerical models often show inconsistent results during verification studies. We address this gap by first introducing novel linearization capabilities within the open-source aero-hydro-servo-elastic framework OpenFAST. Next, a code-to-code benchmark study is presented that compares modal parameters between OpenFAST and HAWCStab2 for a land-based version of the International Energy Agency 15-MW reference wind turbine modeled with quasi-steady aerodynamics. The two solvers are in strong agreement except for discrepancies in the second rotor flapwise modes. The differences are attributed to the torsional flexibility of the tower, which is assumed torsionally stiff in the OpenFAST model. Work is ongoing to close this modeling gap. The aeroelastic stability of a low-specific-power land-based wind turbine is also investigated. The impact of design choices is discussed, high-lighting how narrow the margins are between a stable design and an unstable design.

17 WIND ENERGY↗

Al III, Si IV, and C IV absorption toward zeta Ophiuchi: Evidence for photionized and collisionally ionized gas

We present Goddard High-Resolution Spectrograph observations at 3.5 km/s resolution and signal-to-noise ratios of 30 to 60 for the Al III, Si IV, and N V absorption lines in the far-ultraviolet spectrum of the O9.5 V star zeat Ophiuchi. The measurement reveal three types of highly ionized gas along the 140 pc line of sight. (1) Narrow components of Al III (b = 4.3 km/s, the mean value of (v(helio)) = -7.8 km/s; b = 3.2 km/s, the mean value of (v(sub helio)) = -14.4 km/s) and Si IV (b = 5.3 km/s, the mean value of (v(sub helio)) = -15.0 km/s) trace photionized gas in the expanding H II region surrounding zeta Oph. The observed magnitude and direction of the velocity offset between the Al III and Si IV profiles can be explained by models of H II regions that incorporate expansion. Narrow C IV absorption associated with the H II region is not detected. Predictions of the expected amounts of Si IV and C IV overestimate the column densities of these ions by factors of 30 and more than 10, respectively. The discrepancy may be due to the effects of elemental depletions in the gas and/or to the interaction of the stellar wind with surrounding matter. (2) Broad (b = 15 to 18 km/s) and weak Si IV and C IV absorption components are detected near the mean value of (v(sub helio)) = -26 km/s. The high-ionization species associated with these absorption components are probably produced by electron collisional ionization in a heated gas. This absorption may be physically related to the zeta Oph bow shock ot to a cloud complex situated within the local interstellar medium at d less than 60 pc. The C IV to Si IV column density ratio in this gas is 8, a factor of 6 less than conductive interface models predict, but this discrepancy may be removed by considering the effects of self-photoionization within the cooling gas in the model calculations. (3) A broad (b = 13 km/s) and weak C IV absorption feature detected at the mean value of (v(sub helio)) = -61 km/s is not seen in other species. We tentatively ascribe this absorption to gas in a postshock ragion of an optically thin shock in the zeta Oph stellar wind.

Sembach, Kenneth R.↗

Comparison of Comet Enflow and VA One Acoustic-to-Structure Power Flow Predictions

Comet Enflow is a commercially available, high frequency vibroacoustic analysis software based on the Energy Finite Element Analysis (EFEA). In this method the same finite element mesh used for structural and acoustic analysis can be employed for the high frequency solutions. Comet Enflow is being validated for a floor-equipped composite cylinder by comparing the EFEA vibroacoustic response predictions with Statistical Energy Analysis (SEA) results from the commercial software program VA One from ESI Group. Early in this program a number of discrepancies became apparent in the Enflow predicted response for the power flow from an acoustic space to a structural subsystem. The power flow anomalies were studied for a simple cubic, a rectangular and a cylindrical structural model connected to an acoustic cavity. The current investigation focuses on three specific discrepancies between the Comet Enflow and the VA One predictions: the Enflow power transmission coefficient relative to the VA One coupling loss factor; the importance of the accuracy of the acoustic modal density formulation used within Enflow; and the recommended use of fast solvers in Comet Enflow. The frequency region of interest for this study covers the one-third octave bands with center frequencies from 16 Hz to 4000 Hz.

Grosveld, Ferdinand W.↗

On the Role of Airborne Ice Nucleating Particles in Primary and Secondary Ice Formation Processes in Convective Midlatitude Clouds

Abstract Formation of ice in mixed-phase clouds may proceed via primary nucleation and secondary ice production (SIP). Primary nucleation involves the presence of ice nucleating particles (INPs), while SIP follows primary nucleation, resulting in ice crystal number concentrations N ice exceeding the number of INPs. The ice formation pathways in young congestus clouds are not well understood nor are the environmental conditions favorable for the onset of SIP. Coincident airborne measurements of INPs and N ice collected during the Secondary Production of Ice in Cumulus Experiment (SPICULE) campaign over the U.S. central Great Plains are reported for young congestus clouds. Number concentrations of INPs and N ice at temperatures between −10° and −20°C were used to categorize whether SIP was active in each analyzed cloud pass. Observational data suggested that fragmentation of freezing droplets (FFDs) was occurring and the most likely SIP mechanism for ice enhancement during the early stages of the sampled clouds, where higher cloud-base temperatures and stronger updrafts produced conditions favorable for the onset of FFD. Simple model simulations of a congestus cloud were used to compare predictions of N ice from newly implemented SIP mechanisms (including FFD) alongside those from the existing Hallett–Mossop (HM) mechanism. Although predicted N ice , when all SIP mechanisms were active, was in good agreement with observations, the dominant SIP mechanism in the model was predicted to be HM. The rate of heterogeneous freezing in raindrops was likely inhibiting realistic FFD rates in the simulated cloud. These discrepancies between the observations and simulations underscore the need for extended laboratory studies to formulate improved model representations of SIP in convective clouds.

Patnaude, Ryan J. [Department of Atmospheric Scien↗

Comments concerning magnetic fields as a possible cause of rapid, irregular variability of early-type stars

It is noted that it is helpful to divide the atmosphere of an early-type star into two parts: a photosphere and a mantle. The photosphere can be modeled satisfactorily by normal model-atmosphere procedures. When modeling the mantle, one must take account of the deposition of nonradiative heat and momentum and one should recognize that the mantles of hot stars do not appear to be uniform or spherically symmetric. Five areas of discrepancy between classical theory and observation are noted and a possible way of modeling the mantles of early-type stars is outlined. The model consists of arcades of magnetic loops which form helmet-type structures in the equatorial band of the star, and of coronal-hole-type structures emmanating from weak unipolar regions which are chiefly distributed at polar latitudes.

Underhill, A. B.↗

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗

Further Investigation of the Support System Effects and Wing Twist on the NASA Common Research Model

An experimental investigation of the NASA Common Research Model was conducted in the NASA Langley National Transonic Facility and NASA Ames 11-foot Transonic Wind Tunnel Facility for use in the Drag Prediction Workshop. As data from the experimental investigations was collected, a large difference in moment values was seen between the experiment and computational data from the 4th Drag Prediction Workshop. This difference led to a computational assessment to investigate model support system interference effects on the Common Research Model. The results from this investigation showed that the addition of the support system to the computational cases did increase the pitching moment so that it more closely matched the experimental results, but there was still a large discrepancy in pitching moment. This large discrepancy led to an investigation into the shape of the as-built model, which in turn led to a change in the computational grids and re-running of all the previous support system cases. The results of these cases are the focus of this paper.

Rivers, Melissa B.↗

Sockeye Validation Support Using the SPHERE Facility

The engineering-scale heat pipe application Sockeye, being developed through the Nuclear En- ergy Advanced Modeling and Simulation program is concurrently undergoing code development and verification and validation activities. The Single Primary Heat Extraction and Removal Emu- lator (SPHERE) facility at Idaho National Laboratory has performed numerous sodium heat pipe experiments, one of which is modeled in this work to support Sockeye validation. Heat pipe data produced by the SPHERE facility was found to be very useful not only for Sockeye validation but also as a useful resource for model development. A Sockeye simulation performed using Sockeye’s conduction model modeled an experimental run at SPHERE and was found to be able to accurately reproduce the steady temperature distribution but revealed discrepancies in transient behavior after shutoff of the power. This effort also revealed a number of challenges of modeling an experimental setup such as SPHERE, for example, the importance of details in the setup such as possible gaps in insulation. Also, continuous collaboration between Sockeye developers and SPHERE experimentalists has been established, which allows future experiments to be refined in response to questions and challenges found by simulating earlier experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DSMC Shock Simulation of Saturn Entry Probe Conditions

This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 kms shock waves are obtained at 0.2 and 0.1 Torr respectively and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.

Modeling↗

DSMC Shock Simulation of Saturn Entry Probe Conditions

This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at the NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 km/s shock waves are obtained at 0.2 and 0.1 Torr, respectively, and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.

Higdon, Kyle J.↗