SEARCH · Search NASA

Results for “Contour method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

639 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Methodology and Results of the Near-Earth Object (NEO) Human Space Flight (HSF) Accessible Targets Study (NHATS)

Near-Earth Asteroids (NEAs) have been identified by the current administration as potential destinations for human explorers during the mid-2020s. While the close proximity of these objects' orbits to Earth's orbit creates a risk of highly damaging or catastrophic impacts, it also makes some of these objects particularly accessible to spacecraft departing Earth, and this presents unique opportunities for solar system science and humanity's first ventures beyond cislunar space. Planning such ambitious missions first requires the selection of potentially accessible targets from the growing population of nearly 7,800 NEAs. To accomplish this, NASA is conducting the Near-Earth Object (NEO) Human Space Flight (HSF) Accessible Targets Study (NHATS). Phase I of the NHATS was executed during September of 2010, and Phase II was completed by early March of 2011. The study is ongoing because previously undetected NEAs are being discovered constantly, which has motivated an effort to automate the analysis algorithms in order to provide continuous monitoring of NEA accessibility. The NHATS analysis process consists of a trajectory filter and a minimum maximum estimated size criterion. The trajectory filter employs the method of embedded trajectory grids to compute all possible ballistic round-trip mission trajectories to every NEA in the Jet Propulsion Laboratory (JPL) Small-Body Database (SBDB) and stores all solutions that satisfy the trajectory filter criteria. An NEA must offer at least one qualifying trajectory solution to pass the trajectory filter. The Phase II NHATS filter criteria were purposely chosen to be highly inclusive, requiring Earth departure date between January 1st, 2015 and December 31st, 2040, total round-trip flight time <= 450 days, stay time at the NEA >= 8 days, Earth departure C(sub 3) energy <= 60 km(exp 2)/s(exp 2), total mission delta-v <= 12 km/s (including an Earth departure maneuver from a 400 km altitude circular parking orbit), and a maximum atmospheric re-entry speed of 12 km/s. After determining which NEAs offer at least one trajectory solution meeting the criteria, the estimated size constraint is then imposed whereby those NEAs may only be considered NHATS- qualifying NEAs if their maximum estimated size is >= 30 m. This corresponds to an absolute magnitude H <= 26.5 with an assumed albedo p = 0:05. The following is a brief high-level summary of the Phase II study results. Of the 7,665 NEAs in the SBDB as of February 3rd, 2011, 765 NEAs passed the trajectory filter and yielded a total of 79,157,604 trajectory solutions. The trajectory solutions for each NEA are post-processed into Pork Chop Contour (PCC) plots which show total mission delta-v as a function of Earth departure date and total mission duration. Although the PCC plots necessarily compress a very multi-dimensional design space into a two-dimensional plot, they permit rapid assessment of the breadth and quality of an NEA's available Earth departure season and clearly indicate the regions of the trajectory design space which warrant further analysis and optimization. The PCC plot for the NEA with the greatest number of NHATS-qualifying trajectory solutions, 2000 SG-344, is shown. Of the 765 NEAs which passed the Phase II trajectory filter, a total of 590 NEAs also satisfied the further constraint of maximum estimated size >= 30 m. The distributions of osculating heliocentric orbital semi-major axis (a), eccentricity (e), and inclination (i), for those 590 NEAs are shown. Note that the semi-latus rectum used is equal to alpha (1-e(exp 2)). To further our understanding of round-trip trajectory accessibility dynamics, it is instructive to examine the distribution of the NHATS-Qualifying NEAs according to orbit classification. NEAs are grouped into four orbit families: Atiras (aphelion < 0.983 AU), Atens (aphelion > 0.983 AU, alpha < 1.0 AU), Apollos (perihelion < 1.017 AU, alpha > 1.0 AU), and Amors (1.017 < perihelion < 1.3 AU). Of the 765 NEAhich satisfied the NHATS trajectory criteria, none are Atiras, 193 are Atens (31% of known Atens), 456 are Apollos (11% of known Apollos), and 116 are Amors (4% of known Amors). While Apollos comprise 60% of the NEAs which pass the NHATS trajectory filter and Atens comprise only 25%, the percentages according to orbit family are perhaps more relevant. Note that only 11% of known Apollos passed the trajectory filter while 31% of known Atens passed. These simple statistics alone strongly suggest that Aten orbits possess features which tend to enhance their round-trip trajectory accessibility as compared to Apollos or Amors. This is significant because Atens' orbits cause them to spend considerable time in Earth's daytime sky, making them difficult to discover and track using ground-based observing assets. In this paper we will detail the NHATS analysis algorithms, present and analyze all NHATS results to date, and discuss aspects of HSF mission architecture design for future NEA missions.

