Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHARP”

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 307 records · Page 17

Artificial intelligence for multi-mission planetary operations

A brief introduction is given to an automated system called the Spacecraft Health Automated Reasoning Prototype (SHARP). SHARP is designed to demonstrate automated health and status analysis for multi-mission spacecraft and ground data systems operations. The SHARP system combines conventional computer science methodologies with artificial intelligence techniques to produce an effective method for detecting and analyzing potential spacecraft and ground systems problems. The system performs real-time analysis of spacecraft and other related telemetry, and is also capable of examining data in historical context. Telecommunications link analysis of the Voyager II spacecraft is the initial focus for evaluation of the prototype in a real-time operations setting during the Voyager spacecraft encounter with Neptune in August, 1989. The preliminary results of the SHARP project and plans for future application of the technology are discussed.

Atkinson, David J.↗

Aerothermal study in an axial compression corner with shock impingement at Mach 6

An experimental study of an axial-compression interior-corner region with shock impingement was conducted in the NASA Langley 20-inch Mach 6 Tunnel to obtain aerodynamic heating rate distributions. The model, a generic hypersonic engine inlet, had a sharp leading edge forebody, cowl, and interchangeable struts with sharp and blunt leading edges. Experiments were conducted at free-stream unit Reynolds numbers of 1.68 x 10 to the 6th and 3.35 x 10 to the 6th per foot for the blunt and sharp-strut configurations, respectively. For the sharp-strut configuration, with and without trips, the peak heating occurred on the cowl, near the compression corner, downstream of the forebody oblique shock impingement region. The peak heating levels, for both the tripped and untripped boundary-layer flows were about 45 and 50 times the reference laminar heating levels, respectively. The peak heating for the blunt leading edge strut model, which was 30 times the reference value, occurred on the cowl, away from the corner, where the bow shock from the strut leading edge intersected the cowl surface.

Venkateswaran, S.↗

Overview of the relevant CFD work at Thiokol Corporation

An in-house developed proprietary advanced computational fluid dynamics code called SHARP (Trademark) is a primary tool for many flow simulations and design analyses. The SHARP code is a time dependent, two dimensional (2-D) axisymmetric numerical solution technique for the compressible Navier-Stokes equations. The solution technique in SHARP uses a vectorizable implicit, second order accurate in time and space, finite volume scheme based on an upwind flux-difference splitting of a Roe-type approximated Riemann solver, Van Leer's flux vector splitting, and a fourth order artificial dissipation scheme with a preconditioning to accelerate the flow solution. Turbulence is simulated by an algebraic model, and ultimately the kappa-epsilon model. Some other capabilities of the code are 2-D two-phase Lagrangian particle tracking and cell blockages. Extensive development and testing has been conducted on the 3-D version of the code with flow, combustion, and turbulence interactions. The emphasis here is on the specific applications of SHARP in Solid Rocket Motor design. Information is given in viewgraph form.

Chwalowski, Pawel↗

Gale Crater - Why are We There and What do We Hope to Learn?

The Mars Science Laboratory Rover Curiosity is commencing a two-year investigation of Gale crater and Mt. Sharp, the craters prominent central mound. Gale is a 155 km, late Noachian/early Hesperian impact crater located near the dichotomy boundary separating the southern highlands from the northern plains. The central mound is composed of layered sedimentary rock, with upper and lower mound units separated by a prominent erosional unconformity. The lower mound is of particular interest, as it contains secondary minerals indicative of a striking shift from water-rich to water-poor conditions on early Mars. A key unknown in the history of Gale is the relationship between the sedimentary units in the mound and sedimentary sequences in the surrounding region. We employed orbital remote sensing data to determine if areas within a 1,000 km radius of Gale match the characteristics of sedimentary units in Mt. Sharp. Regions of interest were defined based on: the mound s inferred age, altitude range, and THEMIS nighttime brightness (a proxy for thermal inertia). Using orbital CTX, MOC and HiRISE images we examined all areas within our regions of interest for analogous geomorphic units in the same altitude ranges as the corresponding units in Mt. Sharp. The results are consistent with the hypothesis that sedimentary units in both the upper and lower sections of the Gale mound are related to nearby regional units located along the dichotomy boundary. This relationship supports an inferred geologic history that includes several episodes of widespread sedimentary deposition and erosion in the martian mid-latitudes. In this model Mt. Sharp is the remnant of regional sedimentary deposits that partially or completely filled the crater, became lithified, and were subsequently deeply eroded. Key questions that will be addressed by Curiosity include the compositions of the sediments, the modes of deposition, the mechanisms of lithification, and the nature of the erosion.

