Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 2”

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,297 records · Page 72

Model Rotor Hover Performance at Low Reynolds Number

Hover performance data from four key experiments has been analyzed in detail to shed some light on model rotor hover performance at low Reynolds number. Each experiment used the simplest blade geometry. The blades were constant chord and untwisted. Three experiments used blades with the NACA 0012 airfoil from root to tip. The NACA 0015 was used in the earliest test. The four experiments provide data spanning a Reynolds number range of 136,500 to 548,700. The specific objective of this report is to ask and answer two questions: 1. Does blade aspect ratio influence hover performance or is rotor solidity the fundamental rotor geometry parameter for practical engineering purposes? ANSWER: Rotor solidity is the fundamental rotor geometry parameter for practical engineering purposes. Any effect of blade aspect ratio appears to be such a secondary variable that its effect lies within the range of experimental error. 2. Is Reynolds number a significant factor in scaling up hover performance to full-scale rotor performance? The answer is twofold. ANSWER: (a) Reynolds number effects on the increase of power with thrust do not appear to be a significant factor for practical engineering purposes, and (b) Reynolds number effects on minimum profile power at or very near zero rotor thrust could not be clearly established primarily because the low torque levels could not be accurately measured with the test equipment used. This has led to significant data scatter. A number of other observations can be made based on the analysis provided herein. For instance: 1. The test matrices used in the four key references contained far too few data points. A collective pitch variation of four or five data points is insufficient to establish experimental accuracy and data repeatability. 2. A common property of the power-versus-thrust (raised to the 3/2 exponent) graphs was that this curve was linear below the onset of blade stall. 3. The blade-to-blade interference at or near zero thrust may, in fact, be creating a turbulent flow field such that the effective Reynolds number at a blade element is considerably greater than what theories using two-dimensional (2D) airfoil properties at a blade element would calculate. 4. Definitive experiments answering the two key questions have yet to be made.

Model Rotor Hover Performance↗

Evaluation of single and multiple Doppler lidar techniques to measure complex flow during the XPIA field campaign

Abstract. Accurate three-dimensional information of wind flow fields can be an important tool in not only visualizing complex flow but also understanding the underlying physical processes and improving flow modeling. However, a thorough analysis of the measurement uncertainties is required to properly interpret results. The XPIA (eXperimental Planetary boundary layer Instrumentation Assessment) field campaign conducted at the Boulder Atmospheric Observatory (BAO) in Erie, CO, from 2 March to 31 May 2015 brought together a large suite of in situ and remote sensing measurement platforms to evaluate complex flow measurement strategies. In this paper, measurement uncertainties for different single and multi-Doppler strategies using simple scan geometries (conical, vertical plane and staring) are investigated. The tradeoffs (such as time–space resolution vs. spatial coverage) among the different measurement techniques are evaluated using co-located measurements made near the BAO tower. Sensitivity of the single-/multi-Doppler measurement uncertainties to averaging period are investigated using the sonic anemometers installed on the BAO tower as the standard reference. Finally, the radiometer measurements are used to partition the measurement periods as a function of atmospheric stability to determine their effect on measurement uncertainty. It was found that with an increase in spatial coverage and measurement complexity, the uncertainty in the wind measurement also increased. For multi-Doppler techniques, the increase in uncertainty for temporally uncoordinated measurements is possibly due to requiring additional assumptions of stationarity along with horizontal homogeneity and less representative line-of-sight velocity statistics. It was also found that wind speed measurement uncertainty was lower during stable conditions compared to unstable conditions.

17 WIND ENERGY↗

Environment modelling in near Earth space: Preliminary LDEF results

