Search NASA⌕ Search

SEARCH · Search NASA

Results for “small structures”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Indirectly probing the structural change in ion-irradiated Zr-Based metallic glasses from small scale mechanical tests

Ion irradiation was applied to tailor the structural heterogeneities in Zr-based metallic glasses at room temperature. Experimental methods of X-ray diffraction, nanoindentation, and micropillar compression were conducted to examine the irradiation effects on their structural and mechanical property changes. It is found that the irradiated materials retained amorphous structure after room-temperature Ni ion irradiation. The reduction of elastic modulus and hardness measured by nanoindentation indicated the irradiation-induced mechanical degradation. A unified statistic model was employed to quantitatively predict the density and strength of irradiation defects, although their specific structural and physical nature is not explicitly included in this model. The transition of non-intersecting shear bands to multiple intersecting shear bands was observed on compression tests of irradiated micropillars with the increase of irradiation dose. Furthermore, the analysis of the displacement excursion of micropillar compression tests indicated a different deformation mode from the unirradiated state. These results from nanoindentation pop-in and micro-pillar compression tests suggested that irradiation eventually leads to a new state with different types and characteristics of structural heterogeneities from violent displacement cascades and non-equilibrium energy deposition/dissipation processes, which also proves ion irradiation as an effective method to tune the structure and mechanical properties of metallic glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A twist grain boundary phase in aqueous solutions of the nucleic acid tetramer GTAC

At high concentration, long Watson/Crick (WC) double-helixed DNA forms columnar crystal or liquid crystal phases of linear, parallel duplex chains packed on periodic lattices. This can also be a structural motif of short NA oligomers such as the 5’-GTAC-3’ studied here, which makes four-base WC duplexes having hydrophobic blunt ends. End-to-end aggregation then assembles these duplexes into columns and columnar phases are stabilized, in spite of breaks in the double helix every four bases. But the new degrees of freedom introduced by such breaks also enable opportunities for a more diverse palette of self-assembly modes, producing striking self-assemblies of DNA that would not be achievable with contiguous polymers. These include recently reported three-dimensional (3D) periodic low-density nanoscale networks of GCCG, and the twist grain boundary (TGB) phase presented here. In the TGB, columns of GTAC pairs assemble into monolayer sheets in which the duplex columns are mutually parallel. However, unlike in the columnar crystals, these sheets stack in helical fashion into lamellar arrays in which the column axis of each layer is rotated through a 60° angle with respect to the columns in neighboring layers. This assembly of DNA is unique in that it the fills a 3D volume wherein the major grooves of columns in each layer mutually enter and interlock with the major grooves of columns in neighboring layers. This locking is optimized by small adjustments in structure enabled by the breaks in the duplex backbones.

DNA↗
DOE OSTI2023

Exceptional thermal stability of additively manufactured CoCrFeMnNi high-entropy alloy with cellular dislocation structures

CoCrFeMnNi high-entropy alloy (HEA) was additively manufactured (AM) by laser powder-bed fusion (L-PBF). The AM CoCrFeMnNi has prominent cellular dislocation structures with a small number of Mn-rich oxides. The thermal stability of the AM CoCrFeMnNi was investigated by isochronal annealing treatment at various temperatures from 400 to 1300°C for 1h. Microstructural analysis shows slow dislocation recovery, retarded recrystallization process, and precipitation of additional Cr-Mn based oxides during thermal annealing, resulting in exceptional thermal stability and retained high hardness at elevated temperatures. Further, by thermodynamic calculations, a low stored energy of 1.31 MJ/m 3 and a high activation energy of 353 kJ/mol for recrystallization were estimated for the AM CoCrFeMnNi. The exceptional thermal stability of the AM CoCrFeMnNi HEA is mechanistically attributed to the low crystallographic misorientations across the dislocation cell walls, sluggish atomic diffusion, and the pinning effects of the oxide nanoprecipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Position-dependent power spectrum of the large-scale structure: a novel method to measure the squeezed-limit bispectrum

