SEARCH · Engineering Papers

Results for “sim model”

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

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave–particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel & Engelmann ( Phys. Fluids , vol. 9, 1966, pp. 2377–2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that scatter them out of and into resonance, as for Landau damping as shown by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). As a result, the resonant particle–wave interactions occur in the centre of a narrow collisional boundary when the collision frequency $\unicode[STIX]{x1D708}$ is very small compared to the wave frequency $\unicode[STIX]{x1D714}$ . The diffusive nature of the pitch angle scattering combined with the wave–particle resonance condition enhances the collision frequency by $(\unicode[STIX]{x1D714}/\unicode[STIX]{x1D708})^{2/3}\gg 1$ , resulting in an effective resonant particle collisional interaction time of $\unicode[STIX]{x1D70F}_{\text{int}}\sim (\unicode[STIX]{x1D708}/\unicode[STIX]{x1D714})^{2/3}/\unicode[STIX]{x1D708}\ll 1/\unicode[STIX]{x1D708}$ . A collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel–Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

Physics↗

DOE OSTI2025

Scrambling Dynamics with Imperfections in a Solvable Model

We study how probes of quantum scrambling dynamics respond to two kinds of imperfections -- unequal forward and backward evolutions and decoherence -- in a solvable Brownian circuit model. We calculate a ``renormalized'' out-of-time-order correlator (ROTOC) in the model with $N$ qubits, and we show that the circuit-averaged ROTOC is controlled by an effective probability distribution in operator weight space which obeys a system of $N$ non-linear equations of motion. These equations can be easily solved numerically for large system sizes which are beyond the reach of exact methods. Moreover, for an operator initially concentrated on weight one $w_0=1$, we provide an exact solution to the equations in the thermodynamic limit of many qubits that is valid for all times, all non-vanishing perturbation strengths $p\gtrsim 1/\sqrt{N}$, and all decoherence strengths. We also show that a generic initial condition $w_0 >1$ leads to a metastable state that eventually collapses to the $w_0=1$ case after a lifetime $\sim \log(N/w_0)$. Our results highlight situations where it is still possible to extract the unperturbed chaos exponent even in the presence of imperfections, and we comment on the applications of our results to existing experiments with nuclear spins and to future scrambling experiments.

FOS: Physical sciences↗

DOE OSTI2022

LBNL Superfacility Project Report

The Superfacility model is designed to leverage HPC for experimental science. It is more than simply a model of connected experiment, network, and HPC facilities; it encompasses the full ecosystem of infrastructure, software, tools, and expertise needed to make connected facilities easy to use. The three-year Lawrence Berkeley National Laboratory (LBNL) Superfacility project was initiated in 2019 to coordinate work being performed at LBNL to support this model, and to provide a coherent and comprehensive set of science requirements to drive existing and new work. A key component of the project was the in-depth engagements with eight science teams that represent challenging use cases across the DOE Office of Science.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

The predicted properties of helium-enriched globular cluster progenitors at high redshift

Globular cluster progenitors may have been detected by Hubble Space Telescope, and are predicted to be observable with James Webb Space Telescope (JWST) and ground-based extremely large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched (ΔY = +0.12) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first ~30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being ~0.40 mag brighter at $\lambda \approx 2.0\, {\mu \rm m}$, and to be 0.30-mag redder in the JWST–NIRCam colour (F070W - F200W), and to actually be fainter for $\lambda \lesssim 0.50 \, {\mu \rm m}$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim \!5{{\ \rm per\ cent}}$, though we show that the Lyman-α escape fraction could end up higher for helium-enriched stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

The wave energy converter control competition (WECCCOMP): Wave energy control algorithms compared in both simulation and tank testing