Barbee, Brent↗

NASA NTRS2006

Impact of Probabilistic Weather on Flight Routing Decisions

Flight delays in the United States have been found to increase year after year, along with the increase in air traffic. During the four-month period from May through August of 2005, weather related delays accounted for roughly 70% of all reported delays, The current weather prediction in tactical (within 2 hours) timeframe is at manageable levels, however, the state of forecasting weather for strategic (2-6 hours) timeframe is still not dependable for long-term planning. In the absence of reliable severe weather forecasts, the decision-making for flights longer than two hours is challenging. This paper deals with an approach of using probabilistic weather prediction for Traffic Flow Management use, and a general method using this prediction for estimating expected values of flight length and delays in the National Airspace System (NAS). The current state-of-the-art convective weather forecasting is employed to aid the decision makers in arriving at decisions for traffic flow and flight planing. The six-agency effort working on the Next Generation Air Transportation System (NGATS) have considered weather-assimilated decision-making as one of the principal foci out of a list of eight. The weather Integrated Product Team has considered integrated weather information and improved aviation weather forecasts as two of the main efforts (Ref. 1, 2). Recently, research has focused on the concept of operations for strategic traffic flow management (Ref. 3) and how weather data can be integrated for improved decision-making for efficient traffic management initiatives (Ref. 4, 5). An overview of the weather data needs and benefits of various participants in the air traffic system along with available products can be found in Ref. 6. Previous work related to use of weather data in identifying and categorizing pilot intrusions into severe weather regions (Ref. 7, 8) has demonstrated a need for better forecasting in the strategic planning timeframes and moving towards a probabilistic description of weather (Ref. 9). This paper focuses on. specified probability in a local region for flight intrusion/deviation decision-making. The process uses a probabilistic weather description, implements that in a air traffic assessment system to study trajectories of aircraft crossing a cut-off probability contour. This value would be useful for meteorologists in creating optimum distribution profiles for severe weather, Once available, the expected values of flight path and aggregate delays are calculated for efficient operations. The current research, however, does not deal with the issue of multiple cell encounters, as well as echo tops, and will be a topic of future work.

Sheth, Kapil↗

DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

CFD Predictions of Boiling Regime Transitions during Line Chilldown validated against a 1G LN2 Experiment

Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.

Evaporation Condensation↗

NASA NTRS2016

Effects of Polar Stratospheric Clouds in the Nimbus 7 LIMS Version 6 Data Set

The historic Limb Infrared Monitor of the Stratosphere (LIMS) measurements of 1978-1979 from the Nimbus 7 satellite were re-processed with Version 6 (V6) algorithms and archived in 2002. The V6 data set employs updated radiance registration methods, improved spectroscopic line parameters, and a common vertical resolution for all retrieved parameters. Retrieved profiles are spaced about every 1.6 of latitude along orbits and include the additional parameter of geopotential height. Profiles of O3 are sensitive to perturbations from emissions of polar stratospheric clouds (PSCs). This work presents results of implementing a first-order screening for effects of PSCs using simple algorithms based on vertical gradients of the O3 mixing ratio. Their occurrences are compared with the co-located, retrieved temperatures and related to the temperature thresholds needed for saturation of H2O and/or HNO3 vapor onto PSC particles. Observed daily locations where the major PSC screening criteria are satisfied are validated against PSCs observed with the Stratospheric Aerosol Monitor (SAM) II experiment also on Nimbus 7. Remnants of emissions from PSCs are characterized for O3 and HNO3 following the screening. PSCs may also impart a warm bias in the co-located LIMS temperatures, but by no more than 1-2K at the altitudes of where effects of PSCs are a maximum in the ozone; thus, no PSC screening was applied to the V6 temperatures. Minimum temperatures vary between 187 and 194K and often occur 1 to 2 km above where PSC effects are first identified in the ozone (most often between about 21 and 28 hPa). Those temperature-pressure values are consistent with conditions for the existence of nitric acid trihydrate (NAT) mixtures and to a lesser extent of super-cooled ternary solution (STS) droplets. A local, temporary uptake of HNO3 vapor of order 1-3 ppbv is indicated during mid-January for the 550K surface. Seven-month time series of the distributions of LIMS O3 and HNO3 are shown based on their gridded Level 3 data following the PSC screening. Zonal coefficients of both species are essentially free of effects from PSCs on the 550K surface, based on their average values along PV contours and in terms of equivalent latitude. Remnants of PSCs are still present in O3 on the 450K surface during mid- January. It is judged that the LIMS Level 3 data are of good quality for analyzing the larger-scale, stratospheric chemistry and transport processes during the Arctic winter of 1978- 1979.