Allen, Carlton C.↗

Calculating and Mitigating the Risk of a Cut Glove to a Space Walking Astronaut

One of the high risk operations on the International Space Station (ISS) is conducting a space walk, or an Extra Vehicular Activity (EVA). Threats to the space walking crew include airlock failures, space suit failures, and strikes from micro ]meteoroids and orbital debris (MM/OD). There are risks of becoming untethered from the space station, being pinched between the robotic arm and a piece of equipment, tearing your suit on a sharp edge, and other human errors that can be catastrophic. For decades NASA identified and tried to control sharp edges on external structure and equipment by design; however a new and unexpected source of sharp edges has since become apparent. Until recently, one of the underappreciated environmental risks was damage to EVA gloves during a spacewalk. The ISS has some elements which have been flying in the environment of space for over 14 years. It has and continues to be bombarded with MM/OD strikes that have created small, sharp craters all over the structure, including the dedicated EVA handrails and surrounding structure. These craters are capable of cutting through several layers of the EVA gloves. Starting in 2006, five EVA crewmembers reported cuts in their gloves so large they rendered the gloves unusable and in some cases cut the spacewalk short for the safety of the crew. This new hazard took engineers and managers by surprise. NASA has set out to mitigate this risk to safety and operations by redesigning the spacesuit gloves to be more resilient and designing a clamp to isolate MM/OD strikes on handrails, and is considering the necessity of an additional tool to repair strikes on non ]handrail surfaces (such as a file). This paper will address how the ISS Risk Team quantified an estimate of the MM/OD damage to the ISS, and the resulting likelihood of sustaining a cut glove in order to measure the effectiveness of the solutions being investigated to mitigate this risk to the mission and crew.

Castillo, Theresa↗

Advances and Limitations of Atmospheric Boundary Layer Observations with GPS Occultation over Southeast Pacific Ocean

The typical atmospheric boundary layer (ABL) over the southeast (SE) Pacific Ocean is featured with a strong temperature inversion and a sharp moisture gradient across the ABL top. The strong moisture and temperature gradients result in a sharp refractivity gradient that can be precisely detected by the Global Positioning System (GPS) radio occultation (RO) measurements. In this paper, the Constellation Observing System for Meteorology, Ionosphere & Climate (COSMIC) GPS RO soundings, radiosondes and the high-resolution ECMWF analysis over the SE Pacific are analyzed. COSMIC RO is able to detect a wide range of ABL height variations (1-2 kilometer) as observed from the radiosondes. However, the ECMWF analysis systematically underestimates the ABL heights. The sharp refractivity gradient at the ABL top frequently exceeds the critical refraction (e.g., −157 N-unit per kilometer) and becomes the so-called ducting condition, which results in a systematic RO refractivity bias (or called N-bias) inside the ABL. Simulation study based on radiosonde profiles reveals the magnitudes of the N-biases are vertical resolution dependent. The N-bias is also the primary cause of the systematically smaller refractivity gradient (rarely exceeding −110 N-unit per kilometer) at the ABL top from RO measurement. However, the N-bias seems not affect the ABL height detection. Instead, the very large RO bending angle and the sharp refractivity gradient due to ducting allow reliable detection of the ABL height from GPS RO. The seasonal mean climatology of ABL heights derived from a nine-month composite of COSMIC RO soundings over the SE Pacific reveals significant differences from the ECMWF analysis. Both show an increase of ABL height from the shallow stratocumulus near the coast to a much higher trade wind inversion further off the coast. However, COSMIC RO shows an overall deeper ABL and reveals different locations of the minimum and maximum ABL heights as compared to the ECMWF analysis. At low latitudes, despite the decreasing number of COSMIC RO soundings and the lower percentage of soundings that penetrate into the lowest 500-m above the mean-sea-level, there are small sampling errors in the mean ABL height climatology. The difference of ABL height climatology between COSMIC RO and ECMWF analysis over SE Pacific is significant and requires further studies.