The influence of large-scale density fluctuations on structure formation on small scales is described by the three-point correlation function (bispectrum) in the so-called ''squeezed configurations,'' in which one wavenumber, say k{sub 3}, is much smaller than the other two, i.e., k{sub 3} << k{sub 1} ≈ k{sub 2}. This bispectrum is generated by non-linear gravitational evolution and possibly also by inflationary physics. In this paper, we use this fact to show that the bispectrum in the squeezed configurations can be measured without employing three-point function estimators. Specifically, we use the ''position-dependent power spectrum,'' i.e., the power spectrum measured in smaller subvolumes of the survey (or simulation box), and correlate it with the mean overdensity of the corresponding subvolume. This correlation directly measures an integral of the bispectrum dominated by the squeezed configurations. Measuring this correlation is only slightly more complex than measuring the power spectrum itself, and sidesteps the considerable complexity of the full bispectrum estimation. We use cosmological N-body simulations of collisionless particles with Gaussian initial conditions to show that the measured correlation between the position-dependent power spectrum and the long-wavelength overdensity agrees with the theoretical expectation. The position-dependent power spectrum thus provides a new, efficient, and promising way to measure the squeezed-limit bispectrum from large-scale structure observations such as galaxy redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Small molecule inhibitors of mesotrypsin from a structure-based docking screen

PRSS3/mesotrypsin is an atypical isoform of trypsin, the upregulation of which has been implicated in promoting tumor progression. To date there are no mesotrypsin-selective pharmacological inhibitors which could serve as tools for deciphering the pathological role of this enzyme, and could potentially form the basis for novel therapeutic strategies targeting mesotrypsin. A virtual screen of the Natural Product Database (NPD) and Food and Drug Administration (FDA) approved Drug Database was conducted by high-throughput molecular docking utilizing crystal structures of mesotrypsin. Twelve high-scoring compounds were selected for testing based on lowest free energy docking scores, interaction with key mesotrypsin active site residues, and commercial availability. Diminazene (C1D22956468), along with two similar compounds presenting the bis-benzamidine substructure, was validated as a competitive inhibitor of mesotrypsin and other human trypsin isoforms. Diminazene is the most potent small molecule inhibitor of mesotrypsin reported to date with an inhibitory constant (K i ) of 3.6±0.3 pM. Diminazene was subsequently co-crystalized with mesotrypsin and the crystal structure was solved and refined to 1.25 Å resolution. This high resolution crystal structure can now offer a foundation for structure-guided efforts to develop novel and potentially more selective mesotrypsin inhibitors based on similar molecular substructures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2010

Numerical computations and optical diagnostics of unsteady partially premixed methane/air flames

The structures and dynamics of unsteady laminar partially premixed methane/air Bunsen flames are studied by means of numerical simulations, OH and CH PLIF imaging, and high speed chemiluminescence imaging employing a high framing speed intensified charge coupled device camera. The Bunsen burner has a diameter of 22 mm. Rich methane/air mixtures with an equivalence ratio of 1.5 are injected from the burner into atmosphere at different flow speeds ranging from 0.77 to 1.7 m/s, with Reynolds numbers based on the nozzle flow ranging from 1100 to 2500. The numerical simulations are based on a two-scalar flamelet manifold tabulation approach. Detailed chemistry is used to generate the flamelet manifold tabulation which relates the species concentrations, reaction rates, temperature and density to a distance function G and mixture fraction Z. Two distinct reaction zones are identified using CH and OH PLIF imaging and numerical simulations; one inner reaction zone corresponds to premixed flames on the rich side of the mixture and one outer reaction zone corresponds to mixing controlled diffusion flames on the lean side of the mixture. Under normal gravity conditions both the inner premixed flames and the outer diffusion flames are unsteady. The outer diffusion flames oscillate with a flickering frequency of about 15 Hz, which slightly increases with the burner exit velocity. The inner premixed flames are more random with much more small-scale wrinkling structures. Under zero gravity conditions the outer diffusion flames are stable whereas the inner premixed flames are unstable and highly wrinkled. It appears that the outer diffusion flames are governed by the Rayleigh-Taylor instability whereas the inner premixed flames are dictated by Landau-Darrieus instability. The two-scalar flamelet approach is shown to capture the basic structures and dynamics of the investigated unsteady partially premixed flames. (author)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Observation of spin textures in La{sub 1−x}Sr{sub x}MnO{sub 3} (x = 0.175)

