Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion problems”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Confined Ni-In intermetallic alloy nanocatalyst with excellent coking resistance for methane dry reforming

Carbon dioxide and methane are two main greenhouse gases which are contributed to serious global warming. Fortunately, dry reforming of methane (DRM), a very important reaction developed decades ago, can convert these two major greenhouse gases into value-added syngas or hydrogen. The main problem retarding its industrialization is the seriously coking formation upon the nickel-based catalysts. Herein, a series of confined indium-nickel (In-Ni) intermetallic alloy nanocatalysts (In x Ni@SiO 2 ) have been prepared and displayed superior coking resistance for DRM reaction. The sample containing 0.5 wt.% of In loading (In 0.5 Ni@SiO 2 ) shows the best balance of carbon deposition resistance and DRM reactivity even after 430 h long term stability test. The boosted carbon resistance can be ascribed to the confinement of core–shell structure and to the transfer of electrons from Indium to Nickel in In-Ni intermetallic alloys due to the smaller electronegativity of In. Additionally, both the silica shell and the increase of electron cloud density on metallic Ni can weaken the ability of Ni to activate C–H bond and decrease the deep cracking process of methane. The reaction over the confined InNi intermetallic alloy nanocatalyst was conformed to the Langmuir-Hinshelwood (L-H) mechanism revealed by in situ diffuse reflectance infrared Fourier transform spectroscopy (in-situ DRIFTS). This work provides a guidance to design high performance coking resistance catalysts for methane dry reforming to efficiently utilize these two main greenhouse gases.

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

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus ( T 2 ) and associated systems on the plane ( R 2 ), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Experiments and numerical modeling of surrogate droplet combustion for a 93 octane gasoline calibration fuel

Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.

Guo, Songtao↗
DOE OSTI2023

Development of Spectroscopic Tools for Nuclear Forensics Signatures Development

Diversion of nuclear materialsfrom legitimate commercial processes for use by nefarious actors is a frightening and surreptitious problem. Smuggled nuclear material can be used to manufacture a simple radiological dispersion device (RDD) and potentially incorporated into the development of a nuclear weapon. With some frequency, smuggled nuclear material is interdicted in foreign countries as actors cross border checkpoints. Nuclear material interdicted in this way can be analyzed using a variety of tools. Plutonium-239 is commonly used as the fissile material in nuclear weapon production as an alternative to uranium-235. Plutonium in a production, refining, or finishing facility will exist in many forms including oxide precursors (PuF4, PuF3, Pu oxalate, etc.), oxide, and metal. The ability to identify plutonium in each of these forms and determine properties, characteristics, or history of these materials is crucial for determining if the processes occurring in a nuclear facility are legitimate or not. This presentation will focus on the development and application ofspectroscopic tools developed at SRNL using doubled-walled cells (DWC) to characterize Pu-bearing compounds. Advances in the use of vibrational spectroscopy techniques (Raman and IR), UV-VIS-SWIR diffuse reflectance spectroscopy, gamma emission spectroscopy, luminescence spectroscopy, laser-induced breakdown spectroscopy in conjunction with XRD and SEM in the characterization of Pu-bearing compounds will be presented with emphasis in material processing history.

Villa-Aleman, Eliel↗
DOE OSTI2021

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High-fidelity simulation-driven model development for coarse-grained computational fluid dynamics

Nuclear reactor safety analysis requires analysis of a broad range of accident scenarios and determining their consequences. For a nuclear power plant behavior, it is impossible to obtain experimental data at sufficiently large scales to support calibration and validation of various system models. In single-phase flow convective problems, high-resolution methods of Computational Fluid Dynamics (CFD) such as Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) can provide high fidelity results under conditions where physical measurements are unavailable. However, such high-resolution simulations are computationally expensive and not suitable for simulation of long transient scenarios in nuclear reactor accidents. In this work, we investigate the use of a high fidelity simulation-driven approach to model sub-grid scale (SGS) effect in Coarse Grained Computational Fluid Dynamics CG-CFD. Notably, fine-mesh simulations are used to construct physics-informed statistical surrogate of coarse-mesh SGS model. For an initial analysis, we consider a case of turbulent natural convection in a volumetrically heated fluid layer with a thermally insulated lower boundary and isothermal upper boundary. This scenario of unstable stratification is relevant to turbulent natural convection in a severe nuclear reactor accident, as well as in containment mixing and passive cooling. For this type of flow, the SGS effect is modeled by an added turbulent diffusivity in the energy equation. It is shown that a global correction for the energy equation is sufficient to achieve a significant improvement to the CG-CFD prediction of thermal mixing in the fluid layer. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Patch2Self2: Self-supervised Denoising on Coresets via Matrix Sketching