The wave energy control competition established a benchmark problem which was offered as an open challenge to the wave energy system control community. The competition had two stages: In the first stage, competitors used a standard wave energy simulation platform (WEC-Sim) to evaluate their controllers while, in the second stage, competitors were invited to test their controllers in a real-time implementation on a prototype system in a wave tank. The performance function used was based on converted energy across a range of standard sea states, but also included aspects related to economic performance, such as peak/average power, peak force, etc. This paper compares simulated and experimental results and, in particular, examines if the results obtained in a linear system simulation are borne out in reality. Overall, within the scope of the device tested, the range of sea states employed, and the performance metric used, the conclusion is that high-performance WEC controllers work well in practice, with good carry-over from simulation to experimentation. However, the availability of a good WEC mathematical model is deemed to be crucial.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Deep XMM–Newton observations of the most distant SPT-SZ galaxy cluster

ABSTRACT We present results from a 577 ks XMM–Newton observation of SPT-CL J0459–4947, the most distant cluster detected in the South Pole Telescope 2500 square degree (SPT-SZ) survey, and currently the most distant cluster discovered through its Sunyaev–Zel’dovich effect. The data confirm the cluster’s high redshift, z = 1.71 ± 0.02, in agreement with earlier, less precise optical/IR photometric estimates. From the gas density profile, we estimate a characteristic mass of $M_{500}=(1.8\pm 0.2)\times 10^{14}\, {\rm M}_{\odot }$; cluster emission is detected above the background to a radius of $\sim \!2.2\, r_{500}$, or approximately the virial radius. The intracluster gas is characterized by an emission-weighted average temperature of 7.2 ± 0.3 keV and metallicity with respect to Solar of $Z/\, Z_{\odot }=0.37\pm 0.08$. For the first time at such high redshift, this deep data set provides a measurement of metallicity outside the cluster centre; at radii $r\gt 0.3\, r_{500}$, we find $Z/\, Z_{\odot }=0.33\pm 0.17$ in good agreement with precise measurements at similar radii in the most nearby clusters, supporting an early enrichment scenario in which the bulk of the cluster gas is enriched to a universal metallicity prior to cluster formation, with little to no evolution thereafter. The leverage provided by the high redshift of this cluster tightens by a factor of 2 constraints on evolving metallicity models, when combined with previous measurements at lower redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Optimization of the light detection system of the ICARUS detector

The Short Baseline Neutrino (SBN) Program at Fermilab is designed to investigate short-baseline neutrino oscillations and test the hypothesis of sterile neutrinos, motivated by several experimental anomalies observed over the past decades. Within this program, the ICARUS experiment plays a key role. It employs the world’s largest Liquid Argon Time Projection Chamber (LArTPC) and serves as the farthest and most sensitive SBN detector for studying muon and electron neutrino oscillations. A crucial subsystem of the ICARUS detector is the Light Detection System (LDS), which captures the prompt scintillation light produced by neutrino interactions in the 600-ton active liquid Argon volume. This system provides precise timing information that is essential for event reconstruction, the trigger system, and cosmic background rejection. The LDS is composed of 360 Hamamatsu R5912-MOD 8-inch photomultiplier tubes (PMTs), operating under cryogenic conditions ($\sim 87 \ K$) inside the detector’s cryostats. During the detector’s operation at FNAL, a degradation in PMT gain has been observed, attributed to aging under low-temperature conditions. In collaboration with ICARUS teams from INFN Pavia and Catania, I developed an experimental setup to study the temperature-dependent behavior of the PMTs, performing gain measurements both at room temperature and down to $-70°C$ using a climatic chamber at INFN Catania. The results indicate that while the PMTs maintain stable gain at room temperature, a significant and permanent gain reduction occurs at low temperatures. Although $-70°C$ is still warmer than liquid Argon temperatures, the findings clearly demonstrate a gain-dependent performance degradation. The thesis also discusses mitigation strategies implemented in the ICARUS detector to address this issue and presents a simplified model to describe and simulate the observed behavior.

Saia, Clara [Catania U.] (ORCID:0009000464102417)↗

DOE OSTI2022

HiFiAdapterFilt, a memory efficient read processing pipeline, prevents occurrence of adapter sequence in PacBio HiFi reads and their negative impacts on genome assembly