We have investigated topological spin textures in the ferromagnetic metallic phase of La{sub 0.825}Sr{sub 0.175}MnO{sub 3} with the centrosymmetric crystal structure by small-angle electron diffraction (SmAED) and low-temperature Lorentz transmission electron microscopy (TEM) experiments. In-situ Lorentz TEM and SmAED experiments revealed that type-I and type-II magnetic bubbles evolved from magnetic stripe domains with the Bloch-type domain wall by applying vertical magnetic field. Type-I magnetic bubbles with left-handed and right-handed spin helicity were randomly distributed and simultaneously type-II magnetic bubbles are formed locally. The important point about type-I and type-II magnetic bubbles is that their emergence depends strongly on whether perpendicular magnetic field is applied parallel to the magnetic easy axis along the [001] direction. Our experimental results suggested that the stabilization of magnetic bubbles should originate from the long-range dipole-dipole interactions, as opposed to the Dzyaloshinskii-Moriya interaction in helical magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗
DOE OSTI2014

Dephosphorylation of the Core Clock Protein KaiC in the Cyanobacterial KaiABC Circadian Oscillator Proceeds via an ATP Synthase Mechanism

The circadian clock of the cyanobacterium Synechococcus elongatus can be reconstituted in vitro from three proteins, KaiA, KaiB, and KaiC in the presence of ATP, to tick in a temperature-compensated manner. KaiC, the central cog of this oscillator, forms a homohexamer with 12 ATP molecules bound between its N- and C-terminal domains and exhibits unusual properties. Both the N-terminal (CI) and C-terminal (CII) domains harbor ATPase activity, and the subunit interfaces between CII domains are the sites of autokinase and autophosphatase activities. Hydrolysis of ATP correlates with phosphorylation at threonine and serine sites across subunits in an orchestrated manner, such that first T432 and then S431 are phosphorylated, followed by dephosphorylation of these residues in the same order. Although structural work has provided insight into the mechanisms of ATPase and kinase, the location and mechanism of the phosphatase have remained enigmatic. From the available experimental data based on a range of approaches, including KaiC crystal structures and small-angle X-ray scattering models, metal ion dependence, site-directed mutagenesis (i.e., E318, the general base), and measurements of the associated clock periods, phosphorylation patterns, and dephosphorylation courses as well as a lack of sequence motifs in KaiC that are typically associated with known phosphatases, we hypothesized that KaiCII makes use of the same active site for phosphorylation and dephosphorlyation. We observed that wild-type KaiC (wt-KaiC) exhibits an ATP synthase activity that is significantly reduced in the T432A/S431A mutant. We interpret the first observation as evidence that KaiCII is a phosphotransferase instead of a phosphatase and the second that the enzyme is capable of generating ATP, both from ADP and P{sub i} (in a reversal of the ATPase reaction) and from ADP and P-T432/P-S431 (dephosphorylation). This new concept regarding the mechanism of dephosphorylation is also supported by the strikingly similar makeups of the active sites at the interfaces between {alpha}/{beta} heterodimers of F1-ATPase and between monomeric subunits in the KaiCII hexamer. Several KaiCII residues play a critical role in the relative activities of kinase and ATP synthase, among them R385, which stabilizes the compact form and helps kinase action reach a plateau, and T426, a short-lived phosphorylation site that promotes and affects the order of dephosphorylation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Structural modulation, physical properties, and electronic band structure of the kagome metal UCr 6 ⁢Ge 6

