Search NASA⌕ Search

SEARCH · Search NASA

Results for “Connection length”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,009 records · Page 56

Tuning interfacial ferromagnetism in LaNiO 3 / CaMnO 3 superlattices by stabilizing nonequilibrium crystal symmetry

Perovskite oxide heterostructures offer an important path forward for stabilizing and controlling low-dimensional magnetism. One of the guiding design principles for these materials systems is octahedral connectivity. In superlattices composed of perovskites with different crystal symmetries, variation of the relative ratio of the constituent layers and the individual layer thicknesses gives rise to nonequilibrium crystal symmetries that, in turn, lead to unprecedented control of interfacial ferromagnetism. We have found that in superlattices of CaMnO 3 (CMO) and LaNiO 3 (LNO), interfacial ferromagnetism can be modulated by a factor of 3 depending on LNO and CMO layer thicknesses as well as their relative ratio. Such an effect is only possible due to the nonequilibrium crystal symmetries at the interfaces and can be understood in terms of the anisotropy of the exchange interactions and modifications in the interfacial Ni-O-Mn and Mn-O-Mn bond angles and lengths with increasing LNO layer thickness. Here, these results demonstrate the potential of engineering nonequilibrium crystal symmetries in designing ferromagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tradeoffs and Synergies in Tropical Forest Root Traits and Dynamics for Nutrient and Water Acquisition: Field and Modeling Advances

Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Operational Lessons Learned from the Ares I-X Flight Test

The Ares I-X flight test, launched in 2009, is the first test of the Ares I crew launch vehicle. This development flight test evaluated the flight dynamics, roll control, and separation events, but also provided early insights into logistical, stacking, launch, and recovery operations for Ares I. Operational lessons will be especially important for NASA as the agency makes the transition from the Space Shuttle to the Constellation Program, which is designed to be less labor-intensive. The mission team itself comprised only 700 individuals over the life of the project compared to the thousands involved in Shuttle and Apollo missions; while missions to and beyond low-Earth orbit obviously will require additional personnel, this lean approach will serve as a model for future Constellation missions. To prepare for Ares I-X, vehicle stacking and launch infrastructure had to be modified at Kennedy Space Center's Vehicle Assembly Building (VAB) as well as Launch Complex (LC) 39B. In the VAB, several platforms and other structures designed for the Shuttle s configuration had to be removed to accommodate the in-line, much taller Ares I-X. Vehicle preparation activities resulted in delays, but also in lessons learned for ground operations personnel, including hardware deliveries, cable routing, transferred work and custodial paperwork. Ares I-X also proved to be a resource challenge, as individuals and ground service equipment (GSE) supporting the mission also were required for Shuttle or Atlas V operations at LC 40/41 at Cape Canaveral Air Force Station. At LC 39B, several Shuttle-specific access arms were removed and others were added to accommodate the in-line Ares vehicle. Ground command, control, and communication (GC3) hardware was incorporated into the Mobile Launcher Platform (MLP). The lightning protection system at LC 39B was replaced by a trio of 600-foot-tall towers connected by a catenary wire to account for the much greater height of the vehicle. Like Shuttle, Ares I-X will be stacked on a MLP and rolled out to the pad on a Saturn-era crawler-transporter. While Ares I-X was only held in place by the four hold-down posts on its aft skirt during rollout, a new vehicle stabilization system (VSS) attached to the vertical service structure kept the vehicle from undue swaying prior to launch at the pad, LC 39B. Following the launch, the flight test vehicle first stage was recovered with the aid of new parachutes resized to accommodate the five-segment-long first stage, which had a much greater length and mass than the Shuttle s reusable solid rocket boosters. After splashdown, recovery divers exercised extra care when handling the first stage to ensure that the flight data recorders in the fifth segment simulator were not damaged by exposure to sea water. The data recovered from the Ares I-X flight test will be very valuable in verifying the predicted environments and models used to design the vehicle. Lessons learned from Ares I-X will be shared with the Ares Projects through written and verbal reports and through integration of mission team members into the Project workforce.

Davis, Stephan R.↗

BioWires: Conductive DNA Nanowires in a Computationally-Optimized, Synthetic Biological Platform for Nanoelectronic Fabrication

