Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plot”

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 559 records · Page 31

Preliminary Characterization of IDCSP Spacecrafts Through a Multi-Analytical Approach

Defining the risks present to both crewed and robotic spacecrafts is part of NASA s mission, and is critical to keep these resources out of harm s way. Characterizing orbital debris is an essential part of this mission. We present a proof-of-concept study that employs multiple techniques to demonstrate the efficacy of each approach. The targets of this study are IDCSPs (Initial Defense Communications Satellite Program). 35 of these satellites were launched by the US in the mid-1960s and were the first US military communications satellites in the GEO regime. They were emplaced in slightly sub-synchronous orbits. These targets were chosen for this proof-of-concept study for the simplicity of their observable exterior surfaces. The satellites are 26-sided polygons (86cm in diameter), initially spin-stabilized, and covered on all sides in solar panels. Data presented here include: (a) visible broadband photometry (Johnson/Kron-Cousins BVRI) taken with the 0.9m SMARTs telescope (Small and Medium Aperture Telescopes) at the Cerro Tololo Inter-American Observatory (CTIO) in Chile in April, 2012, (b) laboratory broadband photometry (Johnson/ Bessell BVRI) of solar cells, obtained using the Optical Measurements Center (OMC) at NASA/JSC [1], (c) visible-band spectra taken using the Magellan 6.5m Baade Telescope at Las Campanas Observatory in Chile in May, 2012 [2], and (d) visible-band laboratory spectra of solar cells using an ASD Field Spectrometer. Color-color plots using broadband photometry (e.g. B-R vs. B-V) demonstrate that different material types fall into distinct areas on the plots [1]. Spectra of the same material types as those plotted in the color-color plots each display their own signature as well. Here, we compare lab data with telescopic data, and photometric results with spectroscopic results. The spectral response of solar cells in the visible wavelength regime varies from relatively flat to somewhat older solar cells whose reflectivity can be gently or sharply peaked in the blue. With a target like IDCSPs, the material type is known a priori, aiding in understanding how material type affects one s observations.

Lederer, S. M.↗

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI.

Comparison↗

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗

A Data Exploration Tool for Large Sets of Spectra

We present an exploration tool for very large spectrum data sets such as the SDSS (Sloan Digital Sky Survey), LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope), and 4MOST (4-meter Multi-Object Spectroscopic Telescope) data sets. The tool works in two stages: the first uses batch processing and the second runs interactively. The latter employs the NASA hyperwall, a configuration of 128 workstation displays (8 by 16 array) controlled by a parallelized software suite running on NASA's Pleiades supercomputer. The stellar subset of the Sloan Digital Sky Survey, DR10, was chosen to show how the our tool may be used. In stage one, SDSS files for 569,740 stars are processed through our data pipeline. The pipeline fits each spectrum using an iterative continuum algorithm, distinguishing emission from absorption and handling molecular absorption bands correctly. It then measures 1659 discrete atomic and molecular spectral features that were carefully preselected based on their likelihood of being visible at some spectral type. The depths relative to the local continuum at each feature wavelength are determined for each spectrum: these depths, the local S/N (signal to noise ratio) level, and DR10-supplied variables such as magnitudes, colors, positions, and radial velocities are the basic measured quantities used on the hyperwall. In stage two, each hyperwall panel is used to display a 2-D scatter plot showing the depth of feature A vs the depth of feature B for all of the stars. A and B change from panel to panel. The relationships between the various (A,B) strengths and any distinctive clustering are immediately apparent when examining and inter-comparing the different panels on the hyperwall. The interactive software allows the user to select the stars in any interesting region of any 2-D plot on the hyperwall, immediately rendering the same stars on all the other 2-D plots in a unique color. The process may be repeated multiple times, each selection displaying a distinctive color on all the plots. At any time, the spectra of the selected stars may be examined in detail on a connected workstation display. We illustrate how our approach allows us to quickly isolate and examine such interesting stellar subsets as EMP (Extremely Metal‐Poor) stars, CV (Cataclymic Variable) stars and C (Carbon)-rich stars.

Data Exploration↗

Bringing Analysis Closer to Data: Developing a Visualization Tool for L2 Earth Science Satellite Data