Hypervelocity impacts by space debris cause not only local cratering or penetrations, but also cause large areas of damage in coated, painted or laminated surfaces. Features examined in these analyses display interesting morphological characteristics, commonly exhibiting a concentric ringed appearance. Virtually all features greater than 0.2 mm in diameter possess a spall zone in which all of the paint was removed from the aluminum surface. These spall zones vary in size from approximately 2 - 5 crater diameters. The actual craters in the aluminum substrate vary from central pits without raised rims, to morphologies more typical of craters formed in aluminum under hypervelocity laboratory conditions for the larger features. Most features also possess what is referred to as a 'shock zone' as well. These zones vary in size from approximately 1 - 20 crater diameters. In most cases, only the outer-most layer of paint was affected by this impact related phenomenon. Several impacts possess ridge-like structures encircling the area in which this outer-most paint layer was removed. In many ways, such features resemble the lunar impact basins, but on an extremely reduced scale. Overall, there were no noticeable penetrations, bulges or spallation features on the backside of the tray. On Row 12, approximately 85 degrees from the leading edge (RAM direction), there was approximately one impact per 15 cm(exp 2). On the trailing edge, there was approximately one impact per 72 cm(exp 2). Currently, craters on four aluminum experiment trays from Bay E09, directly on the leading edge are being measured and analyzed. Preliminary results have produced more than 2200 craters on approximately 1500 cm(exp 2) - or approximately 1 impact per 0.7 cm(exp 2).

Coombs, C. R.↗

Model‐Based Approach to Selecting Pegfilgrastim Dose for Pharmacokinetic and Pharmacodynamic Similarity Studies in Biosimilar Development

This study applied modeling and simulation (M&S) approaches to evaluate the sensitivity of pegfilgrastim pharmacokinetics (PKs) and pharmacodynamics (PDs) to changes in dose amount, and linear or nonlinear clearance (CL) over pegfilgrastim subcutaneous dose of 2–6 mg. A previously published model was adapted to better describe pegfilgrastim PK and PD data in healthy subjects and used in simulation. Nonlinear CL accounts for 98% and 77%, respectively, of the total CL of pegfilgrastim at 2 and 6 mg. The sensitivity analyses showed: (i) PK of 2 and 6 mg doses are similarly sensitive to detect differences for a 5% change in dose; (ii) PK of 2 mg dose is more sensitive to changes in receptor binding affinity, a model parameter for nonlinear CL, and a product quality attribute characterized with orthogonal methods as part of demonstrating analytical similarity between products; (iii) PK of approved 6 mg dose is more sensitive to changes in linear CL, which has not been associated with any specific product quality attributes, and (iv) the PDs are not sensitive to changes in linear or nonlinear CL. Taken together, our analyses support that the approved pegfilgrastim dose of 6 mg is appropriate for detecting differences between a biosimilar and the reference products in pegfilgrastim PK and PD similarity studies. The described M&S approaches can be adopted to support dose selection for biosimilars with nonlinear PK and complex PK–PD interplay.

Pharmacology & Pharmacy↗

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Angular decay coefficients of J/ψ mesons at forward rapidity from p+p collisions at √s = 510 GeV

In this paper, we report the first measurement of the full angular distribution for inclusive J/ψ → μ + μ - decays in p+p collisions at √s = 510 GeV. The measurements are made for J/ψ transverse momentum 2 < p T < 10 GeV/c and rapidity 1.2 < y < 2.2 in the Helicity, Collins-Soper, and Gottfried-Jackson reference frames. In all frames the polar coefficient λ θ is strongly negative at low p T and becomes close to zero at high p T , while the azimuthal coefficient λ Φ is close to zero at low p T , and becomes slightly negative at higher p T . The frame-independent coefficient λ ~ is strongly negative at all p T in all frames. Finally, the data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗

Impact of Zr Content on Radionuclide Release Rates in Metal Waste Forms

Reference alloy waste forms (RAW) underwent hybrid corrosion/immersion testing to parameterize the ANL analytical oxidative-dissolution model to enable the calculation of fractional release rates. The waste forms tested at SRNL were variations of the RAW-3 formulation that uses HT9 as the main alloy component, and are meant to enable evaluation of the impact of Zr:U ratio on the release rates of actinides and Tc-99. The test solutions were air saturated alkaline and acidic brines, representing potential repositories with those conditions. The testing approach consisted of 4 major steps; 1) bare surface corrosion measurements, 2) hybrid potentiostatic hold/exposure measurements, 3) measurement of radionuclide concentrations and relations to anodic current from potentiostatic holds, and 4) identification of corroding phases using SEM/EDS of electrodes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