DNA is an ideal template for a biological nanowire-it has a linear structure several atoms thick; it possesses addressable nucleobase geometry that can be precisely defined; and it is massively scalable into branched networks. Until now, the drawback of DNA as a conducting nanowire been, simply put, its low conductance. To address this deficiency, we extensively characterize a chemical variant of canonical DNA that exploits the affinity of natural cytosine bases for silver ions. We successfully construct chains of single silver ions inside double-stranded DNA, confirm the basic dC-Ag+-dC bond geometry and kinetics, and show length-tunability dependent on mismatch distribution, ion availability and enzyme activity. An analysis of the absorbance spectra of natural DNA and silver-binding, poly-cytosine DNA demonstrates the heightened thermostability of the ion chain and its resistance to aqueous stresses such as precipitation, dialysis and forced reduction. These chemically critical traits lend themselves to an increase in electrical conductivity of over an order of magnitude for 11-base silver-paired duplexes over natural strands when assayed by STM break junction. We further construct and implement a genetic pathway in the E. coli bacterium for the biosynthesis of highly ionizable DNA sequences. Toward future circuits, we construct a model of transcription network architectures to determine the most efficient and robust connectivity for cell-based fabrication, and we perform sequence optimization with a genetic algorithm to identify oligonucleotides robust to changes in the base-pairing energy landscape. We propose that this system will serve as a synthetic biological fabrication platform for more complex DNA nanotechnology and nanoelectronics with applications to deep space and low resource environments.

DNA↗

Two alternative binding mechanisms connect the protein translocation Sec71-Sec72 complex with heat shock proteins

The biosynthesis of many eukaryotic proteins requires accurate targeting to and translocation across the endoplasmic reticulum membrane. Post-translational protein translocation in yeast requires both the Sec61 translocation channel, and a complex of four additional proteins: Sec63, Sec62, Sec71, and Sec72. The structure and function of these proteins are largely unknown. This pathway also requires the cytosolic Hsp70 protein Ssa1, but whether Ssa1 associates with the translocation machinery to target protein substrates to the membrane is unclear. Here, we use a combined structural and biochemical approach to explore the role of Sec71-Sec72 subcomplex in post-translational protein translocation. To this end, we report a crystal structure of the Sec71-Sec72 complex, which revealed that Sec72 contains a tetratricopeptide repeat (TPR) domain that is anchored to the endoplasmic reticulum membrane by Sec71. We also determined the crystal structure of this TPR domain with a C-terminal peptide derived from Ssa1, which suggests how Sec72 interacts with full-length Ssa1. Surprisingly, Ssb1, a cytoplasmic Hsp70 that binds ribosome-associated nascent polypeptide chains, also binds to the TPR domain of Sec72, even though it lacks the TPR-binding C-terminal residues of Ssa1. We demonstrate that Ssb1 binds through its ATPase domain to the TPR domain, an interaction that leads to inhibition of nucleotide exchange. Taken together, our results suggest that translocation substrates can be recruited to the Sec71-Sec72 complex either post-translationally through Ssa1 or co-translationally through Ssb1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical timescale effects on detonation convergence

Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.

Adaptive mesh refinement↗

T-1072: Beam test for muon strips

Beam test for muon strips is an experiment to measure the efficiency, time resolution and longitudinal position resolution of scintillator strips for a muon system at future colliders. The muon system at future colliders is envisioned as several layers of position-sensitive detectors embedded in the iron flux-return yoke of the solenoidal magnet. The role of the muon system is primarily the identification of muons and track matching to the central tracker, besides serving as the tail catcher for the hadronic showers that penetrate beyond the hadron calorimetry. Previous tests performed at Fermilab with cosmic muons have shown that scintillators with wavelength-shifting (WLS) fibers and a SiPM readout offer promising performance (NIM A 823, 2016, pp.120-125). The purpose of the present beam tests is to determine the performance parameters with better accuracy. The measurement setup is schematically presented in Fig. 1. The length of the tested strips will vary between 1 m and 2 m. The scintillator strip to be tested will be placed horizontally perpendicular to the beam on a platform in the area behind the beam dump of MT6.2. Two fast scintillators S1 and S2 with vacuum PMT readout will be placed in front and after the tested strip in the beam. The distance between S1 and S2 is 1 m. One fast scintillator strip S3 will be placed in beam immediately in front of the tested strip, perpendicular to the tested strip. A rack with a CAMAC crate and two NIM crates will be used for the data acquisition. The data collection will be performed and monitored from a PC with USB connection to the CAMAC Crate controller.