The chemical flexibility of the 𝑅⁢𝑀 6 ⁢𝑋 6 stoichiometry, where an 𝑓-block element is intercalated in the CoSn structure type, allows for the tuning of flatbands associated with kagome lattices to the Fermi level and for emergent phenomena due to interactions between the 𝑓- and 𝑑-electron lattices. Yet, 5⁢𝑓 members of the “166” compounds are underrepresented compared with 4⁢𝑓 members. In this work, we report single-crystal growth of UCr 6 ⁢Ge 6 , which crystallizes in a monoclinically distorted Y 0.5⁢ Co 3 ⁢Ge 3 -type structure. The real-space character of the modulation, which is unique within the 𝑅⁢𝑀 6 ⁢𝑋 6 family, is approximated by a 3×1×2 supercell of the average monoclinic cell. The compound has kagome-lattice flatbands near the Fermi level and a moderately enhanced electronic heat capacity, as evidenced by its low-temperature Sommerfeld coefficient (𝛾=86.5 mJ mol −1 K −2 ) paired with band structure calculations. The small, isotropic magnetization and featureless resistivity of UCr 6 ⁢Ge 6 suggest itinerant uranium 5⁢𝑓 electrons and Pauli paramagnetism. Angle-resolved photoemission spectroscopy results provide evidence for uranium 5⁢𝑓 weight at the Fermi level and for a flatband near the Fermi level associated with the chromium 3⁢𝑑 kagome lattice. The isotropic magnetic behavior of the uranium 5⁢𝑓 electrons starkly contrasts with localized behavior in other uranium 166 compounds, highlighting the high tunability of the magnetic ground state across the material family.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparison of the static structure factor at long wavelengths for a dusty plasma liquid and other liquids

Especially small values of the static structure factor S(k) at long wavelengths, i.e., small k, were obtained in an analysis of experimental data, for a two-dimensional dusty plasma in its liquid state. For comparison, an analysis of S(k) data was carried out for many previously published experiments with other liquids. Here, the latter analysis indicates that the magnitude of S(k) at small k is typically in a range 0.02–0.13. In contrast, the corresponding value for a dusty plasma liquid was found to be as small as 0.0139. Another basic finding for the dusty plasma liquid is that S(k) at small k generally increases with temperature, with its lowest value, noted above, occurring near the melting point. Simulations were carried out for the dusty plasma liquid, and their results are generally consistent with the experiment. Since a dusty plasma has a soft interparticle interaction, our findings support earlier theoretical suggestions that a useful design strategy for creating materials having exceptionally low values of S(0), so-called hyperuniform materials, is the use of a condensed material composed of particles that interact softly at their periphery.

2-dimensional systems↗
DOE OSTI2018

A generative graph model for electrical infrastructure networks

We propose a generative graph model for electrical infrastructure networks that accounts for heterogeneity in both node and edge type. To inform the design of this model, we analyze the properties of power grid graphs derived from the U.S. Eastern Interconnection, Texas Interconnection, and Poland transmission system power grids. Across these datasets, we find subgraphs induced by nodes of the same voltage level exhibit shared structural properties atypical to small-world networks, including low local clustering, large diameter and large average distance. On the other hand, we find subgraphs induced by transformer edges linking nodes of different voltage types contain a more limited structure, consisting mainly of small, disjoint star graphs. The goal of our proposed model is to match both these inter and intra-network properties by proceeding in two phases: we first generate subgraphs for each voltage level and then generate transformer edges that connect these subgraphs. The first phase of the model adapts the Chung-Lu random graph model, taking desired vertex degrees and desired diameter as inputs, while the second phase of the model is based on a simpler random star graph generation process. We test the performance of our model by comparing its output across many runs to the aforementioned real data. In nearly all categories tested, we find our model is more accurate in reproducing the unusual mixture of properties apparent in the data than the Chung-Lu model. We also include graph visualization comparisons, as well as a brief analysis of edge-deletion resiliency. Ultimately, our model may be used to generate synthetic graph data, test hypotheses and algorithms at different scales, and serve as a baseline model on top of which further electrical network properties, device models, and interdependencies to other networks, may be appended.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2013