HFTS-1 Natural Joint Data and Engineering Summary

This report provides an analysis of engineering and geologic data collected from the Hydraulic Fracturing Test Site #1 (HTSF-1) project in the southern Midland Basin, Reagan County, Texas. The site is being studied as part of the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative at the National Energy Technology Laboratory (NETL). The data collected is intended to provide a basis of understanding of the site and to construct simulation models. This report provides a summary analysis of the engineering and geologic data collected from the project to provide a basis of understanding of the site and to construct simulation models. The report examines various possible correlations in engineering properties based on natural fracture data from the program, which consisted of 11 horizontal wells, one vertical well, and one slant well. The horizontal wells are in two horizons: 1) Upper and 2) Middle Wolfcamp formations. Program data include: fracture frequency and fracture orientation data from four core runs in a slant well; fracture spacing and orientation from a vertical pilot well; laboratory triaxial testing and mineralogical determinations; and porosity results from magnetic resonance analyses from various wells, together with observations based on the data collection. In addition, available references were reviewed on the site for additional insights. As the focus of the report is on the natural system, hydraulic fracture data from the site were not examined in detail in this report. Data variability is the chief observation in examination of the database. Fracture frequency in the slant well can range from sections with values as high as five fractures per ft to sections up to 100+ ft in length with no natural observed fractures. Fracture spacing across is typically less than 10 ft, but can range up to hundreds of feet. Apparent fracturing shows the trends in two predominate orientations, E-SW and WNW-ESE, but minor variations exist. The rock units vary across the site from siliceous mudstones to calcareous mudstones, showing a general layering with depth. The laboratory properties such as strength and modulus show no apparent trend with depth, but appear to correlate with rock mineralogy with high strength and modulus values where calcium content is high. In addition, an attempt to examine variability and mineralogy on a larger scale was made using a color-coded system based on gamma ray measurements. A staged colored approach was adopted, presuming that lower gamma ray values indicate higher value of calcium content (blue scale) and that higher gamma ray values indicate higher clay mineral content (orange scale). As provided in report appendices, the system correlated well with visual examination of the slant core and the petrofabric analyses of the vertical pilot well. The results showed large variability in mineral content along the horizontal plane across the site. The change in mineral content was also rapid, on a scale less than that of the average hydraulic fracture stage length of about 180 ft.

42 ENGINEERING↗

MELCOR Accident Progression and Source Term Demonstration Calculations for a FHR

MELCOR is an integrated thermal hydraulics, accident progression, and source term code for reactor safety analysis that has been developed at Sandia National Laboratories for the United States Nuclear Regulatory Commission (NRC) since the early 1980s. Though MELCOR originated as a light water reactor (LWR) code, development and modernization efforts have expanded its application scope to include non-LWR reactor concepts. Current MELCOR development efforts include providing the NRC with the analytical capabilities to support regulatory readiness for licensing non-LWR technologies under Strategy 2 of the NRC's near- term Implementation Action Plans. Beginning with the Next Generation Nuclear Project (NGNP), MELCOR has undergone a range of enhancements to provide analytical capabilities for modeling the spectrum of advanced non-LWR concepts. This report describes the generic plant model developed to demonstrate MELCOR capabilities to perform fluoride-salt-cooled high-temperature reactor (FHR) safety evaluations. The generic plant model is based on publicly-available FHR design information. For plant aspects (e.g., reactor building leak rate and details of the cover-gas system) that are not described in the FHR references, the analysts made assumptions needed to construct a MELCOR full-plant model. The FHR model uses a TRi-structural ISOtropic (TRISO)-particle fuel pebble-bed reactor with a primary system rejecting heat to two coiled tube air heat ex changers. Three passive direct reactor auxiliary cooling systems provide heat removal to supplement or replace the emergency secondary system heat removal during accident conditions. Surrounding the reactor vessel is a low volume reactor cavity that insulates the reactor with fire bricks and thick concrete walls. A refractory reactor liner system provides water cooling to reduce the concrete wall temperatures. Example calculations are performed to show the plant response and MELCOR capabilities to characterize a range of accident conditions. The accidents selected for evaluation consider a range of degraded and failed modes of operation for key safety functions providing reactivity control, the primary system decay heat removal and also a piping leak of the line to the coolant drain tank.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Zero Power Warming (ZPW) Chamber Prototype Measurements, Utqiagvik (Barrow), Alaska, 2016