43 PARTICLE ACCELERATORS↗

Time-Variable Gravity Signal due to Extratropic Pacific Water Mass Redistribution

Using the satellite-laser-ranging (SLR) data, Cox and Chao [2002] reported the detection of a large post-1998 anomaly (in the form of a positive jump) in the time series of Earth s lowest-degree gravity harmonic 52, or the dynamic oblateness. Among several groups now examining the mass redistribution in the global geophysical fluids in search of the cause(s), we report here a temporally coinciding anomalies found in the extratropic north + south Pacific basins. Clearly seen in the leading EOFPC mode for extratropic Pacific, these anomalies occurred in sea-surface height, sea-surface temperature, and temperature- and salinity-depth profiles. We based our analysis on two different data sources: TOPEX/Poseidon altimetry, and the ECCO ocean general circulation model output assimilating T/P data. The magnitude of these changes, when converted to equivalent J2 change, appears to be a few times too small to explain the observed J2 directly. These findings, and the fact that the anomalies occurred following the strong 1997-98 El Nino, suggest strong geophysical connection of the interannual-to-decadal variation of 52 with the Pacific Decadal Oscillation (PDO) and the ultimate global-change processes that cause PDO. More work is underway, and additional independent data sources are examined, paying close attention to the fact that the J2 anomaly has been reversing back to normal since 2001. These include: (1) cryospheric contributions (melting of glaciers and ice sheets); (2) land hydrological contributions; (3) polar sea influences ( e g , via deep flow); (4) fluid flow in Earth's core; (5) time-variable gravity signals from SLR in higher harmonic degree/order, including J3,J4, (2,1), and (2,2) coefficients, considering their lower signal-to-noise ratios; (6) Earth rotation data in terms of length-of-day and polar motion.

Chao, B. F.↗

Bridging Atomic Solvation Environment with Electrochemical Properties for the Bis(trifluoromethylsulfonyl)imide-Based Divalent Cation Electrolytes for the Next-Generation Energy Storage Systems

A deep molecular-level understanding of the multivalent electrolyte and its correlation with the electrochemical properties is crucial for designing optimized electrolytes for next-generation rechargeable batteries. Comprehensive knowledge of the atomic level of the solvation structure and its connection with electrochemical stability and ion transport properties is especially critical. However, the interaction of these three components coupled with clear atomistic insights is lacking in the literature. Here, our current contribution evaluates representative electrolytes with the bis(trifluoromethanesulfonyl)imide (TFSI) anions for multivalent cations of Mg, Ca, and Zn, at different ionic conditions with and without a cosolvated environment in ether-based solvent. Two critical problems are investigated: first, resolving the solvation structures in the electrolyte solutions as a function of concentrations through pair distribution function analysis and the corresponding electrochemical transport properties; second, unmasking the quantitative correlation of the atomistic environment with both electrochemical kinetics and cation dependence. We discovered that the magnesium- and calcium-based electrolytes display versatile coordination lengths but poor average anodic stability due to ion pairing with TFSI - . On the contrary, the zinc-based electrolytes show the shortest solvent coordination lengths, shielding the Zn cation from rigid solvent interactions and resulting in the highest anodic stabilities. Calcium-based electrolytes exhibit the longest and most concentration-independent coordination lengths. This work provides valuable insights into the molecular structural and electrochemical features of diverse multivalent electrolyte systems with cations in various solvation environments, emphasizing the importance of the solvation structure and construction in designing high-performance electrolytes.

cation coordination↗

Designing Peptide Fossils That Model the Evolution of the Bacterial Ferredoxin Fold