Edible oil structures at low and intermediate concentrations. II. Ultra-small angle X-ray scattering of in situ tristearin solids in triolein

Ultra-small angle X-ray scattering has been used for the first time to elucidate, in situ, the aggregation structure of a model edible oil system. The three-dimensional nano- to micro-structure of tristearin solid particles in triolein solvent was investigated using 5, 10, 15, and 20% solids. Three different sample preparation procedures were investigated: two slow cooling rates of 0.5°/min, case 1 (22 days of storage at room temperature) and case 2 (no storage), and one fast cooling of 30°/min, case 3 (no storage). The length scale investigated, by using the Bonse-Hart camera at beamline ID-15D at the Advanced Photon Source, Argonne National Laboratory, covered the range from 300 Å to 10 μm. The unified fit and the Guinier-Porod models in the Irena software were used to fit the data. The former was used to fit 3 structural levels. Level 1 structures showed that the primary scatterers were essentially 2-dimensional objects for the three cases. The scatterers possessed lateral dimensions between 1000 and 4300 Å. This is consistent with the sizes of crystalline nanoplatelets present which were observed using cryo-TEM. Level 2 structures were aggregates possessing radii of gyration, R{sub g2} between 1800 Å and 12000 Å and fractal dimensions of either D{sub 2}=1 for case 3 or 1.8≤D{sub 2}≤2.1 for case 1 and case 2. D{sub 2} = 1 is consistent with unaggregated 1-dimensional objects. 1.8 ≤ D{sub 2} ≤ 2.1 is consistent with these 1-dimensional objects (below) forming structures characteristic of diffusion or reaction limited cluster-cluster aggregation. Level 3 structures showed that the spatial distribution of the level 2 structures was uniform, on the average, for case 1, with fractal dimension D{sub 3}≈3 while for case 2 and case 3 the fractal dimension was D{sub 3}≈2.2, which suggested that the large-scale distribution had not come to equilibrium. The Guinier-Porod model showed that the structures giving rise to the aggregates with a fractal dimension given by D{sub 2} in the unified fit level 2 model were cylinders described by the parameter s≈1 in the Guinier-Porod model. The size of the base of these cylinders was in agreement with the cryo-TEM observations as well as with the results of the level 1 unified fit model. By estimating the size of the nanoplatelets and understanding the structures formed via their aggregation, it will be possible to engineer novel lipids systems that embody desired functional characteristics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The C‐Terminal Domain of α‐Synuclein Confers Steric Stabilization on Synaptic Vesicle‐Like Surfaces

Abstract While α‐synuclein, an intrinsically disordered protein linked to Parkinson's disease, has been shown to associate with membrane organelles, its overall cellular function remains nebulous. α‐Synuclein binds to membranes through its amino‐terminal domain (first ≈100 residues), but there is no consensus on the biophysical function of the carboxyl‐terminal domain (last ≈40 residues) due, in part, to its lack of strong interaction partners and persisting intrinsic disorder even when membrane bound. Here, by directly applying force on α‐synuclein bound to spherical nanoparticle‐supported lipid bilayers (SSLBs) and tracking higher‐order structural changes through small‐angle X‐ray scattering, strong evidence is presented that α‐synuclein sterically stabilizes membrane surfaces through its carboxyl‐terminal domain. Full‐length α‐synuclein dramatically increases the critical osmotic pressure at which SSLBs cluster ( P C ≈ 1.3 × 10 5 Pa) compared to α‐synuclein without the carboxyl‐terminal domain ( P C ≈ 1.9 × 10 4 Pa) at physiological salt and temperature conditions. This clustering of α‐synuclein‐bound SSLBs is shown to be reversible and sensitive to monovalent/divalent salt, both features of grafted polyelectrolyte‐mediated steric stabilization. In elucidating the biophysical function of α‐synuclein in the framework of polymer science, it is demonstrated that the carboxyl‐terminal domain can potentially utilize its persisting intrinsic disorder to functionalize membrane surfaces.