COSMIC RO↗

Mineralogy of Mudstone at Gale Crater, Mars: Evidence for Dynamic Lacustrine Environments

The Mars Science Laboratory Curiosity rover landed in Gale crater in August 2012 to assess the habitability of sedimentary deposits that show orbital evidence for diverse ancient aqueous environments. Gale crater contains a 5 km high mound of layered sedimentary rocks in its center, informally named Mount Sharp. The lowermost rocks of Mount Sharp contain minerals that are consistent with a dramatic climate change during Mars' early history. During the rover's traverse across the Gale crater plains to the base of Mount Sharp, Curiosity discovered sedimentary rocks consistent with a fluviolacustrine sequence. Curiosity studied ancient lacustrine deposits at Yellowknife Bay on the plains of Gale crater and continues to study ancient lacustrine deposits in the Murray formation, the lowermost unit of Mount Sharp. These investigations include drilling into the mudstone and delivering the sieved less than 150 micrometers fraction to the CheMin XRD/XRF instrument inside the rover. Rietveld refinement of XRD patterns measured by CheMin generates mineral abundances with a detection limit of 1-2 wt.% and refined unit-cell parameters of minerals present in abundances greater than approximately 5 wt.%. FULLPAT analyses of CheMin XRD patterns provide the abundance of X-ray amorphous materials and constrain the identity of these phases (e.g., opal-A vs. opal-CT). At the time of writing, CheMin has analyzed 14 samples, seven of which were drilled from lacustrine deposits. The mineralogy from CheMin, combined with in-situ geochemical measurements and sedimentological observations, suggest an evolution in the lake waters through time, including changes in pH and salinity and transitions between oxic and anoxic conditions. In addition to a geochemically dynamic lake environment, the igneous minerals discovered in the lake sediments indicate changes in source region through time, with input from mafic and silicic igneous sources. The Murray formation is predominantly comprised of lacustrine mudstone and is 150-200 m thick, suggesting long history of lake environments in Gale crater. Curiosity has traversed through the lowermost approximately 30 m of the Murray formation, and each additional sample provides clues about the climate on early Mars.

Rampe, E. B.↗

Shockwave/Boundary-Layer Interaction Studies Performed in the NASA Langley 20-Inch Mach 6 Air Tunnel