Earth Science satellite missions provide a unique opportunity for scientists to visualize complex and multifaceted observations projected geospatially across maps of the Earth. While visualization tools can help scientists comprehend, analyze, and share data, visualizing Level-2 Earth Sciences data poses its own specific set of challenges. Since the geospatial information in Level-2 data files is stored as independent variables, the plotting process involves matching dimensional information from latitude and longitude with a desired variable. Variables are stored in different ways across various Earth Science data file formats, which complicates the process of extracting data and plotting variables from a given file without requiring extensive user input and prerequisite familiarity with the file type variable structure. In coordination with NASA’s Goddard Earth Sciences Data Information Services Center (GES DISC), the team developed a Level-2 Earth Science data visualization tool that aims to address some of the complexities associated with plotting Level-2 data. This tool offers command-line and user interface support for file and variable selection to accommodate varying use cases and degrees of user familiarity with the structure of a given file. The visualization tool is written in Python 3 and utilizes a modular approach to facilitate continued expansion and reuse. In addressing some common complications involved in plotting Level-2 Earth Sciences data, the tool aims to help to link the process of analysis more directly with data acquisition and visualization, bringing analysis closer to data across levels of processing.

Li, Angela W.↗

Assessing Alaskan boreal forest landcover affected by climate-wildfire interactions from ground truth surveys and NASA airborne remote sensing

Alaska’s boreal forest is facing unprecedented challenges under rapid climate warming (increasingly severe fires, droughts, pest/disease outbreaks) that may destabilize its function as a global carbon sink. Forests near Fairbanks may be especially vulnerable, impacting air quality and ecosystem services. We combined GT (ground truthing) with Airborne Visible InfraRed Imaging Spectrometer (AVIRIS-NG) images collected by the NASA Arctic-Boreal Vulnerability Experiment (ABoVE) program (2017-2019) to assess landcover change at five recently burned sites (2001-2019) of different fire severities and moisture regimes within 30 miles of Fairbanks. GT included tree seedling counts, understory % cover and >50% leaf canopy color assessment. 36 circular plots (1/30 ha radius) including 6 moderate to severely burned plots were selected across sites. 31 additional sites including 12 burned sites were geotagged in photos. AVIRIS images were processed from 29 spectral bands selected to identify changes in chlorophyll and water content. Images were segmented into natural boundaries (polygons) using ENVI 5.5 software. A spectral library of 8 AVIRIS bands with high between-class/low within-class variation was used in two random forest models to predict vegetation classes (model 1: 12 classes, model 2: 14 classes) in each AVIRIS scene, using 20% of the data as training data. Model 2 classified 20% more polygons overall, but only 42% of GT/geotagged polygons were correctly classified by both models. More forest sites were correctly classified (63%) than open vegetation (32%) or post-fire sites (46%). 50% of aspen forest and post-fire polygons were misclassified as shrubland. GT revealed that post-fire plots supported 134,000 (± 48,000) tree seedlings and saplings ha-1 (0.2 - 4 m height, 64% deciduous) versus 2500 (± 2100) shrubs ha-1 (1-6 m height). > 50% canopy browning was observed in conifer forest (8 plots) with no signs of insect infestation. Canopy herbivory > 50% (leaf miner, leaf beetle) and moose herbivory of tree bark was seen across aspen sites. Our study suggests: 1) low canopy vegetation presents challenges for improved landcover classification, and 2) aspen forest should be differentiated in vegetation maps which would aid in tracking herbivory.

Alaska↗

NASA GeneLab Multi-study Visualization Portal

NASA GeneLab has helped advance the field of Space Biology by providing a public repository where researchers can store, share, analyze and visualize the results of space flight related omics experiments. The GeneLab data visualization portal allows any user, regardless of bioinformatics knowledge or access to computational resources, to interact with the experimental data, draw their own conclusions, and gain insights about the effects of space on living systems. These tools help democratize scientific research and foster the NASA Open Science initiative. The new multi-study feature of the GeneLab visualization platform allows users to mine study metadata from RNA sequencing (RNA-seq) experiments to identify samples of interest by filtering datasets based on organism, tissue, assay technology type, and/or factor. Once samples are selected from multiple datasets, users can combine and normalize the sample data, then utilize the visualization displays, including Principal Component Analysis (PCA) plots, to assess sample distributions. Finally, users can perform differential gene expression analysis on the combined data and visualize the results through PCA plots, Volcano plots, Pair plots, Heatmap, Ideogram and Gene Set Enrichment Analysis. All user-generated results and visualizations will be available for download. Here, we present a biological study using samples from multiple GeneLab RNA-seq datasets and analyzed using the multi-study visualization platform to demonstrate inter- and intra-study variability, as well as commonly differentially expressed genes between spaceflight and ground control conditions across datasets. This new feature opens a wide range of possibilities and opportunities for further development including combining other assay technology types and integration with batch effect correction techniques and machine learning applications. Overall, this tool allows users to increase the statistical power of individual experiments, validate hypothesis, identify patterns, and opens the door to new and exciting research.