Data were collected during one season of prototyping associated with the development of a passive warming technology. An experimental chamber, the Zero Power Warming (ZPW) chamber, was fitted with apparatus to modulate venting of a field enclosure and enhance elevation of air temperature by solar radiation. The ZPW chamber was compared with a control chamber (Control) and an ambient open air plot (Ambient). The control chamber was identical to the ZPW chamber but lacked the apparatus necessary to modulate venting, the chamber vents in the control chamber were fixed open for the majority of the trial period. The three plots were located over Carex aquatilis growing in an area of moderately degraded permafrost. Chambers were placed on the same footprints that were used for a similar exercise in 2015 (no data) and therefore those plots had experienced some thaw and degradation prior to 2016. The following data were collected for 80 days at 1 minute intervals from within two chambers and an ambient plot: solar input, chamber venting, air temperature, relative humidity, soil temperature (at 5, 10 and 15 cm), soil moisture, downward and upward NIR. CO2 concentration within the chambers and of the ambient air plot was collected at 1 minute intervals over a 62 hour period in the middle of the growing season. This dataset contains six csv files and one pdf file.The documentation may refer to a ZPW-metadata.xlsx file. This file was split into two csv files: ZPW-metadata.csv and ZPW-metadata-instrumentation.csv.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

F2 peak electron density at Millstone Hill and Hobart: Comparsion of theory and measurement at solar maximum

This paper compares the observed behavior of the (F2) layer of the ionosphere at Millstone Hill and Hobart with calculations from the field line interhemispheric plasma (FLIP) model for solar maximum, solstice conditions in 1990. During the study period the daily F(sub 10.7) index varied by more than a factor of 2 (123 to 280), but the 81-day mean F(sub 10.7) (F(sub 10.7 A)) was almost constant near 190. Calculations were performed with and without the effects of vibrationally excited N2 (N(sup *)(sub 2) which affects the loss rate of atomic oxygen ions. In the case without N(sup *)(sub 2) there is generally good agreement between the model and measurement for the daytime, peak density of the F region (NmF2). Both the model and the measurement show a strong seasonal anomaly with the winter noon densities a factor of 3 to 4 greater than the summer noon densities at Millstone Hill and a factor of 2 greater at Hobart. The seasonal anomaly in the model is caused by changes in the neutral composition as given by the mass spectrometer and incoherent scatter (MSIS) 86 neutral density model. There is generally little or no increase in the observed noon NmF2 as a function of daily F(sub 10.7) except at Millstone Hill in winter. In contrast to the generally good agreement between model and data at noon, the model badly underestimates the density at night at Millstone Hill at all seasons. At Hobart the model reproduces the nighttime density variations well in both winter and summer. The international reference ionosphere (IRI) model generally provides a good representation of the average behavior of noon NmF2 and hmF2 but because the data show a lot of day-to-day variability, there are often large differences. The FLIP model is able to reproduce this variability when hmF2 is specified. The IRI model peak densities are better than the FLIP densities at night, but the IRI model does not represent the Millstone Hill summer data very well at night in 1990.

Richards, P. G.↗

Cross Calibration of TOMS, SBUV/2 and Sciamachy Radiances from Ground Observations

Verification of a stratospheric ozone recovery remains a high priority for environmental research and policy definition. Models predict an ozone recovery at a much lower rate than the measured depletion rate observed to date. Therefore improved precision of the satellite and ground ozone observing systems are required over the long term to verify recovery. We have shown that validation of radiances is the most effective means for correcting absolute accuracy and long term drifts of backscatter type satellite measurements. This method by-passes the algorithms used for both satellite and ground based measurements which are normally used to validate and correct the satellite data. Validation of radiances will also improve all higher level data products derived from the satellite observations. Backscatter algorithms suffer from several errors such as unrepresentative a-priori data and air mass factor corrections. Radiance comparisons employ forward models but are inherently more accurate and than inverse (retrieval) algorithms. A new method for satellite validation is planned which will compliment measurements from the existing ground-based networks. This method will employ very accurate comparisons between ground based zenith sky radiances and satellite nadir radiances. These comparisons will rely heavily on the experience derived from the Shuttle SBUV (SSBUV) program which provided a reference standard of radiance measurements for SBUV/2, TOMS, and GOME. This new measurement program, called "Skyrad", employs two well established capabilities at the Goddard Space Flight Center, 1) the SSBUV calibration facilities and 2) the radiative transfer codes used for the TOMS and SBUV/2 algorithms and their subsequent refinements. Radiative transfer calculations show that ground based zenith sky and satellite nadir backscatter ultraviolet comparisons can be made very accurately under certain viewing conditions. The Skyrad instruments (SSBUV, Brewer spectrophotometers, and possibly others) will be calibrated and maintained to a precision of a few tenths of a percent. Skyrad data will then enable long term calibration of upcoming satellite instruments such as QuickTOMS. SBUV/2s and SCIAMACHY with a high degree of precision. This technique can be further employed to monitor the performance of future instruments such as GOME-2, OMI, and OMPS. Initial ground observations taken from Goddard Space Flight Center compared with radiative transfer calculations has indicated the feasibility of this method.