This paper highlights results from a collaborative study performed by The University of Tennessee Space Institute (UTSI) and NASA Langley Research Center on the Shockwave/Boundary-Layer Interaction (SWBLI) generated by a cylindrical protuberance on a flat plate in a Mach 6 flow. The study was performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center and consisted of two separate entries. In the first entry, simultaneous high-speed schlieren and high-speed pressure-sensitive paint (PSP) imaging – which was performed for the first time in the 20-Inch Mach 6 facility at NASA Langley – were performed as well as simultaneous high-speed schlieren and oil-flow imaging. In the second entry, the model configuration was modified to increase the size of the interaction region. High-speed schlieren and infrared thermography (IR) surface imaging were performed in this second entry. The goal of these tests was to characterize the SBLI in the presence of a laminar, transitional, and turbulent boundary layer using high-speed optical imaging techniques. AoA = sting angle-of-attack (°) dcylinder = cylinder diameter (mm) dtrip = cylindrical tripping element diameter (mm) Δshock = shock stand-off distance (mm) hcylinder = cylinder height (mm) htrip = cylindrical tripping element height (mm) HSS = high-speed schlieren M∞ = freestream Mach number PSP = pressure-sensitive paint Re∞ = freestream unit Reynolds number (m-1) SWBLI = shockwave/boundary-layer interaction θplate = model plate angle (°) Introduction his paper highlights two experimental entries performed in the 20-Inch Mach 6 Air Blowdown Tunnel at NASA Langley Research Center in collaboration with The University of Tennessee Space Institute (UTSI). The purpose of these entries was to characterize the dynamic shockwave/boundary-layer interaction (SWBLI) between a vertical cylinder on a flat plate and laminar, transitional (XSWBLI), and turbulent (SWTBLI) boundary layers with a freestream Mach number of 6 using non-intrusive optical diagnostics. Experiments performed by Murphree et al.1,2 were among the first to specifically characterize XSWBLI induced by a vertical cylinder on a flat plate geometry using several optical measurement techniques. Recent optical studies of XSWBLI phenomenon have been performed by UTSI at Mach 2 in their low-enthalpy blow wind tunnel3-8 and by Texas A&M University and UTSI at Mach numbers of 6 and 7 in their Adjustable Contour Expansion wind tunnel.9 The experiments described in this paper were intended to complement previous studies by expanding the freestream unit Reynolds number range, Re∞, over which the XSWBLI phenomena has been observed. Additionally these experiments, made possible under NASA’s new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project, promoted collaboration between university and NASA researchers. The initial entry in the 20-Inch Mach 6 Air Tunnel at NASA Langley occurred in December of 2016. Originally, testing was to occur in November of 2016 in the 31-Inch Mach 10 Air Tunnel at NASA Langley. This facility was chosen so that the XSWBLI phenomenon could be observed at much higher Mach numbers than had previously been attempted in ground test experiments. The model selected for this experiment, a 10° half-angle wedge with a sharp leading edge (described in detail in section II.B), had previously been used by Danehy et al. [10] for boundary layer transition studies using the nitric oxide planar laser-induced fluorescence (NO PLIF) flow visualization technique. In that work, it was determined that transition could be induced downstream of a single htrip = 1-mm tall, dtrip = 4-mm diameter cylindrical tripping element and that the streamwise location of the transition could be changed for a single Re∞ by changing the model angle-of-attack (AoA) (see Fig. A3 in Ref. [10] for more details). Based on the findings of that work, a decision was made to use the wedge model with the cylindrical tripping element to trip the boundary layer flow ahead of a cylindrical protuberance in order to achieve a XSWBLI. Unfortunately, the 31-Inch Mach 10 facility had been taken offline for repairs in October of 2016 and a decision was made to move the test to the 20-Inch Mach 6 facility. Since the behavior of the boundary layer with the chosen model configuration had not been studied before in that facility and the available test time was limited, the entry was considered to be exploratory and was used to collect spatially-resolved and time-resolved flow and surface visualization data that would be used to inform a second entry. Test techniques included simultaneous high-speed schlieren (HSS) captured at 160 kHz and high-speed pressure sensitive paint captured at 10 kHz as well as oil flow visualization, captured at 750 Hz. The second entry in the 20-Inch Mach 6 facility occurred in June and July of 2017. In this follow-on test, modifications to the wind tunnel model were made based on observations made during the first entry and included removing the cylindrical tripping element, increasing the size of the cylinder used to induce the SWBLI to increase the size of the interaction while simultaneously improving spatial resolution, and using a swept ramp array, similar to that described in Ref. [11], to trip the flow to turbulence. Simultaneous HSS (captured at 140 kHz, 100 kHz, and 40 kHz) and conventional IR thermography (captured at 30 Hz) imaging were performed simultaneously in this follow-on entry. This paper is intended to serve as a summary of the work performed during these two entries, to detail lessons learned from each entry, and to highlight some of the datasets acquired. Details on the experimental setup, model configuration, and techniques used are provided. Papers providing a more rigorous analysis of data acquired during the second entry, including statistical, spectral, and modal decomposition methods, can be found in Refs. [12,13]. An entry examining XSWBLI in the 31-Inch Mach 10 Blowdown Wind Tunnel facility is currently planned for mid-to-late calendar year 2019, pending the success of facility repairs. The work performed and described in this paper and the upcoming entry in the 31-Inch Mach 10 facility at NASA Langley have been made possible by NASA’s new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project. Wind Tunnel Facility All experiments discussed in this paper were performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center. Specific details pertaining to this facility can be found in Refs. [14,15], with only a brief description of the facility provided here. For both entries, the nominal freestream unit Reynolds number was varied between 1.8×106 m-1 (0.5×106 ft-1) and 26.3×106 m-1 (8×106 ft-1). The nominal stagnation pressure was varied between 0.21 MPa and 3.33 MPa and the nominal stagnation temperature was varied between 480 K and 520 K to achieve the desired Re∞ condition. For all runs, the nominal freestream Mach number was 6. The nearly square test section is 520.7-mm (20.5-inches) wide by 508-mm (20-inches) high. Two 431.8-mm (17-inch) diameter windows made of Corning 7940, Grade 5F schlieren-quality glass serve as the side walls of the tunnel and provide optical access for the high-speed schlieren measurements. A rectangular window made of the same material as the side windows served as the top wall of the test section and provided optical access for the high-speed PSP and oil flow measurements. For the second entry, this top window was replaced with a Zinc Selenide (ZnSe) window with an anti-reflection coating capable of passing IR wavelengths between 8μm and 12μm with greater than 98% transmittance. The model was sting supported by a strut attached to a hydraulic system that allows for the model pitch angle to be adjusted between -5° to +55°. For the first entry, an initial pitch/pause sweep of the model AoA was performed to observe the resulting SWBLI. Ultimately, however, the sting pitch angle for this entry was fixed at +10.0° so that the angle of the top surface of the wedge relative to the streamwise axis of the tunnel (referred to herein as the plate angle, θplate), was θplate = 0°. For the second entry, θplate = 0° and θplate = -13.25° were initially tested with the swept ramp array (discussed in the following section) to determine which orientation produced conditions most favorable for XSWBLI to occur based on the heating signatures observed over the top surface of the model in the IR thermography images. Based on these initial tests, θplate = -13.25° was set for the remainder of the runs in the second entry. For both entries, any model changes were performed in a housing located beneath the closed test section. Prior to performing a run of the tunnel, the housing was sealed and the tunnel started. Once the appropriate freestream conditions were achieved, the model was injected into the test section using a hydraulic injection system. B. Model Geometry For all runs, a 10° half-angle (20° full-angle) wedge model with a sharp leading edge was used. The model is described in detail in Refs. [10,16]. The top surface of the sharp leading edge of the model extended 47.8 mm from its upstream-most edge to a junction with the upstream edge of a stainless steel top plate that then extended an (a) (c) (b) Fig. 1 (a) Schematic of top surface of wedge model with gas seeding insert, (b) perspective view of the model in the 20-Inch Mach 6 tunnel with centerline pressure orifices on sharp leading edge, and (c) a perspective view of the model with stainless steel (top) and SLA middle insert (bottom) during the first entry. Flow occurs from left to right.