space biology↗

Open-source Numerical Modeling of Solidification Cracking Susceptibility: Application to Refractory Alloy Systems

Introduction. Alloys such as aluminum, nickel-base, and austenitic stainless steels are susceptible to solidification cracking during welding and 3D printing. Compositional optimization is one method used to effectively mitigate solidification cracking of those alloy systems. With the surge in hypersonic and in-space propulsion activities, refractory metals (Nb, Mo, Ta, W, and Re) and their alloy derivatives are increasing in importance due to their extreme high melting point and retention of high-temperature strength; however, their chemistry was most typically optimized to promote ductility during mechanical operations such as drawing and forming. Welding of such alloys has been a challenge due to a number of issues including solidification cracking, atmospheric contamination (O, C, and N), as well as a shift in ductile-to-brittle transition to higher temperature following grain growth induced by welding. Compositional optimization of refractory alloys for solidification cracking resistance in particular is desirable as their usage increases with the advent of advanced manufacturing methods such as 3D printing. This work evaluates the effect of compositional variation in refractory metal systems on the solidification cracking susceptibility with the goals of optimizing existing alloys and joining process techniques, and formulating new alloys with increased solidification cracking resistance. Experimental Procedures. A python code was developed in a Jupyter notebook environment (Michael and Sowards, 2023) to facilitate the calculation of crack susceptibility index proposed by Kou (2015). Composition is entered as a single point, or as a 1-D or 2-D array. The notebook calls pycalphad (Otis and Liu, 2017 and Bocklund et al, 2020) to calculate the evolution of fraction solid as a function of temperature (under either Scheil or equilibrium assumptions) and then evaluates steepness of the fraction solid curve near the terminal stage of solidification to predict solidification cracking resistance. Open source thermodynamic databases available at online repositories are used (van de Walle). The process is setup in an automated fashion to generate plots that show variation in solidification cracking susceptibility according to composition on 1-D line plots or 2-D contour plots. The Jupyter notebook and crack susceptibility algorithm was also integrated with a widely used commercial CALPHAD code for validation and alloy exploration. Results and Discussion. The crack susceptibility model was first validated against a series of refractory alloy compositions evaluated in past work which utilized a specialized Varestraint test built inside a vacuum chamber environment (Lessman and Gold, 1971). The alloys tested in the Varestraint apparatus included T-111 (Ta-8W-2Hf), ASTAR-811C (Ta-8W-1Re-0.7Hf-0.025C), FS-85 (Nb-27Ta-10W-1Zr), T-222 (Ta-9.6W-2.4Hf-0.01C), Ta-10W, B-66 (Nb-5Mo-5V-1Zr), and SCb-291 (Nb-10W-10Ta). The initial test of the model showed a strong correlation with empirical Varestraint data, i.e., a Spearman rank correlation between model predictions and hot cracking measurements was observed to be greater than 0.8. Following the validation, a set of refractory metal binary mixtures was investigated to evaluate sensitivity of Nb, Mo, W, and Ta to C, N, and O content. A series of plots were produced that suggest ppmw ranges of C, N, and O where solidification cracking increases significantly and reaches a maximum. Also comparative ranking of each primary refractory metal to each interstitial was produced. For example C produces greater cracking response in Mo whereas O produces greater cracking response in Ta and Nb. Such compositional values have utility in setting limits on pickup of these interstitial elements during welding and printing rather than using a one-size-fits-all approach. Furthermore, the results have use in determining additive powder recycling requirements, which is especially pertinent for refractory metal powders due to their high cost compared to conventional alloys. Another application created thousands of hypothetical alloys within the nominal specified composition range of two widely used refractory alloys C103 (Nb-10Hf-1Ti) and TZM (Mo-0.5Ti-0.1Zr). The cracking index was calculated for the alloys and results were fed into machine learning regression techniques including Multiple Linear Regression, Ridge Regression, and Lasso Regression to determine relative potency each alloying element had on computed solidification cracking index. A series of linear equations were produced that relate composition of C103 and TZM to solidification cracking index. The crack susceptibility of C103 for example is described by an equation of the form: cracking index ~ O + 0.667*C + 0.635*N + 0.00037*Ta – 0.0008*Hf (in wt.%) From that equation, it is clear that O has strong propensity to induce solidification cracking. Interestingly, Hf is shown to reduce calculated cracking response. Finally, realizing the potential of this method to discover new refractory alloy formulations across the period table that have low solidification cracking sensitivity, the code was applied to new untested alloy systems including W-Zr-C, W-Ta-C, and others. Conclusions. In summary, an open source numerical method has been developed using Python code to calculate Kou’s crack susceptibility index. The method was applied to refractory metals which are inherently difficult to study from a weldability testing standpoint since inert shielding gas is not sufficient and welding is typically done in vacuum, especially in light of findings presented here where oxygen has profound influence on solidification cracking. This work revealed the effect of compositional variations on a series of refractory metals and showed the framework defined here will be useful in 1) the development of new alloys that have improved weldability and 3D printability, 2) placing compositional limits on existing alloys, and 3) ensuring adequate controls of manufacturing processes such as 3D printing where powder reuse is critical. Keywords. pycalphad; Python; refractory metals; solidification cracking. References. B. Bocklund et. al. (2020) http://doi.org/10.5281/zenodo.3630657. S. Kou. (2015) https://doi.org/10.1016/j.actamat.2015.01.034. G.G. Lessmann and R.E. Gold. Welding Journal, issue 1, pp. 1-s – 8-s (1971). F.N. Michael and J.W. Sowards. NASA/TM-20230002218 (2023). R. Otis and Z.-K. Liu. (2017) http://doi.org/10.5334/jors.140. A. Van de Wallle et. al. (2018) https://doi.org/10.1016/j.calphad.2018.04.003.