Hillsenrath, Ernest↗

International Energy Agency Ocean Energy Systems Task 10 Wave Energy Converter Modeling Verification and Validation: Preprint

This is the first joint reference paper for the Ocean Energy Systems (OES) Task 10 Wave Energy Converter modeling verification and validation group. The group is established under the OES Energy Technology Network program under the International Energy Agency. OES was founded in 2001 and Task 10 was proposed by Bob Thresher (National Renewable Energy Laboratory) in 2015 and approved by the OES Executive Committee EXCO in 2016. The kickoff workshop took place in September 2016, wherein the initial baseline task was defined. Experience from similar offshore wind validation/verification projects (OC3-OC5 conducted within the International Energy Agency Wind Task 30) [1], [2] showed that a simple test case would help the initial cooperation to present results in a comparable way. A heaving sphere was chosen as the first test case. The team of project participants simulated different numerical experiments, such as heave decay tests and regular and irregular wave cases. The simulation results are presented and discussed in this paper.

16 TIDAL AND WAVE POWER↗

On-the-fly generation of differential resonance scattering probability distribution functions for Monte Carlo codes

Current Monte Carlo codes use one of three models to model neutron scattering in the epithermal energy range: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S({alpha},{beta}) model, depending on the neutron energy and the specific Monte Carlo code. The free gas scattering model assumes the scattering cross section is constant over the neutron energy range, which is usually a good approximation for light nuclei, but not for heavy nuclei where the scattering cross section may have several resonances in the epithermal region. Several researchers in the field have shown that using the free gas scattering model in the vicinity of the resonances in the lower epithermal range can under-predict resonance absorption due to the up-scattering phenomenon. Existing methods all involve performing the collision analysis in the center-of-mass frame, followed by a conversion back to the laboratory frame. In this paper, we will present a new sampling methodology that (1) accounts for the energy-dependent scattering cross sections in the collision analysis and (2) acts in the laboratory frame, avoiding the conversion to the center-of-mass frame. The energy dependence of the scattering cross section was modeled with even-ordered polynomials to approximate the scattering cross section in Blackshaw's equations for the moments of the differential scattering PDFs. These moments were used to sample the outgoing neutron speed and angle in the laboratory frame on-the-fly during the random walk of the neutron. Results for criticality studies on fuel pin and fuel assembly calculations using these methods showed very close comparison to results using the reference Doppler-broadened rejection correction (DBRC) scheme. (authors)

97 MATHEMATICS AND COMPUTING↗

