Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective medium approximations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ruling out Color Transparency in quasi-elastic 12C(e,e′p) up to Q2=14 (GeV/c)2*

Color transparency (CT) is a fundamental phenomenon of QCD postulating that at high momentum transfer, one can preferentially measure hadrons that fluctuate to a small color neutral transverse size in the nucleus, and final state interactions within the nuclear medium are suppressed. This talk will discuss the recent quasi-elastic 12C(?,???) scattering measurement in Hall C at momentum transfer squared ?2=8,9.4,11.4, and 14.2 (GeV/c)2, the highest ever achieved to date. Nuclear transparency for this reaction was extracted by comparing the measured yield to that expected from a plane-wave impulse approximation calculation without any final state interactions. The measured transparency was observed to be independent of ?2, ruling out the quantum chromodynamics effect of color transparency at such momentum scales. These new results impose strict constraints on models of color transparency for protons.

Szumila-Vance, Holly↗
DOE OSTI2026

Computational schemes for the Magnus expansion of the in-medium similarity renormalization group

The in-medium similarity renormalization group (IMSRG) is a popular many-body method used for computations of nuclei. It solves the many-body Schrödinger equation through a continuous unitary transformation of the many-body Hamiltonian. The IMSRG transformation is typically truncated at the normal-ordered two-body level, the IMSRG(2), but recently several approaches have been developed to capture the effects of normal-ordered three-body operators, the IMSRG(3). In particular, a factorized approximation to the IMSRG(3) proposes to capture the leading effects of three-body operators at the same computational cost as the IMSRG(2) approximation. This approach often employs an approximate scheme for solving the IMSRG equations, the so-called hunter-gatherer scheme. In this work, I study the uncertainty associated with this scheme. I find that the hunter-gatherer scheme differs by up to 7MeV for ground-state energies and 0.5MeV for excitation energies from standard IMSRG(2) approaches. These differences are in some cases comparable to the expected size of IMSRG(3) corrections.

39 ≤ A ≤ 58↗
DOE OSTI2022

Metastable phase of UTe 2 formed under high pressure above 5 GPa

Uranium ditelluride (UTe 2 ) has attracted recent interest due to its unique superconducting properties, which include the potential for a topological odd-parity superconducting state. Recently, ac-calorimetry measurements under pressure indicate a change in the ground state of UTe 2 from superconducting to antiferromagnetic at 1.4 GPa. In this work, we investigate the effect of pressure on the crystal structure of UTe 2 up to 25 GPa at room temperature using x-ray diffraction. We find that UTe 2 , which at ambient conditions has an orthorhombic (Immm) structure, transforms to a body-centered tetragonal (I4/mmm) structure at 5 GPa in a quasihydrostatic neon (Ne) pressure-transmitting medium. In the absence of a pressure-transmitting medium, this transformation occurs between 5 and 8 GPa. The data were fit with a third-order Birch-Murnaghan equation of state resulting in values of B 0 = 46.0 ± 0.6 GPa , B' = 9.3 ± 0.5 (no pressure medium), and B 0 = 42.5 ± 2.0 GPa , B' = 9.3 (fixed) (neon pressure medium) for the Immm phase. For the I4/mmm phase, B 0 = 78.9 ± 0.5 GPa GPa and B' = 4.2 ± 0.1 (no pressure-transmitting medium), and B 0 = 70.0 ± 1.1 GPa and B' = 4.1 ± 0.2 (neon pressure medium). The high-pressure tetragonal phase is retained after decompression to ambient pressure, with approximately 30% remaining after 2 days. We argue that the observed phase transition into a higher-symmetry structure at P ~ 5 GPa (orthorhombic to tetragonal) is accompanied by an increase in the shortest distance between uranium atoms from 3.6 Å (orthorhombic) to 3.9 Å (tetragonal), which suggests localization of the 5ƒ electrons, albeit with a 10.7% decrease in volume.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Probing the Impact of Radio-mode Feedback on the Properties of the Cool Circumgalactic Medium