Diffusion MRI (dMRI) non-invasively maps brain white matter yet necessitates denoising due to low signal-to-noise ratios. Patch2Self (P2S) employing self-supervised techniques and regression on a Casorati matrix effectively denoises dMRI images and has become the new de-facto standard in this field. P2S however is resource intensive both in terms of running time and memory usage as it uses all voxels (n) from all-but-one held-in volumes (d-1) to learn a linear mapping Phi : \mathbb R ^ n x(d-1) \mapsto \mathbb R ^ n for denoising the held-out volume. The increasing size and dimensionality of higher resolution dMRI acquisitions can make P2S infeasible for large-scale analyses. This work exploits the redundancy imposed by P2S to alleviate its performance issues and inspect regions that influence the noise disproportionately. Specifically this study makes a three-fold contribution: (1) We present Patch2Self2 (P2S2) a method that uses matrix sketching to perform self-supervised denoising. By solving a sub-problem on a smaller sub-space so called coreset we show how P2S2 can yield a significant speedup in training time while using less memory. (2) We present a theoretical analysis of P2S2 focusing on determining the optimal sketch size through rank estimation a key step in achieving a balance between denoising accuracy and computational efficiency. (3) We show how the so-called statistical leverage scores can be used to interpret the denoising of dMRI data a process that was traditionally treated as a black-box. Experimental results on both simulated and real data affirm that P2S2 maintains denoising quality while significantly enhancing speed and memory efficiency achieved by training on a reduced data subset.

Fadnavis, Shreyas↗
DOE OSTI2015

First and second order operator splitting methods for the phase field crystal equation

In this paper, we present operator splitting methods for solving the phase field crystal equation which is a model for the microstructural evolution of two-phase systems on atomic length and diffusive time scales. A core idea of the methods is to decompose the original equation into linear and nonlinear subequations, in which the linear subequation has a closed-form solution in the Fourier space. We apply a nonlinear Newton-type iterative method to solve the nonlinear subequation at the implicit time level and thus a considerably large time step can be used. By combining these subequations, we achieve the first- and second-order accuracy in time. We present numerical experiments to show the accuracy and efficiency of the proposed methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Experimental study on magnetically insulated transmission line electrode surface evolution process under MA/cm current density

The design of high-current density magnetically insulated transmission line (MITL) is a difficult problem of current large-scale Z-pinch device. In particular, a thorough understanding of the MITL electrode surface evolution process under high current density is lacking. On the “QiangGuang-I” accelerator, the load area possesses a low inductance short-circuit structure with a diameter of 2.85 mm at the cathode, and three reflux columns with a diameter of 3 mm and uniformly distributed circumference at the anode. The length of the high density MITL area is 20 mm. A laser interferometer is used to assess and analyze the state of the MITL cathode and anode gap, and their evolution process under high current density. Experimental results indicate that evident current loss is not observed in the current density area at pulse leading edge, and peak when the surface current density reaches MA/cm. Analysis on electrode surface working conditions indicates that when the current leading edge is at 71.5% of the peak, the total evaporation of MITL cathode structure can be realized by energy deposition caused by ohmic heating. The electrode state changes, and diffusion conditions are reflected in the laser interferometer image. The MITL cathode area mainly exists in metal vapor form. The metal vapor density in the cathode central region is higher than the upper limit of laser penetration density (∼4 × 10{sup 21}/cm{sup 3}), with an expansion velocity of ∼0.96 km/s. The metal vapor density in the electrode outer area may lead to evident distortion of fringes, and its expansion velocity is faster than that in the center area (1.53 km/s).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2019

Algebraic Hybridization and Static Condensation with Application to Scalable $H$(div) Preconditioning