Electron transfer coupled to redox chemistry is at the heart of metabolism. The proteins responsible for moving electrons (protein electron carriers) must have emerged at the origin of life. The small iron–sulfur-binding bacterial ferredoxins were likely among these first proteins. Embedded within the ferredoxin sequence and structure is a symmetry that points to an ancient gene duplication event. Little is understood about the nature of ferredoxins prior to this duplication event or what environmental factors may have driven the selection for more complex forms. The deep-time molecular history of ferredoxins goes back billions of years and cannot be reconstructed by phylogenetic analyses based on amino acid sequences. Here, we use structure-guided protein design to model a fossil half-ferredoxin stage in the evolution of this fold, the semidoxins, and their symmetric full-length counterparts, the symdoxins. Semidoxin designs homodimerize, exhibiting structural, thermodynamic, and electrochemical behaviors in most cases identical to cognate symdoxins. However, the semi- and symdoxin fossil stages behave differently when incorporated into an in vivo electron transfer complementation assay. Both can support bacterial growth dependent on protein expression. Growth rates of bacteria expressing the semidoxins are much more sensitive to oxygen than those of bacteria expressing symdoxins. Motivated by the in vivo functionality of designed semidoxins, we identified putative naturally occurring semidoxins in extant anaerobic microorganisms. This is consistent with the observed in vivo oxygen sensitivity of the semidoxin designs. One natural semidoxin is shown to be folded and redox active. However, it exists as a mixture of monomers and dimers, suggesting a potential connection between semidoxins and even simpler single iron–sulfur cluster-binding peptides.

59 BASIC BIOLOGICAL SCIENCES↗

Understanding interfacial segregation in polymer blend films with random and mixed side chain bottlebrush copolymer additives

Bottlebrush polymers are complex macromolecules with tunable physical properties dependent on the chemistry and architecture of both the side chains and the backbone. Prior work has demonstrated that bottlebrush polymer additives can be used to control the interfacial properties of blends with linear polymers but has not specifically addressed the effects of bottlebrush side chain microstructures. Here, using a combination of experiments and self-consistent field theory (SCFT) simulations, we investigated the effects of side chain microstructures by comparing the segregation of bottlebrush additives having random copolymer side chains with bottlebrush additives having a mixture of two different homopolymer side chain chemistries. Specifically, we synthesized bottlebrush polymers with either poly(styrene-ran-methyl methacrylate) side chains or with a mixture of polystyrene (PS) and poly(methyl methacrylate) (PMMA) side chains. Furthermore, the bottlebrush additives were matched in terms of PS and PMMA compositions, and they were blended with linear PS or PMMA chains that ranged in length from shorter to longer than the bottlebrush side chains. Experiments revealed similar behaviors of the two types of bottlebrushes, with a slight preference for mixed side-chain bottlebrushes at the film surface. SCFT simulations were qualitatively consistent with experimental observations, predicting only slight differences in the segregation of bottlebrush additives driven by side chain microstructures. Specifically, these slight differences were driven by the chemistries of the bottlebrush polymer joints and side chain end-groups, which were entropically repelled and attracted to interfaces, respectively. Using SCFT, we also demonstrated that the interfacial behaviors were dominated by entropic effects with high molecular weight linear polymers, leading to enrichment of bottlebrush near interfaces. Surprisingly, the SCFT simulations showed that the chemistry of the joints connecting the bottlebrush backbones and side chains played a more significant role compared with the side chain end groups in affecting differences in surface excess of bottlebrushes with random and mixed side chains. This work provides new insights into the effects of side chain microstructure on segregation of bottlebrush polymer additives.

36 MATERIALS SCIENCE↗

Factors in the Design of Centrifugal Type Injection Valves for Oil Engines

This research was undertaken in connection with a general study of the application of the fuel injection engine to aircraft. The purpose of the investigation was to determine the effect of four important factors in the design of a centrifugal type automatic injection valve on the penetration, general shape, and distribution of oil sprays. The general method employed was to record the development of single sprays by means of special high-speed photographic apparatus capable of taking 25 consecutive pictures of the moving spray at a rate of 4,000 per second. Investigations were made concerning the effects on spray characteristics, of the helix angle of helical grooves, the ratio of the cross-sectional area of the orifice to that of the grooves, the ratio of orifice length to diameter, and the position of the seat. Maximum spray penetration was obtained with a ratio of orifice length to diameter of about 1.5. Slightly greater penetration was obtained with the seat directly before the orifice.