Remsberg, Ellis↗

NASA NTRS2015

On Hilbert-Huang Transform Based Synthesis of a Signal Contaminated by Radio Frequency Interference or Fringes

Norden E. Huang et al. had proposed and published the Hilbert-Huang Transform (HHT) concept correspondently in 1996, 1998. The HHT is a novel method for adaptive spectral analysis of non-linear and non-stationary signals. The HHT comprises two components: - the Huang Empirical Mode Decomposition (EMD), resulting in an adaptive data-derived basis of Intrinsic Mode functions (IMFs), and the Hilbert Spectral Analysis (HSA1) based on the Hilbert Transform for 1-dimension (1D) applied to the EMD IMF's outcome. Although paper describes the HHT concept in great depth, it does not contain all needed methodology to implement the HHT computer code. In 2004, Semion Kizhner and Karin Blank implemented the reference digital HHT real-time data processing system for 1D (HHT-DPS Version 1.4). The case for 2-Dimension (2D) (HHT2) proved to be difficult due to the computational complexity of EMD for 2D (EMD2) and absence of a suitable Hilbert Transform for 2D spectral analysis (HSA2). The real-time EMD2 and HSA2 comprise the real-time HHT2. Kizhner completed the real-time EMD2 and the HSA2 reference digital implementations respectively in 2013 & 2014. Still, the HHT2 outcome synthesis remains an active research area. This paper presents the initial concepts and preliminary results of HHT2-based synthesis and its application to processing of signals contaminated by Radio-Frequency Interference (RFI), as well as optical systems' fringe detection and mitigation at design stage. The Soil Moisture Active Passive (SMAP mission (SMAP) carries a radiometer instrument that measures Earth soil moisture at L1 frequency (1.4 GHz polarimetric - H, V, 3rd and 4th Stokes parameters). There is abundant RFI at L1 and because soil moisture is a strategic parameter, it is important to be able to recover the RFI-contaminated measurement samples (15% of telemetry). State-of-the-art only allows RFI detection and removes RFI-contaminated measurements. The HHT-based analysis and synthesis facilitates recovery of measurements contaminated by all kinds of RFI, including jamming [7-8]. The fringes are inherent in optical systems and multi-layer complex contour expensive coatings are employed to remove the unwanted fringes. HHT2-based analysis allows test image decomposition to analyze and detect fringes, and HHT2-based synthesis of useful image.

Spacecraft Electronics Equipment↗

NASA NTRS2019

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.↗

DOE OSTI2023

IceCat-1: The IceCube Event Catalog of Alert Tracks

We present a catalog of likely astrophysical neutrino track-like events from the IceCube Neutrino Observatory. IceCube began reporting likely astrophysical neutrinos in 2016, and this system was updated in 2019. The catalog presented here includes events that were reported in real time since 2019, as well as events identified in archival data samples starting from 2011. We report 275 neutrino events from two selection channels as the first entries in the catalog, the IceCube Event Catalog of Alert Tracks, which will see ongoing extensions with additional alerts. The Gold and Bronze alert channels respectively provide neutrino candidates with a 50% and 30% probability of being astrophysical, on average assuming an astrophysical neutrino power-law energy spectral index of 2.19. For each neutrino alert, we provide the reconstructed energy, direction, false-alarm rate, probability of being astrophysical in origin, and likelihood contours describing the spatial uncertainty in the alert's reconstructed location. We also investigate a directional correlation of these neutrino events with gamma-ray and X-ray catalogs, including 4FGL, 3HWC, TeVCat, and Swift-BAT.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2005

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