We propose here a unified algebraic approach for the practical application and preconditioning of static condensation and hybridization, two popular techniques in finite element discretizations. We demonstrate the use of this algebraic framework for the construction of scalable solvers for problems involving $H$(div)-spaces discretized by conforming (Raviart--Thomas) elements of arbitrary order. We illustrate through numerical experiments the relative performance of the two (in some sense dual) techniques in comparison with a state-of-the-art parallel solver, ADS [T. V. Kolev and P. S. Vassilevski, SIAM J. Sci. Comput., 34 (2012), pp. A3079--A3098], available at http://www.llnl.gov/casc/hypre and http://mfem.org. Based on these results, we recommend the use of the hybridization technique in practice, due to its clearly demonstrated superior performance with increased benefit for higher-order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Floating Island International PHASE I FINAL TECHNICAL REPORT (Advancing the Development of Floating Solar-Powered Nanobubble Aeration Systems for Use with Floating Treatment Wetlands in Natural and Man-Made Waterbodies)

More and more freshwater lakes are suffering from algae blooms, which deprive water of oxygen and lead to widespread loss of aquatic life and production of dangerous toxins. When oxygen is lacking, methane is generated in the sediments as algae is decomposed; it has been calculated that more than half of global methane emissions come from nutrient-impaired freshwater systems (Beaulieu, 2021). Nutrients, mainly from agricultural run-off, feed algae blooms and are exacerbated by warming related to climate change. Artificial aeration is frequently used to restore oxygen to a waterbody. But diffuser aeration is inefficient and expensive, needs on-site grid power, and is failing to keep pace with the demands of water as it gets warmer with climate change. In the last five years, nanobubble aeration has been introduced into freshwater applications and shows promise for rapidly increasing dissolved oxygen effectively throughout the water column. Nanobubbles deliver oxygen in bubbles that have no buoyancy, so they stay in water longer and release their oxygen more fully, compared to large bubbles that rise and burst at the surface. The prospect of this new technology answering the deficiencies of diffuser aeration drove us to initiate the current project. During our Phase I period, we have successfully tested a nanobubble aeration system that can super-saturate oxygen levels. It is operated on solar power, and both the nanobubbler and solar array are mounted on a proprietary floating island platform, that also performs biological nutrient recycling. A rudimentary system was assembled and tested on a 6.5-acre research lake at our headquarters in Montana, where it was subjected to a summer drought that reduced water levels significantly, periods of severe cold in winter, and spring conditions that included heavy rainfall and violent thunderstorms. Our twice-weekly sampling throughout the project showed that dissolved oxygen levels rose rapidly and were maintained throughout winter and spring, well into June. The nanobubbler was able to run on solar power for long periods. Methane levels were tested periodically and found to decrease as oxygen increased. The nanobubbles appeared to have no adverse impact on fish or other aquatic life exposed to them. The stability of the installation survived the weather conditions. The many problems we encountered taught us what we need to improve in the next iteration. Towards the end of our Phase I, we were awarded a supplementary state grant that enabled us to acquire a new nanobubble system for testing that runs on DC, is extremely efficient and has few moving parts. We plan to team this with low-profile solar panels, lithium batteries and a controller that provides real-time data. The ultimate goal of this project is to commercialize an affordable and scalable lake management system that will oxygenate water from the surface down to the sludge. It will enliven fisheries, prevent toxic algae blooms and – most importantly for the fate of our planet - inhibit the production and release of methane from the sediments. We anticipate that as carbon credits expand to include methane, the cost of remediating waterbodies that produce methane will be offset. We firmly believe that this is a practical technology that works and will make a big difference when fully implemented.

14 SOLAR ENERGY↗
DOE OSTI2017

Helium interactions with alumina formed by atomic layer deposition show potential for mitigating problems with excess helium in spent nuclear fuel

Helium gas accumulation from alpha decay during extended storage of spent fuel has potential to compromise the structural integrity the fuel. We report results obtained with surrogate nickel particles which suggest that alumina formed by atomic layer deposition can serve as a low volume-fraction, uniformly-distributed phase for retention of helium generated in fuel particles such as uranium oxide. Thin alumina layers may also form transport paths for helium in the fuel rod, which would otherwise be impermeable. Micron-scale nickel particles, representative of uranium oxide particles in their low helium solubility and compatibility with the alumina synthesis process, were homogeneously coated with alumina approximately 3–20 nm by particle atomic layer deposition (ALD) using a fluidized bed reactor. Particles were then loaded with helium at 800 °C in a tube furnace. Subsequent helium spectroscopy measurements showed that the alumina phase, or more likely a related nickel/alumina interface structure, retains helium at a density of at least 10 17 atoms/cm 3 . High resolution transmission electron microscopy revealed that the thermal treatment increased the alumina thickness and generated additional porosity. Results from Monte Carlo simulations on amorphous alumina predict the helium retention concentration at room temperature could reach 10 21 atoms/cm 3 at 400 MPa, a pressure predicted by others to be developed in uranium oxide without an alumina secondary phase. This concentration is sufficient to eliminate bubble formation in the nuclear fuel for long-term storage scenarios, for example. Finally, measurements by others of the diffusion coefficient in polycrystalline alumina indicate values several orders of magnitude higher than in uranium oxide, which then can also allow for helium transport out of the spent fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Extended hierarchical solvent perturbations from curved surfaces of mesoporous silica particles in a deep eutectic solvent