We explore the influence of radio-mode feedback on the properties of the cool circumgalactic medium (CGM). To this end, we assemble a statistical sample of approximately 30,000 radio galaxies with background quasars by combining optical spectroscopic measurements of luminous red galaxies and quasars from the year 1 data set of the Dark Energy Spectroscopic Instrument and radio sources from the LOw-Frequency ARray Two-metre Sky Survey (LoTSS) DR2 catalog and the Very Large Array Sky Survey (VLASS) quick-look catalog. Galaxies with similar optical properties but with no radio counterparts in LoTSS and VLASS are selected as the control group. We measure the cool CGM properties of radio galaxies and their control samples traced by Mg ii absorption lines, including covering fraction, rest equivalent width, and gas kinematics. Our results show no significant difference in the properties of gas around radio galaxies and their control sample, indicating that the operating radio-mode feedback of massive galaxies does not produce detectable effects on the properties of the cool CGM. Finally, we show that the CGM of radio galaxies contains a nonnegligible amount of cool gas with approximately 10 10 M ⊙ . This abundance can place a stringent constraint on the radio-mode feedback models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The impact of environmental effects on active galactic nuclei: A decline in the incidence of ionized outflows

Active galactic nuclei (AGNs) have generally been considered to be less frequent in denser environments due to a lower number of galaxy-galaxy interactions and/or the removal of their gas-rich reservoirs by the dense intergalactic medium. However, recent observational and theoretical works suggest that the effect of ram-pressure stripping acting on galaxies in dense environments might reduce the angular momentum of their gas, causing it to infall towards the super massive black hole at their centre, activating the AGN phase. In this work we explore the connection between environment and nuclear activity by evaluating the variation in the incidence of ionized outflows, a common phenomenon associated with nuclear activity, in AGNs across different environments. We select a sample of approximately 3300 optical AGNs from the Sloan Digital Sky Survey Data Release 13, which we match with a group catalogue to identify galaxies in isolation or residing in groups. We further probe their environments through the projected distance to the central galaxy of the group or cluster and the projected surface density to the fifth neighbour (δ 5 ). The presence of ionized outflows is determined through the modelling of the [O III] λ5007 emission line. We find that at lower masses (< 10 10.3 M ⊙ ), the fraction of ionized outflows is significantly lower in satellite AGNs (~7%) than in isolated (~22%) AGNs, probably due to their different AGN luminosity, L [O III] , in this stellar mass range. The fraction of outflows decreases at distances closer to the central galaxy of the group or cluster for all satellite AGNs; however, only the lower-mass ones (10 9 - 10 10.3 M ⊙ ) display a significant decline with δ 5 . Although this study does not include AGNs in the densest regions of galaxy clusters, our findings suggest that AGNs in dense environments accrete less gas than those in the field, potentially due to the removal of the gas reservoirs via stripping or starvation, leading to a negative connection between environment and AGN activity. Based on our results, we propose that the observed change in the incidence of outflows when moving towards denser regions of groups and clusters could contribute to the higher gas metallicities of cluster galaxies compared to field galaxies, especially at lower masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Low-energy enhancement in the magnetic dipole γ-ray strength functions of heavy nuclei

A low-energy enhancement (LEE), which was observed experimentally in the gamma-ray strength function (γSF) describing the decay of compound nuclei, would have profound effects on r-process nucleosynthesis if it persists in heavy neutron-rich nuclei. The LEE was shown to be a feature of the magnetic dipole (M1) strength function in configuration-interaction shell-model calculations in medium-mass nuclei. However, its existence in heavy open-shell nuclei remains an open question. Here, using a combination of many-body methods, we identify a LEE in the M1 γSFs of heavy samarium nuclei. In particular, we use the static-path plus random-phase approximation (SPA+RPA), which includes static and small-amplitude quantal fluctuations beyond the mean field. Using the SPA+RPA strength as a prior, we apply the maximum-entropy method (MEM) to obtain finite-temperature M1 γSFs from exact imaginary-time response functions calculated with the shell model Monte Carlo (SMMC) method. We find that the slope of the LEE in samarium isotopes is roughly independent of the average initial energy over a wide range below the neutron separation energy. As the neutron number increases, strength transfers to a low-energy excitation, which we interpret as the scissors mode built on top of excited states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Date Release Report for Large Surface Explosion Coupling Experiment (LSECE) Nevada National Security Site