Joachim, W F↗

Thermal compensator for closed-cycle helium refrigerator

The wave length of an infrared, semiconductor laser diode having an output frequency that is dependent on the diode temperature is maintained substantially constant by maintaining the diode temperature constant. The diode is carried by a cold tip of a closed cycle helium refrigerator. The refrigerator has a tendency to cause the temperature of the cold tip to oscillate. A heater diode and a sensor diode are placed on a thermal heat sink that is the only highly conductive thermal path between the laser diode and the cold tip. The heat sink has a small volume and low thermal capacitance so that the sensing diode is at substantially the same temperature as the heater diode and substantially no thermal lag exists between them. The sensor diode is connected in a negative feedback circuit with the heater diode so that the tendency of the laser diode to thermally oscillate is virtually eliminated.

Jennings, D. E.↗

Extending unbiased stereology of brain ultrastructure to three-dimensional volumes

OBJECTIVE: Analysis of brain ultrastructure is needed to reveal how neurons communicate with one another via synapses and how disease processes alter this communication. In the past, such analyses have usually been based on single or paired sections obtained by electron microscopy. Reconstruction from multiple serial sections provides a much needed, richer representation of the three-dimensional organization of the brain. This paper introduces a new reconstruction system and new methods for analyzing in three dimensions the location and ultrastructure of neuronal components, such as synapses, which are distributed non-randomly throughout the brain. DESIGN AND MEASUREMENTS: Volumes are reconstructed by defining transformations that align the entire area of adjacent sections. Whole-field alignment requires rotation, translation, skew, scaling, and second-order nonlinear deformations. Such transformations are implemented by a linear combination of bivariate polynomials. Computer software for generating transformations based on user input is described. Stereological techniques for assessing structural distributions in reconstructed volumes are the unbiased bricking, disector, unbiased ratio, and per-length counting techniques. A new general method, the fractional counter, is also described. This unbiased technique relies on the counting of fractions of objects contained in a test volume. A volume of brain tissue from stratum radiatum of hippocampal area CA1 is reconstructed and analyzed for synaptic density to demonstrate and compare the techniques. RESULTS AND CONCLUSIONS: Reconstruction makes practicable volume-oriented analysis of ultrastructure using such techniques as the unbiased bricking and fractional counter methods. These analysis methods are less sensitive to the section-to-section variations in counts and section thickness, factors that contribute to the inaccuracy of other stereological methods. In addition, volume reconstruction facilitates visualization and modeling of structures and analysis of three-dimensional relationships such as synaptic connectivity.

NASA Discipline Neuroscience↗

Interplanetary Magnetic Field Guiding Relativistic Particles

The origin and the propagation of relativistic solar particles (0.5 to few Ge V) in the interplanetary medium remains a debated topic. These relativistic particles, detected at the Earth by neutron monitors have been previously accelerated close to the Sun and are guided by the interplanetary magnetic field (IMF) lines, connecting the acceleration site and the Earth. Usually, the nominal Parker spiral is considered for ensuring the magnetic connection to the Earth. However, in most GLEs the IMF is highly disturbed, and the active regions associated to the GLEs are not always located close to the solar footprint of the nominal Parker spiral. A possible explanation is that relativistic particles are propagating in transient magnetic structures, such as Interplanetary Coronal Mass Ejections (ICMEs). In order to check this interpretation, we studied in detail the interplanetary medium where the particles propagate for 10 GLEs of the last solar cycle. Using the magnetic field and the plasma parameter measurements (ACE/MAG and ACE/SWEPAM), we found widely different IMF configurations. In an independent approach we develop and apply an improved method of the velocity dispersion analysis to energetic protons measured by SoHO/ERNE. We determined the effective path length and the solar release time of protons from these data and also combined them with the neutron monitor data. We found that in most of the GLEs, protons propagate in transient magnetic structures. Moreover, the comparison between the interplanetary magnetic structure and the interplanetary length suggest that the timing of particle arrival at Earth is dominantly determined by the type of IMF in which high energetic particles are propagating. Finally we find that these energetic protons are not significantly scattered during their transport to Earth.