Slip Tendency Analysis of Fracture Networks to Determine Suitability of Candidate Testbeds for the EGS Collab Hydroshear Experiment

Experiment 2 of the EGS Collab project is aimed at testing stimulation by hydro-shearing of existing natural fractures, versus Experiment 1, which is focused on hydraulic fracturing a rock mass to enhance permeability. The main criterion for the testbed selection in Experiment 2 is the presence of an interconnected network of fractures, at least 10 meters in extent, and with orientations such that shear slip can be induced at injection pressures less than the minimum horizontal stress (Shmin). The feasibility analysis for this experiment requires a well-constrained stress state along with well-characterized fracture networks. The fracture systems at two candidate locations in the Sanford Underground Research Facility (SURF), the 4850 Level (number refers to depth below ground surface in feet) and the 4100 Level, have been characterized to different extents. The 4850 Level has a well-characterized Discrete Fracture Network (DFN) from borehole (drilled for an experimental test bed) and drift observations. In contrast, the 4100 level has a fracture network characterized only by observations from the drift wall, which provides little constraint on fracture extents. This paper will present assessments for the interpreted slip potential for the natural fractures at the two locations. Data uncertainties are addressed by performing a probabilistic analysis that takes into account the uncertainty in the stress state, uncertainty in the fracture properties, and preliminary borehole locations. There is strong evidence of shear stimulation of some natural fractures in Experiment 1 on the 4850 level. We use examples of natural fractures strongly linked to shear stimulation in Experiment 1 from multiple monitoring indicators to test the consistency of the stress model and guide the uncertainty interpretation. This analysis is intended to guide the site selection process for Experiment 2 by highlighting the fracture orientations that are likely to be shear-stimulated in a majority of modeled realizations.

15 GEOTHERMAL ENERGY↗

Slip Tendency Analysis of Fracture Networks to Determine Suitability of Candidate Testbeds for the EGS Collab Hydroshear Experiment

Experiment 2 of the EGS Collab project is aimed at testing stimulation by hydro-shearing of existing natural fractures, versus Experiment 1, which is focused on hydraulic fracturing a rock mass to enhance permeability. The main criterion for the testbed selection in Experiment 2 is the presence of an interconnected network of fractures, at least 10 meters in extent, and with orientations such that shear slip can be induced at injection pressures less than the minimum horizontal stress (Shmin). The feasibility analysis for this experiment requires a well-constrained stress state along with well-characterized fracture networks. The fracture systems at two candidate locations in the Sanford Underground Research Facility (SURF), the 4850 Level (number refers to depth below ground surface in feet) and the 4100 Level, have been characterized to different extents. The 4850 Level has a well-characterized Discrete Fracture Network (DFN) from borehole (drilled for an experimental test bed) and drift observations. In contrast, the 4100 level has a fracture network characterized only by observations from the drift wall, which provides little constraint on fracture extents. This paper will present assessments for the interpreted slip potential for the natural fractures at the two locations. Data uncertainties are addressed by performing a probabilistic analysis that takes into account the uncertainty in the stress state, uncertainty in the fracture properties, and preliminary borehole locations. There is strong evidence of shear stimulation of some natural fractures in Experiment 1 on the 4850 level. We use examples of natural fractures strongly linked to shear stimulation in Experiment 1 from multiple monitoring indicators to test the consistency of the stress model and guide the uncertainty interpretation. This analysis is intended to guide the site selection process for Experiment 2 by highlighting the fracture orientations that are likely to be shear-stimulated in a majority of modeled realizations.

40 EE - Geothermal Technologies Office (EE-4G)↗

Global Reference Atmospheric Model (GRAM) Advancements and Additions