The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.

58 GEOSCIENCES↗
DOE OSTI2025

Thermodynamics and collisionality in firehose-susceptible high- β plasmas

We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma β, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-β plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time τ exceeds the ion-Larmor frequency by a β-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is ν eff ∼ β/τ, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.

astrophysical plasmas↗
DOE OSTI2020

TeraChem Cloud: A High-Performance Computing Service for Scalable Distributed GPU-Accelerated Electronic Structure Calculations

The encapsulation and commoditization of electronic structure arise naturally as interoperability, and the use of nontraditional compute resources (e.g., new hardware accelerators, cloud computing) remains important for the computational chemistry community. Here, we present TERACHEM CLOUD, a high-performance computing service (HPCS) that offers on-demand electronic structure calculations on both traditional HPC clusters and cloud-based hardware. The framework is designed using off-the-shelf web technologies and containerization to be extremely scalable and portable. Within the HPCS model, users can quickly develop new methods and algorithms in an interactive environment on their laptop while allowing TERACHEM CLOUD to distribute ab initio calculations across all available resources. This approach greatly increases the accessibility of hardware accelerators such as graphics processing units (GPUs) and flexibility for the development of new methods as additional electronic structure packages are integrated into the framework as alternative backends. Cost-performance analysis indicates that traditional nodes are the most cost-effective long-term solution, but commercial cloud providers offer cutting-edge hardware with competitive rates for short-term large-scale calculations. We demonstrate the power of the TERACHEM CLOUD framework by carrying out several showcase calculations, including the generation of 300,000 density functional theory energy and gradient evaluations on medium-sized organic molecules and reproducing 300 fs of nonadiabatic dynamics on the B800-B850 antenna complex in LH2, with the latter demonstration using over 50 Tesla V100 GPUs in a commercial cloud environment in 8 h for approximately $1250.

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

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hydrogen and Silicon Effects on Hexagonal Close Packed Fe Alloys at High Pressures: Implications for the Composition of Earth's Inner Core

Hexagonal close-packed (hcp) structured Fe-Ni alloy is believed to be the dominant phase in the Earth's inner core. This phase is expected to contain 4%–5% light elements, such as Si and H. While the effects of individual light element candidates on the equation of state (EoS) of the hcp Fe metal have been studied, their combined effects remain largely unexplored. Here, in this study, we report the equations of state for two hcp-structured Fe-Si-H alloys, namely Fe 0.83 Si 0.17 H 0.07 and Fe 0.83 Si 0.17 H 0.46 , using synchrotron X-ray diffraction measurements up to 125 GPa at 300 K. These alloys were synthesized by cold compression of Fe-9wt%Si in either pure H 2 or Ar-H 2 mixture medium in diamond-anvil cells. The volume increase caused by a H atom in hcp Fe-Si-H alloys is approximately eight times greater than that by a Si atom. We used the improved data set to develop a composition-dependent EoS that covers a wide range of compositions. Our calculated density and bulk sound velocity of hcp Fe-Si-H alloys suggest a large trade-off between Si and H contents in fitting the seismic properties of the inner core. Combining our new EoS with geophysical and geochemical constraints, we propose 1.6–3 wt% Si and 0.15–0.6 wt% H in the Earth's inner core.

58 GEOSCIENCES↗
DOE OSTI2022