Masson, S.↗

RivGraph: Automatic extraction and analysis of river and delta channel network topology

River networks sustain life and landscapes by carrying and distributing water, sediment, and nutrients throughout ecosystems and communities. At the largest scale, river networks drain continents through tree-like tributary networks. At typically smaller scales, river deltas and braided rivers form loopy, complex distributary river networks via avulsions and bifurcations.In order to model flows through these networks or analyze network structure, the topology, or connectivity, of the network must be resolved. Additionally, morphologic properties of each river channel as well as the direction of flow through the channel inform how fluxes travel through the network’s channels. Riv Graphis a Python package that automates the extraction and characterization of river channel networks from a user-provided binary image, or mask, of a channel network (Fig. 1). Masks may be derived from (typically remotely-sensed) imagery, simulations, or even hand-drawn. RivGraph will create explicit representations of the channel network by resolving river centerlines as links, and junctions as nodes. Flow directions are solved for each link of the network without using auxiliary data, e.g., a digital elevation model (DEM). Morphologic properties are computed as well, including link lengths, widths, sinuosities, branching angles,and braiding indices. If provided,RivGraph will preserve georeferencing information of the mask and will export results as ESRI shapefiles, GeoJSONs, and GeoTIFFs for easy import into GIS software.RivGraph can also return extracted networks as networkx objects for convenient interfacing with the full-featured networkx package (Hagberg et al., 2008). Finally, RivGraph offers a suite of topologic metrics that were specifically designed for river channel network analysis (Tejedor et al., 2015b).

54 ENVIRONMENTAL SCIENCES↗

Penalty-Based Interface Technology for Prediction of Delamination Growth in Laminated Structures

An effective interface element technology has been developed for connecting and simulating crack growth between independently modeled finite element subdomains (e.g., composite plies). This method has been developed using penalty constraints and allows coupling of finite element models whose nodes do not necessarily coincide along their common interface. Additionally, the present formulation leads to a computational approach that is very efficient and completely compatible with existing commercial software. The present interface element has been implemented in the commercial finite element code ABAQUS as a User Element Subroutine (UEL), making it easy to test the approach for a wide range of problems. The interface element technology has been formulated to simulate delamination growth in composite laminates. Thanks to its special features, the interface element approach makes it possible to release portions of the interface surface whose length is smaller than that of the finite elements. In addition, the penalty parameter can vary within the interface element, allowing the damage model to be applied to a desired fraction of the interface between the two meshes. Results for double cantilever beam DCB, end-loaded split (ELS) and fixed-ratio mixed mode (FRMM) specimens are presented. These results are compared to measured data to assess the ability of the present damage model to simulate crack growth.

Averill, Ronald C.↗

Flow Channeling in Fracture Networks: Characterizing the Effect of Density on Preferential Flow Path Formation

Flow channelization is a commonly observed phenomenon in fractured subsurface media where the flow of fluids is restricted primarily to highly transmissive fracture networks surrounded by a low-permeability rock matrix. The multiscale structural heterogeneity of these networks results in multiscale flow channelization where preferential flow paths form at length scales ranging from the entire system down to the subfracture size. We present an analysis of how one of the largest scales in fractured media, the network density, influences the degree of flow channeling that occurs using an ensemble of semigeneric three-dimensional discrete fracture network (DFN) simulations. We construct 10 DFNs, whose fracture lengths follow a power law distribution, at four densities for a total of 40 networks. We characterize their structure in terms of the network topology and geometry. Eulerian and Lagrangian observations of the steady-state flow fields obtained within the networks are used to quantify the degree of flow channelization at the network scale. We introduce a measure for the importance-ranking/hierarchy of different flow paths in the network using graph-based analysis of Lagrangian transport by which the degree of flow channeling between networks is compared. These flow observations are then linked to the structural properties of the networks. In general, network-scale flow channeling decreases as the network density increases. However, at low densities, there is more uniform flow within the entire connected network than in high-density networks. We also demonstrate how standard transport observables can be used to infer the degree of flow channelization occurring within a fracture network.

58 GEOSCIENCES↗