Bathel, Brett F.↗

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗

New Mastcam Multispectral Rock Classes in Sulfate-Bearing Strata, Gale Crater, Mars

Orbital data of Gale Crater, Mars has identified a transition in the stratigraphy of Mt. Sharp indicating an environmental change from a wetter one that accommodated clay mineral formation to a drier environment that led to the precipitation of sulfates. Mt. Sharp is not the only location on Mars for which this transition has been observed which suggests a global environmental change around the Noachian-Hesperian transition. Also visible from orbit is a dark-toned marker band that lies within the sulfate-bearing strata hypothesized to be a volcanic ash deposit. The Mars Science Laboratory (MSL) Curiosity rover is currently exploring this region of Mt. Sharp. In-situ analyses by the rover show changes in morphology and chemistry between the sulfate-bearing layers below and above the marker band and the marker band itself. The MSL instrument payload includes the Mast Camera (Mastcam), which collects multispectral images to provide context for other instruments and broad mineralogic interpretations. In addition to mafic Fe-bearing minerals, Mastcam is sensitive to oxides and some hydrated minerals, important indicators of alteration. In our previous work, we compiled a database of Mastcam spectra through Curiosity’s exploration of Vera Rubin ridge (sols 0-2302), from which 9 spectral classes of rocks were identified (Fig. 1D). Each class has diagnostic spectral features that reflect a common mineralogic interpretation. More recently, as Curiosity traversed through the clay-rich region, Glen Torridon (GT), there was some variation in the rock spectral classes potentially due to nontronite and other phyllosilicates, but not enough to warrant new spectral classes (Fig. 1B). We hypothesize that new classes will appear in accordance with sulfate-bearing strata and the marker band, specifically appearing as a ‘downturn’ in the last two wavelength filters that is inherent, but not unique, to polyhydrated Mg-sulfates. Here, we expand upon our analysis of Mastcam multispectral observations through Curiosity’s second encounter with the marker band (up to sol 3672), comparing the spectral diversity within the sulfate-bearing unit and the marker band to the spectral classes encountered previously in the traverse. These analyses will provide a basis for comparison and aid in selecting targets for in-situ investigations as Curiosity continues its ascent of Mt. Sharp.