Introduction:The Global Reference Atmospheric Model (GRAM) is one of the most widely used engineering models of the atmosphere. GRAM development and maintenance has been led by NASA Marshall Space Flight Center (MSFC). The NASA Science Mission Directorate (SMD) has provided funding support to upgrade the GRAMs since Fiscal Year 2018. NASA Langley Research Center has been working with MSFC on the upgrades.This presentation will provide details regarding the upgrades that have been made to the existing GRAMs, the development of new GRAMs, as well as the ongoing objectives, tasks, and milestones re-lated to the GRAM upgrades funded by NASA SMD. GRAM: The GRAMs are engineering-oriented atmospheric models that estimate mean values and statistical variations of the atmospheric properties for numerous planetary destinations.They provide mean values and variability for any point in the atmosphere as well as seasonal, geographic, and altitude variations. GRAM outputs include atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path.Theyhave been widely used by the engineering community because of their ability to create realistic dispersions. GRAMs have been integrated into high fidelity flight dynamic simulations of launch, entry, descent and landing (EDL), aerobraking and aerocapture. MSFC has been developing and updating GRAMs since 1974; GRAMs are currently available for Earth, Mars, Venus, Neptune, and Titan. GRAM Upgrade Status: Code Moderization.The planetary GRAMs are being rearchitected from Fortran to a common object-ori-entedC++ framework called the GRAM Suite. This new architecture creates a common GRAM library of data models and utilities. The first C++ releases of the existing planetary GRAMs (Mars, Venus, Neptune, and Titan-GRAM) in the GRAM Suite are straight conversions from the latest Fortran version. Model Upgrades.The focus of the model upgrade task is to improve the atmosphere models in the existing GRAMs and to establish a foundation for developing GRAMs for additional destinations. The GRAM ephemerishas been upgraded to the NASA Navigation and Ancillary Information Facility (NAIF) SPICE toolkit (version N0066). The calculation of the speed of sound has also been improved in the GRAMs. The GRAM team has received updated Mars General Circulation Model (MGCM) datasets from NASA Ames Research Center. Mars Global Ionosphere-Thermosphere Model (M-GITM) data is being obtained toreplace the Mars Thermospheric General Circulation Model (MTGCM) data in legacy Mars-GRAM. M-GITM and updated MGCM data will be incorporated into a future GRAM Suite release.Twoprojects that will improve the atmospheric model data in the GRAMs have been funded by the GRAM team since Fiscal Year 2020. Sanjay Limaye and Patrick Fry at the University of Wisconsin are reanalyzing the Venus Express radio occultation observations and analyzing the Akatsukiradio occultation observations.This will lead tothe calculation of number density, temperature, and pressure profiles for the 40-90 km altitude range. Kunio Sayanagi, Justin Garland, and Ryan McCabeat Hampton University are developing empirical global models for Venus,Jupiter, Saturn, Uranus, Neptune, and Titan that incorporates the latest data available for each of these planetary destinations. Upgraded GRAM Releases. GRAM Suite Version 1.0 was released in May 2020 and contains the rearchitected Neptune-GRAM, including the common GRAM framework and planet–specific code. GRAM Suite Version 1.1 was released in September 2020 and addedthe rearchitected Titan-GRAM to the GRAM Suite. A User Guide and Programmer’s Manualarereleased with all GRAMs. The rearchitected Mars and Venus-GRAMs will be released in upcomingversions of the GRAM Suite. New GRAM Releases. New GRAMs have been de-veloped for Uranus and Jupiter. Uranus-GRAM is based on an individual profile generated by Gary Allen (ARC) from Voyager 2 occultation data and will be released in GRAM Suite Version 1.2. Jupiter-GRAM is based on individual profile produced from Al Seiff’s Ju-piter model[1].Jupiter-GRAM will be released in GRAM Suite Version 1.3. Saturn-GRAM is currently under development and will be released in a future version of the GRAM Suite. Conclusions: GRAMs are vital and frequently used toolsets. Releases of the GRAM Suite, upgrades of the existing planetary GRAMs, and development of new planetary GRAMs are ongoing. NASA SMD funding has been essential to addressing current limitations and accomplishing GRAM developmental goals. Continua-tion ofSMD fundingwill ensure the development, up-grades, and maintenance of the GRAMs. References: [1] Seiff, A., et al. (1998) JGR, 103, 22,857-22,889. Acknowledgments: The authors gratefully acknowledge support from the NASA SMD

atmospheric models↗