Abstract Background Pacific Biosciences HiFi read technology is currently the industry standard for high accuracy long-read sequencing that has been widely adopted by large sequencing and assembly initiatives for generation of de novo assemblies in non-model organisms. Though adapter contamination filtering is routine in traditional short-read analysis pipelines, it has not been widely adopted for HiFi workflows. Results Analysis of 55 publicly available HiFi datasets revealed that a read-sanitation step to remove sequence artifacts derived from PacBio library preparation from read pools is necessary as adapter sequences can be erroneously integrated into assemblies. Conclusions Here we describe the nature of adapter contaminated reads, their consequences in assembly, and present HiFiAdapterFilt, a simple and memory efficient solution for removing adapter contaminated reads prior to assembly.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Detecting Anomalies in Time Series Using Kernel Density Approaches

This paper introduces a novel anomaly detection approach tailored for time series data with exclusive reliance on normal events during training. Our key innovation lies in the application of kernel-density estimation (KDE) to scrutinize reconstruction errors, providing an empirically derived probability distribution for normal events post-reconstruction. This non-parametric density estimation technique offers a nuanced understanding of anomaly detection, differentiating it from prevalent threshold-based mechanisms in existing methodologies. In post-training, events are encoded, decoded, and evaluated against the estimated density, providing a comprehensive notion of normality. In addition, we propose a data augmentation strategy involving variational autoencoder-generated events and a smoothing step for enhanced model robustness. The significance of our autoencoder-based approach is evident in its capacity to learn normal representation without prior anomaly knowledge. Through the KDE step on reconstruction errors, our method addresses the versatility of anomalies, departing from assumptions tied to larger reconstruction errors for anomalous events. Our proposed likelihood measure then distinguishes normal from anomalous events, providing a concise yet comprehensive anomaly detection solution. The extensive experimental results support the feasibility of our proposed method, yielding significantly improved classification performance by nearly 10% on the UCR benchmark data.

Frehner, Robin↗

DOE OSTI2024

Life-cycle analysis of offshore macroalgae production systems in the United States

Offshore macroalgae production offers the potential to provide valuable biomass for food, energy, and higher value products without the use of land or freshwater while using excess nutrients and carbon dioxide. To realize this potential, the Macroalgae Research Inspiring Novel Energy Resources program of the Advanced Research Projects Agency-Energy has initiated projects to develop advanced cultivation technologies that enable the cost- and energy-efficient production of macroalgal biomass. Here, this study addresses the life-cycle greenhouse gas emissions and energy return on investment for five U.S. offshore macroalgae production systems designed for deployment at the thousand-hectare scale using a detailed module developed within the GREET life-cycle analysis model for this study. The carbon intensity of macroalgae production system designs, expressed as kg of carbon dioxide equivalent per dry metric ton of algae harvested, vary widely from 49 to 220 and confirm that biomass productivity has the highest degree of sensitivity across the model parameters tested. Regardless of the system designs, the upstream and combustion emissions from fuel use are the key contributor (over 45 %) to carbon intensity, indicating that the use of low-carbon fuels (e.g., renewable diesel) could further reduce greenhouse gas emissions. Further studies need to specify the market opportunity and specific product slates for macroalgae to provide a complete picture of the environmental impacts of macroalgal feedstock.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

The Alcock–Paczynski effect from Lyman- α forest correlations: analysis validation with synthetic data