A M Eng↗

Separation and Transition on a CCF: Experimental Campaigns

Several experimental campaigns have been conducted across a number of Mach-6 facilities on the hypersonic flow around a CCF (CCF) geometry with a 5° half-angle cone and a 12° half-angle flare. These experiments were conducted as part of the NATO STO Research Task Group AVT-346, which is focused on predicting hypersonic boundary-layer transition on complex geometries. Two conventional wind tunnels (AFRL M6LT and ONERA R2Ch) and one quiet tunnel (Purdue BAM6QT) were used to test the same CCF geometry and to study the instabilities in both the boundary layer in the attached parts of the flow and the shear layer above the axisymmetric separation bubble near the cylinder-cone junction. Two nosetip radii (one nominally sharp and one blunt with a 5 mm radius) were tested. For the sharp nosetip case, there was a great deal of agreement between the measurements of both second-mode and shear-layer instabilities across the two conventional facilities. However, the measured spectra for the blunt nosetip case showed more significant differences between the two tunnels, potentially due to an alternate dominant instability mechanism coupled with the variations in the freestream noise spectra. The quiet facility resulted in a flow that remained laminar to much higher freestream unit Reynolds numbers, as well as in instability measurements that had more distinct spectral peaks for the sharp tip case and broadband rises for the blunt one. The instability mechanisms at play in the sharp quiet case were found to be the same as those in the conventional facilities.

Boundary Layer Transition↗

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗

Angular-spectral filtering of recoil protons for optimization of fast neutron imaging employing proton converters

Fast neutron imaging is an important capability for diverse applications such as inertial confinement fusion diagnostics, cargo security, nuclear nonproliferation and arms control, and industrial inspection. Traditional phosphor image plates can be enhanced for fast neutron imaging using hydrogenous plastic converters which allow fast neutrons to scatter off hydrogen nuclei to produce energetic protons that can be recorded by the image plate. However, protons emitted by image plates are not constrained in their emission angle, which contributes to the blur of the resulting image. Here, we investigate two methods that can alter the spatial extent of converted protons that deposit energy in the image plate: reducing the converter thickness, and introducing a proton filter between the plastic converter and image plate to reduce the contribution of lower-energy, off-axis protons to the image. Here we determine the optimal plastic converter thickness for maximizing the signal intensity to be 2–3 mm through Monte Carlo simulations, and we benchmark this result against experimental measurements with a deuterium-tritium (DT) neutron generator. Next, we evaluate the image smearing and signal loss for various converters to show that solely reducing the converter thickness has the expected effect of reducing the blur from proton image smearing of the sharpness of an edge recorded on the image plate at the cost of reducing the signal intensity. The use of a proton filter is shown to achieve a similar improvement of edge sharpness as reducing the converter thickness while also sacrificing the signal intensity. We conclude that the use of proton energy filtering can improve the sharpness of fast neutron images in situations where the converter thickness cannot be reduced below some practical minimum. For more intense neutron sources, the signal intensity is of less concern, and optimizing the resolution of the image plate and therefore of the imaging system could have greater value. In these applications, proton filters may allow for improved fast neutron imaging measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