Chung, Peter J.↗
DOE OSTI2026

Effect of Thermodynamic and Environmental Factors on Crystallization of DNA‐Origami Superlattices

The directed self‐assembly of nanoscale materials into ordered superlattices presents a powerful strategy for creating next‐generation materials with programmable mechanical, optical, and photonic properties. Deoxyribonucleic acid (DNA) origami has emerged as a versatile scaffold for encoding nanoscale geometry and guiding the crystallization of complex 3D architectures. However, a systematic understanding of the parameters that govern the efficiency and quality of superlattice formation remains limited. In this study, we utilize octahedral DNA nanoscale frames as a model system to investigate the relative influence of key factors, including buffer composition, ionic strength, frame concentration, and thermal annealing protocols, on the size, order, and reproducibility of the resulting superlattices. Our findings provide a quantitative framework to rationally optimize DNA‐based assembly pathways. Structural characterization via small‐angle x‐ray scattering (SAXS), scanning electron microscopy (SEM), and optical microscopy validates the quality and fidelity of the assembled lattices. Moreover, by templating these DNA frameworks into inorganic replicas, we establish general design principles that extend beyond biomolecular systems, providing a foundation for the synthesis of programmable materials in broader nanofabrication contexts.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Analytic atomic gradients in the fermi‐löwdin orbital self‐interaction correction

We derived, implemented, and thoroughly tested the complete analytic expression for atomic forces, consisting of the Hellmann‐Feynman term and the Pulay correction, for the Fermi‐Löwdin orbital self‐interaction correction (FLO‐SIC) method. Analytic forces are shown to be numerically accurate through an extensive comparison to forces obtained from finite differences. Using the analytic forces, equilibrium structures for a small set of molecules were obtained. This work opens the possibility of routine self‐interaction free geometrical relaxations of molecules using the FLO‐SIC method. © 2018 Wiley Periodicals, Inc.

Trepte, Kai↗
DOE OSTI2017

Inverse modeling and uncertainty analysis of potential groundwater recharge to the confined semi-fossil Ohangwena II Aquifer, Namibia

The identification of potential recharge areas and estimation of recharge rates to the confined semi-fossil Ohangwena II Aquifer (KOH-2) is crucial for its future sustainable use. The KOH-2 is located within the endorheic transboundary Cuvelai-Etosha-Basin (CEB), shared by Angola and Namibia. The main objective was the development of a strategy to tackle the problem of data scarcity, which is a well-known problem in semi-arid regions. In a first step, conceptual geological cross sections were created to illustrate the possible geological setting of the system. Furthermore, groundwater travel times were estimated by simple hydraulic calculations. A two-dimensional numerical groundwater model was set up to analyze flow patterns and potential recharge zones. The model was optimized against local observations of hydraulic heads and groundwater age. The sensitivity of the model against different boundary conditions and internal structures was tested. Parameter uncertainty and recharge rates were estimated. Results indicate that groundwater recharge to the KOH-2 mainly occurs from the Angolan Highlands in the northeastern part of the CEB. The sensitivity of the groundwater model to different internal structures is relatively small in comparison to changing boundary conditions in the form of influent or effluent streams. Uncertainty analysis underlined previous results, indicating groundwater recharge originating from the Angolan Highlands. The estimated recharge rates are less than 1% of mean yearly precipitation, which are reasonable for semi-arid regions.

58 GEOSCIENCES↗
DOE OSTI2015

Particle tracking approach for transport in three-dimensional discrete fracture networks: Particle tracking in 3-D DFNs

The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates mass balance-related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces, and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. As a result, we demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.

58 GEOSCIENCES↗