The three-dimensional distribution of the Ly α forest has been extensively used to constrain cosmology through measurements of the baryon acoustic oscillations (BAO) scale. However, more cosmological information could be extracted from the full shapes of the Ly α forest correlations through the Alcock–Paczynski (AP) effect. In this work, we prepare for a cosmological analysis of the full shape of the Ly α forest correlations by studying synthetic data of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We use a set of 100 eBOSS synthetic data sets in order to validate such an analysis. These mocks undergo the same analysis process as the real data. We perform a full-shape analysis on the mean of the correlation functions measured from the 100 eBOSS realizations, and find that our model of the Ly α correlations performs well on current data sets. We show that we are able to obtain an unbiased full-shape measurement of D M /D H (z eff ), where D M is the transverse comoving distance, D H is the Hubble distance, and z eff is the effective redshift of the measurement. We test the fit over a range of scales, and decide to use a minimum separation of r min = 25 h –1 Mpc. Here, we also study and discuss the impact of the main contaminants affecting Ly α forest correlations, and give recommendations on how to perform such analysis with real data. While the final eBOSS Ly α BAO analysis measured D M /D H (z eff = 2.33) with 4 per cent statistical precision, a full-shape fit of the same correlations could provide an $\sim 2~{{\ \rm per\ cent}}$ measurement.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

DESI DR2 Galaxy luminosity functions

We present galaxy luminosity functions (LFs) for the Dark Energy Spectroscopic Instrument (DESI) DR2 Bright Galaxy Survey (BGS) in the g, r, z, and $w1$ bands over $0.002\lt z\lt 0.6$. Our analysis uses updated k-corrections and evolutionary corrections, including new polynomial kcorrection fits derived from BGS Year 1 data that supersede earlier GAMA-based prescriptions. Exploiting the statistical power of DESI, we measure LFs to very faint magnitudes, reaching $^{0.1}M_r-5\log h\sim -10$. Independent measurements from the North and South survey regions show excellent agreement around the LF knee, but the very small statistical uncertainties reveal that simple analytic forms fail to capture the full LF shape. The bright end departs from a pure exponential decline, while the faint end exhibits complex, non-powerlaw behaviour, including a pronounced upturn at $^{0.1}M_r-5\log h\gtrsim -15$, which is stronger for red galaxies than for blue. We show that our LFs are largely complete for galaxies with surface brightness $\mu _{50}\lt 25$, and that an apparent steepening fainter than $-13$ is driven primarily by local overdensity and fragmentation of large galaxies. A systematic North–South offset at the brightest magnitudes is traced to red galaxies and may reflect shallower North photometry underestimating extended earlytype profiles, although this remains inconclusive. We therefore also provide LFs based on model Petrosian magnitudes. Redshift splitting reveals small but significant residuals, indicating limitations of a simple global evolutionary model. Using the redshift limits of J. Loveday et al. (2012), we find excellent agreement with GAMA, with substantially reduced statistical errors. These measurements provide a precise reference for studies of environmental and population-dependent LFs and for testing galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Colour gradients of low-redshift galaxies in the DESI Legacy Imaging Survey

ABSTRACT Radial colour gradients within galaxies arise from gradients of stellar age, metallicity, and dust reddening. Large samples of colour gradients from wide-area imaging surveys can complement smaller integral-field spectroscopy data sets and can be used to constrain galaxy formation models. Here, we measure colour gradients for low-redshift galaxies (z < 0.1) using photometry from the DESI Legacy Imaging Survey DR9. Our sample comprises ∼93 000 galaxies with spectroscopic redshifts and ∼574 000 galaxies with photometric redshifts. We focus on gradients across a radial range 0.5Reff to Reff, which corresponds to the inner disc of typical late-type systems at low redshift. This region has been the focus of previous statistical studies of colour gradients and has recently been explored by spectroscopic surveys such as MaNGA. We find that the colour gradients of most galaxies in our sample are negative (redder towards the centre), consistent with the literature. We investigate empirical relationships between colour gradient, average g − r and r − z colour, Mr, M⋆, and sSFR. Trends of gradient strength with Mr (M⋆) show an inflection around Mr ∼ −21 ($\log _{10} \, M_\star /\mathrm{M_\odot }\sim 10.5$). Below this mass, colour gradients become steeper with increasing M⋆, whereas colour gradients in more massive galaxies become shallower. We find that positive gradients (bluer stars at smaller radii) are typical for galaxies of $M_{\star }\sim 10^{8}\, \mathrm{M_\odot }$. We compare our results to age and metallicity gradients in two data sets derived from fits of different stellar population libraries to MaNGA spectra, but find no clear consensus explanation for the trends we observe. Both MaNGA data sets seem to imply a significant contribution from dust reddening, in particular, to explain the flatness of colour gradients along the red sequence.