pycalphad↗

The impact of tillage row orientation on physical and chemical sediment enrichment

Abstract This study aimed to better understand how tillage row orientation with respect to dominant flow‐pathway along hillslope impacts runoff and the transport of different sediment size fractions. Experimental plots were constructed in contour ridge till (CRT) and parallel ridge till (PRT) sites to monitor runoff and sediment fluxes. Particle size fractions of the sediment, along with organic C and N contents were measured to quantify physical and chemical enrichment ratios for the two tillage orientations. In the CRT plot, tillage produced large oriented roughness elements along the contours, which acted as little check dams, while in the PRT site, the roughness elements helped confine and concentrate the runoff. In the CRT, the “check dams” resulted in a runoff coefficient of .03, while in the PRT, the flow confinement between rows produced a runoff coefficient of .82. Moreover, the erosion rates at the CRT site were 97% less than those in the PRT plot. Large contours produced finer sediment fractions due to the selective sorting in ponded furrows. More aggregate sediment fractions were present in the PRT site, which was dominated by rill erosion. Physical enrichment revealed selective entrainment of finer sediment particles during erosion. Finer‐grain particles with higher specific surface areas, attached more organic C resulting in chemical enrichment. Physical and chemical enrichment methods were found to be in good agreement. These findings suggest transport models that can simulate size fraction updates to the soil active layer can be used to estimate SOC redistribution and hence more accurate C budgets.

Wacha, Kenneth M.↗

PAH/PAH(CF 3 ) n Donor/Acceptor Charge‐Transfer Complexes in Solution and in Solid‐State Co‐Crystals

Abstract A solution, solid‐state, and computational study is reported of polycyclic aromatic hydrocarbon PAH/PAH(CF 3 ) n donor/acceptor (D/A) charge‐transfer complexes that involve six PAH(CF 3 ) n acceptors with known gas‐phase electron affinities that range from 2.11(2) to 2.805(15) eV and four PAH donors, including seven CT co‐crystal X‐ray structures that exhibit hexagonal arrays of mixed π‐stacks with 1/1, 1/2, or 2/1 D/A stoichiometries (PAH=anthracene, azulene, coronene, perylene, pyrene, triphenylene; n= 5, 6). These are the first D/A CT complexes with PAH(CF 3 ) n acceptors to be studied in detail. The nine D/A combinations were chosen to allow several structural and electronic comparisons to be made, providing new insights about controlling D/A interactions and the structures of CT co‐crystals. The comparisons include, among others, CT complexes of the same PAH(CF 3 ) n acceptor with four PAH donors and CT complexes of the same donor with four PAH(CF 3 ) n acceptors. All nine CT complexes exhibit charge‐transfer bands in solution with λ max between 467 and 600 nm. A plot of E ( λ max ) versus [ IE (donor)− EA (acceptor)] for the nine CT complexes studied is linear with a slope of 0.72±0.03 eV eV −1 . This plot is the first of its kind for CT complexes with structurally related donors and acceptors for which precise experimental gas‐phase IE s and EA s are known. It demonstrates that conclusions based on the common assumption that the slope of a CT E ( λ max ) versus [ IE − EA ] plot is unity may be incorrect in at least some cases and should be reconsidered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tropical understory herbaceous community responds more strongly to hurricane disturbance than to experimental warming