Here, many applications of deep eutectic solvents (DES) rely on exploitation of their unique yet complex liquid structures. Due to the ionic nature of the DES components, their diffuse structures are perturbed in the presence of a charged surface. We hypothesize that it is possible to perturb the bulk DES structure far (> 100 nm) from a curved, charged surface with mesoscopic dimensions. We performed in situ, synchrotron-based ultra-small angle X-ray scattering (USAXS) experiments to study the solvent distribution near the surface of charged mesoporous silica particles (MPS) (≈ 0.5 µm in diameter) suspended in both water and a common type of DES (1:2 choline Cl-:ethylene glycol). A careful USAXS analysis reveals that the perturbation of electron density distribution within the DES extends ≈ 1 μm beyond the particle surface, and that this perturbation can be manipulated by the addition of salt ions (AgCl). The concentration of the pore-filling fluid is greatly reduced in the DES. Notably, we extracted the real-space structures of these fluctuations from the USAXS data using a simulated annealing approach that does not require a priori knowledge about the scattering form factor, and can be generalized to a wide range of complex small-angle scattering problems.

36 MATERIALS SCIENCE↗
DOE OSTI2017

grim: A Flexible, Conservative Scheme for Relativistic Fluid Theories

Hot, diffuse, relativistic plasmas such as sub-Eddington black-hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects such as heat conduction and viscosity can be important in a collisionless plasma and will potentially alter the dynamics and radiative properties of the flow from that in ideal fluid models; we refer to models that include these processes as Extended GRMHD. Here we describe a new conservative code, grim, that enables all of the above and additional physics to be efficiently incorporated. grim combines time evolution and primitive variable inversion needed for conservative schemes into a single step using an algorithm that only requires the residuals of the governing equations as inputs. This algorithm enables the code to be physics agnostic as well as flexibility regarding time-stepping schemes. grim runs on CPUs, as well as on GPUs, using the same code. We formulate a performance model and use it to show that our implementation runs optimally on both architectures. grim correctly captures classical GRMHD test problems as well as a new suite of linear and nonlinear test problems with anisotropic conduction and viscosity in special and general relativity. As tests and example applications, we resolve the shock substructure due to the presence of dissipation, and report on relativistic versions of the magneto-thermal instability and heat flux driven buoyancy instability, which arise due to anisotropic heat conduction, and of the firehose instability, which occurs due to anisotropic pressure (i.e., viscosity). Finally, we show an example integration of an accretion flow around a Kerr black hole, using Extended GRMHD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Optimization of electrical treatment strategy for surface roughness reduction in conducting thin films

The surface roughness of deposited conducting thin films is responsible for various materials reliability problems in nanoelectronics and nanofabrication technologies. In this work, we report a modeling and simulation study that aims at optimizing the electrical surface treatment of deposited conducting thin films as a physical processing strategy for their surface roughness reduction. Our research is based on a continuum model of film surface morphological evolution that accounts for the residual stress in the deposited conducting thin film, the film’s wetting of the substrate layer that it is deposited on, film texture and surface diffusional anisotropy, and surface electromigration. Through systematic linear stability analysis and dynamical simulation protocols, we examine in detail the effects of film surface crystallographic orientation and applied electric field direction toward minimizing the electric field strength required for film surface smoothening. We believe that the critical electric field strength requirement for surface roughness reduction on {110}, {100}, and {111} surfaces of face-centered cubic crystalline conducting thin films exhibits a very strong dependence on the applied electric field direction, expressed as the electric field misalignment with respect to the principal residual stress directions in the film and the fast surface diffusion directions. Due to these findings, we optimize the electrical treatment strategy for surface roughness reduction of conducting thin films with respect to all relevant processing and material parameters.

42 ENGINEERING↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