Astronomy & Astrophysics↗

DOE OSTI2024

Investigation of the $\gamma $-ray Properties of the $2^+$ States in $^{14}$C

The properties of the 2$^+_1$ and 2$^+_2$ excited states in 14 C were studied in an experiment conducted at Argonne National Laboratory. A 9 Be( 6 Li,pγ) fusion-evaporation reaction and the GRETINA-ORRUBA setup were employed to populate states of 14 C and detect γ-particle coincidence events. Finally, the precise determination of the 2$^+_1$ level energy, complemented by the estimation of the γ-ray branch of the 2$^+_2$ near-threshold state, will serve as a benchmark to test the Shell Model Embedded in the Continuum calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Ultrafast Carrier–Lattice Interactions and Interlayer Modulations of Bi 2 Se 3 by X-ray Free-Electron Laser Diffraction

We report as a 3D topological insulator, bismuth selenide (Bi 2 Se 3 ) has potential applications for electrically and optically controllable magnetic and optoelectronic devices. Understanding the coupling with its topological phase requires studying the interactions of carriers with the lattice on time scales down to the sub-picosecond regime. Here we investigate the ultrafast carrier-induced lattice contractions and interlayer modulations in Bi 2 Se 3 thin films by time-resolved diffraction using an X-ray free-electron laser. The lattice contraction depends on the carrier concentration and is followed by an interlayer expansion accompanied by oscillations. Using density functional theory (DFT) and the Lifshitz model, the initial contraction can be explained by van der Waals force modulation of the confined free carrier layers. Our theory suggests that band inversion, related to a topological phase transition, is modulated by the expansion of the interlayer distance. These results provide insights into the topological phase control by light-induced structural change on ultrafast timescales.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2020

From core collapse to superluminous: the rates of massive stellar explosions from the Palomar Transient Factory

ABSTRACT We present measurements of the local core-collapse supernova (CCSN) rate using SN discoveries from the Palomar Transient Factory (PTF). We use a Monte Carlo simulation of hundreds of millions of SN light-curve realizations coupled with the detailed PTF survey detection efficiencies to forward model the SN rates in PTF. Using a sample of 86 CCSNe, including 26 stripped-envelope SNe (SESNe), we show that the overall CCSN volumetric rate is $r^\mathrm{CC}_v=9.10_{-1.27}^{+1.56}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$ at 〈z〉 = 0.028, and the SESN volumetric rate is $r^\mathrm{SE}_v=2.41_{-0.64}^{+0.81}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$. We further measure a volumetric rate for hydrogen-free superluminous SNe (SLSNe-I) using eight events at z ≤ 0.2 of $r^\mathrm{SLSN-I}_v=35_{-13}^{+25}\, \text{SNe yr}^{-1}\text{Gpc}^{-3}\, h_{70}^{3}$, which represents the most precise SLSN-I rate measurement to date. Using a simple cosmic star formation history to adjust these volumetric rate measurements to the same redshift, we measure a local ratio of SLSN-I to SESN of ${\sim}1/810^{+1500}_{-94}$, and of SLSN-I to all CCSN types of ${\sim}1/3500^{+2800}_{-720}$. However, using host galaxy stellar mass as a proxy for metallicity, we also show that this ratio is strongly metallicity dependent: in low-mass (logM* < 9.5 M⊙) galaxies, which are the only environments that host SLSN-I in our sample, we measure an SLSN-I to SESN fraction of $1/300^{+380}_{-170}$ and $1/1700^{+1800}_{-720}$ for all CCSN. We further investigate the SN rates a function of host galaxy stellar mass, and show that the specific rates of all CCSNe decrease with increasing stellar mass.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