The effects of climate change on tropical forests may have global consequences due to the forests’ high biodiversity and major role in the global carbon cycle. In this study, we document the effects of experimental warming on the abundance and composition of a tropical forest floor herbaceous plant community in the Luquillo Experimental Forest, Puerto Rico. This study was conducted within Tropical Responses to Altered Climate Experiment (TRACE) plots, which use infrared heaters under free-air, open-field conditions, to warm understory vegetation and soils + 4°C above nearby control plots. Hurricanes Irma and María damaged the heating infrastructure in the second year of warming, therefore, the study included one pretreatment year, one year of warming, and one year of hurricane response with no warming. We measured percent leaf cover of individual herbaceous species, fern population dynamics, and species richness and diversity within three warmed and three control plots. Results showed that one year of experimental warming did not significantly affect the cover of individual herbaceous species, fern population dynamics, species richness, or species diversity. In contrast, herbaceous cover increased from 20% to 70%, bare ground decreased from 70% to 6%, and species composition shifted pre to posthurricane. The negligible effects of warming may have been due to the short duration of the warming treatment or an understory that is somewhat resistant to higher temperatures. Our results suggest that climate extremes that are predicted to increase with climate change, such as hurricanes and droughts, may cause more abrupt changes in tropical forest understories than longer-term sustained warming.

54 ENVIRONMENTAL SCIENCES↗

TATB Sensitivity to Shocks from Electrical Arcs

Use of insensitive high explosives (IHEs) has significantly improved ammunition safety because of their remarkable insensitivity to violent cook–off, shock and impact. Triamino–trinitrobenzene (TATB) is the IHE used in many modern munitions. Previously, lightning simulations in different test configurations have shown that the required detonation threshold for standard density TATB at ambient and elevated temperatures (250 °C) has a sufficient margin over the shock caused by an arc from the most severe lightning. In this paper, the Braginskii model with Lee–More channel conductivity prescription is used to demonstrate how electrical arcs from lightning could cause detonation in TATB. The steep rise and slow decay in typical lightning pulse are used in demonstrating that the shock pressure from an electrical arc, after reaching the peak, falls off faster than the inverse of the arc radius. For detonation to occur, two necessary detonation conditions must be met: the Pop–Plot criterion and minimum spot size requirement. The relevant Pop–Plot for TATB at 250 °C was converted into an empirical detonation criterion, which is applicable to explosives subject to shocks of variable pressure. The arc cross–section was required to meet the minimum detonation spot size reported in the literature. Furthermore, one caveat is that when the shock pressure exceeds the detonation pressure the Pop–Plot may not be applicable, and the minimum spot size requirement may be smaller.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Theoretical Study of Superconducting Gap Parameters, Density of States, and Condensation Energy of Two-Band Iron-Based Superconductor BaFe{sub 2}(As{sub 1−x}P{sub x}){sub 2}