Liquid Hydrogen Heavy-Duty Vehicle Safety Review and Refueling Facility Design

Liquid hydrogen (LH 2 ) used as a fuel onboard a heavy-duty vehicle can result in increased storage capacity and faster refueling relative to compressed gas. However, there are concerns about hydrogen losses from boil-off, potential safety issues, gaps in codes and standards for cryogenic hydrogen fuel, and technical challenges with LH 2 systems for widespread transportation applications. A failure modes and effects analysis (FMEA), a safety codes and standards review, and a design review of the onboard liquid hydrogen system for a heavy-duty vehicle identified some of these potential safety issues and gaps in the codes and standards. The FMEA identified some medium and low risk failure points of the conceptual design, and the design review identified how carefully pressure relief needs to be considered for LH 2 systems. In addition, a conceptual design for a LH 2 refueling station was developed. Rough capital costs for the refueling station design were $\$1 million$ and the layout occupied approximately 13,000 ft 2 . These results can be used to inform future designs and analyses for LH 2 heavy-duty vehicles.

08 HYDROGEN↗
DOE OSTI2022

The LPM effect in sequential bremsstrahlung: 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ corrections

An important question concerning in-medium high-energy parton showers in a quark-gluon plasma or other QCD medium is whether consecutive splittings of the partons in a given shower can be treated as quantum mechanically independent, or whether the formation times for two consecutive splittings instead have significant overlap. Various previous calculations of the effect of overlapping formation times have either (i) restricted attention to a soft bremsstrahlung limit, or else (ii) used the large-N c limit (where N c = 3 is the number of quark colors). In this paper, we make a first study of the accuracy of the large-N c limit used by those calculations of overlap effects that avoid a soft bremsstrahlung approximation. Specifically, we calculate the 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ correction to previous N c = ∞ results for overlap g → gg → ggg of two consecutive gluon splittings g → gg. At order 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ , there is interesting and non-trivial color dynamics that must be accounted for during the overlap of the formation times.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling

SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.

58 GEOSCIENCES↗
DOE OSTI2025

A buffering heat exchanger/thermal energy storage system for desalination applications

A buffering heat exchanger/thermal energy storage (BHXTES) system was designed using a fast-response phase-change material/graphite foam (PCM/GF) medium for desalination and potentially other industrial applications. Unlike a store-now-and-use-later TES system, this study focused on developing a dynamically balanced, continuously used, system that optimized three functions. A lab-scale prototype was designed and fabricated with a PureTemp 151/GF storage medium. The experimental data showed excellent thermal performance repeatability indicating minimal effects of PCM expansion and contraction on the GF as well as negligible PCM redistribution effects. The experimental data also served as a validation of the numerical model of the BHXTES system, after which optimization of the BHXTES system was conducted through numerical simulations for both PureTemp 151/GF and solar salt/GF as the storage medium. The results show that, with a Therminol 55 heat transfer fluid, continuous operation of the system can be realized with a solar energy source through cycles of 8-h changing for desalination and energy storage with the supplied heat followed by 16-h discharging for desalination with the stored energy. As a result, the working fluid outlet temperatures are generally in the ranges of approximately 128–160 °C for a PureTemp 151/GF medium and 198–220 °C for a solar salt/GF medium.

Buffering↗
DOE OSTI2020

Incorporating Electronic Information into Machine Learning Potential Energy Surfaces via Approaching the Ground-State Electronic Energy as a Function of Atom-Based Electronic Populations

Machine learning (ML) approximations to density functional theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present, they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new machine learning approach called "BpopNN" for obtaining efficient approximations to DFT PESs. Conceptually, the methodology is based on approaching the true DFT energy as a function of electron populations on atoms; in practice, this is realized with available density functionals and constrained DFT (CDFT). Additionally, the new approach creates approximations to this function with neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on LinHn clusters.

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

Modeling Geologic Waste Repository Systems Below Residual Saturation

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