This research work focuses on the theoretical study of superconducting gap parameters, density of states, and condensation energy of two-band iron-based superconductor BaFe{sub 2}(As{sub 1−x}P{sub x}){sub 2}. By developing a model Hamiltonian for the given system and by using the double time temperature-dependent Green’s function formalism, we obtained mathematical expressions for superconducting order parameters for the electron intra-band (Δ{sub e}), hole intra-band (Δ{sub h}), inter-band between the two bands (Δ{sub eh}) and superconducting transition temperature (T{sub C}). Furthermore, we obtained the density of states for each intra-band, N(ε) and condensation energy (E{sub c}) for BaFe{sub 2}(As{sub 1−x}P{sub x}){sub 2}. By using the experimental values and by considering some plausible approximations of the parameters in the obtained expressions, phase diagrams of superconducting order parameters for Δ{sub e}, Δ{sub h}, and Δ{sub eh} versus temperature and superconducting transition temperature (T{sub C}) versus the inter-band pairing potential (V{sub eh}) are plotted. We also plotted the density of states for the electron intra-band, N{sub e}(ε) and hole intra-band, N{sub h}(ε) versus excitation energy (ε) at T = 0K, and the variation of density of states with temperature. Similarly, phase diagrams of condensation energy (E{sub c}) versus temperature, inter-band pairing potential, and superconducting transition temperature are plotted for the material. Our results are in agreement with previous findings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nanoindentation Stress–Strain for Fracture Analysis and Computational Modeling for Hardness and Modulus

Stress–Strain plots based on nanoindentation load–depth curves were obtained to study phenomena like internal fracture and ductile to brittle transitions. Fracture phenomena during the indentation process were analyzed based on the stress–strain plots. A transition from ductile to brittle fracture was observed on increasing the depth or load of indenter penetration. A new approach with shape factors in the fracture studies based on radial crack branching and micro-cracking was done. Hardness and modulus plots were fitted with polynomials. The fitting parameters were varied to obtain different hardness and modulus responses.

36 MATERIALS SCIENCE↗

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗

Linking soil biodiversity and ecosystem function in a Neotropical savanna

Conserving the remaining savanna ecosystems in the Brazilian savanna (Cerrado) — a global biodiversity hotspot that stores carbon and provides water to a large portion of South America — requires understanding the ecological processes maintaining their function. Nutrient cycling supports savanna function via plant litter production and decomposition by soil fauna, releasing nutrients for plant and soil organism uptake. Soil biodiversity and biogeochemistry linkages with litter dynamics in Neotropical savannas under a changing climate are poorly understood. Here, in this study, we combined two years of rainfall seasonality, leaf and wood litter production and decomposition with soil epigeic fauna abundance — the number of ground-surface dwelling invertebrates collected through pitfall traps — taxa richness, Shannon's diversity and Pielou's evenness, and 16 soil biogeochemical variables measured in 12 plots of preserved savanna. Rainfall seasonality modulated the mean soil epigeic fauna diversity and evenness across all plots, which were highest in the rainy season, in contrast to litterfall rates, which peaked in the dry season. In the dry season (April to September), the Formicidae family was the most abundant with 50% of all individuals, while in the rainy season (October to March), the Isoptera order was the most abundant with approximately 39% of individuals. Wood litter decomposition grouped with annual Hemiptera abundance, co-varying with soil epigeic fauna diversity and evenness per plot and against soil fertility variables. Leaf litter decomposition co-varied with the total epigeic fauna abundance and soil pH. We speculate that the specific need to decompose wood litter may be associated with a greater need for diversity than an abundance of soil epigeic fauna. Our work highlights the role of rainfall seasonality on soil biodiversity and physicochemistry, which is also tightly linked with litter production and decomposition. This study advances our understanding of the mechanisms governing nutrient cycling in savanna ecosystems on nutrient-impoverished soils, with implications for achieving sustainable conservation and restoration goals.

59 BASIC BIOLOGICAL SCIENCES↗

PyProcar: A Python library for electronic structure pre/post-processing

The PyProcar Python package plots the band structure and the Fermi surface as a function of site and/or s,p,d,f - projected wavefunctions obtained for each k-point in the Brillouin zone and band in an electronic structure calculation. This can be performed on top of any electronic structure code, as long as the band and projection information is written in the PROCAR format, as done by the VASP and ABINIT codes. PyProcar can be easily modified to read other formats as well. This package is particularly suitable for understanding atomic effects into the band structure, Fermi surface, spin texture, etc. PyProcar can be conveniently used in a command line mode, where each one of the parameters define a plot property. In the case of Fermi surfaces, the package is able to plot the surface with colors depending on other properties such as the electron velocity or spin projection. Furthermore, the mesh used to calculate the property does not need to be the same as the one used to obtain the Fermi surface. A file with a specific property evaluated for each k-point in a k-mesh and for each band can be used to project other properties such as electron–phonon mean path, Fermi velocity, electron effective mass, etc. Another existing feature refers to the band unfolding of supercell calculations into predefined unit cells.

97 MATHEMATICS AND COMPUTING↗