Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods in transport”

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 523 records · Page 29

Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics

We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.

97 MATHEMATICS AND COMPUTING↗

Sparge Sampling of Molten Salts for Online Monitoring via Laser-Induced Breakdown Spectroscopy

A method was developed to sample molten salts by sparging to generate and transport aerosols to an isolated instrument for compositional analysis by laser-induced breakdown spectroscopy (LIBS). Real-time monitoring of molten salt composition is critical to developing molten salt nuclear reactors, which offer enhanced safety and efficiency. In this article, the sparge sampling method is described and compared with sampling using a Collison nebulizer. The size distribution and transport of aerosols produced from molten eutectic NaNO 3 –KNO 3 salt were compared for multiple gas flow rates (75–1200 mL min –1 ) and transport distances (0.68–2.61 m). Both methods produced aerosols ranging from 0.5 to 5.0 μm determined using a cascade impactor. Aerosols were effectively transported without pre- or trace-heating of gas lines, but transport efficiency was reduced by the formation of agglomerates. Sparge sampling was found to use less sample and less gas than a Collison nebulizer while producing a more concentrated aerosol stream (up to 5 μg L –1 ). The effects of laser energy and delay time on the signal quality of LIBS measurements of these aerosols were also studied. High energy and short delay times were found to enhance signal and repeatability, whereas signal-to-background and signal-to-noise ratios were highest at low energy and longer delay times. The capabilities of this system for online monitoring of molten salts were demonstrated with calibrations for Sr and Li with relative standard deviations of 2.6% and 1.5% and limits of detection of 380 and 180 μg g –1 , respectively.

Aerosols↗

LANL Critical Benchmark Comparison Study and Subsequent Revision

As part of an international collaboration within the DOE Nuclear Criticality Safety Program (NCSP), LANL is involved in a comparison study to quantify differences in k-effective results from neutron transport simulations of critical benchmark experiments. The DOE NCSP Mission and Vision details the activity in which the French Institut De Radioprotection et De Sûreté Nucléaire (IRSN) leads the study with LANL and in conjunction with ORNL and LLNL to compare results of various neutron transport codes and nuclear data libraries to compute k-effective for ICSBEP benchmarks held in common by the entities. This report documents results obtained through partial completion of the overall effort with a focus on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. The feedback received through participation in the comparison collaboration has prompted an effort to review particular input files for benchmarks and revise when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, this effort prompted a new collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Noninvasive measurements of spin transport properties of an antiferromagnetic insulator

Antiferromagnetic insulators (AFIs) are of substantial interest because of their potential in the development of next-generation spintronic devices. One major effort in this emerging field is to harness AFIs for long-range spin information communication and storage. Here, we report a noninvasive method to optically access the intrinsic spin transport properties of an archetypical AFI α-Fe 2 O 3 via nitrogen-vacancy (NV) quantum spin sensors. By NV relaxometry measurements, we successfully detect the frequency-dependent dynamic fluctuations of the spin density of α-Fe 2 O 3 along the Néel order parameter, from which an intrinsic spin diffusion constant of α-Fe 2 O 3 is experimentally measured in the absence of external spin biases. Our results highlight the significant opportunity offered by NV centers in diagnosing the underlying spin transport properties in a broad range of high-frequency magnetic materials such as two-dimensional magnets, spin liquids, and magnetic Weyl semimetals, which are challenging to access by the conventional measurement techniques.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Integral Kernel Methods for Nonlinear Parabolic-Elliptic Systems

Nonlinear parabolic-elliptic systems arise in many physical, biological, and chemical phenomena such as chemotaxis, ion transport, self-gravitating particles, and Brownian vortices. Existing methods struggle with the strong coupling and high nonlinearity and nonlocality of some of these systems, especially the ill-conditioned, convection-dominated problems. To overcome numerical difficulties, current approaches rely on initial guesses, preconditioning, or iterative techniques with no convergence guarantees. They might suffer from poor scalability, large memory usage, and difficulty to parallelize. Inspired by the connection of parabolic-elliptic systems to stochastic processes, we introduce a novel meshless, monolithic, and fully explicit method that naturally encapsulates the elliptic and parabolic operators into a single step which updates each node deterministically with global information. By being fully quadrature-based, it avoids solving systems of discretized equations and does not utilize initial guesses or preconditioning, while requiring little memory and being easy to parallelize. We first derive the method in an integral kernel formulation with quadratic complexity in the number of integration nodes and then leverage kernel-independent fast multipole methods (FMM) to present a scalable algorithm with linear complexity. We provide numerical examples for the Poisson-Nernst-Planck equations in one, two, and three dimensions, together with the derivation of the integral kernel for each case. Furthermore, the examples demonstrate the fast convergence and scalability of the FMM-accelerated algorithm, as well as its suitability for convection-dominated problems, making it competitive against traditional PDE solvers.

PDE systems↗

MFiX: Fractional-Step Method Implementation

A comprehensive, multiphase computational fluid dynamics (CFD) simulation solves several coupled transport equations including continuity, momentum, species, and energy. Chemical reactions further couple these equations through heats of reaction and rates of formation of products and rates of destruction of reactants. A fractional-step method separates changes attributed to chemical reactions from transport phenomena like convection and diffusion. When the governing equations are split into the transport and reacting components, efficient and independent methodologies can be exploited to solve the different systems. Specifically, discretization of field variable transport equations results in large, sparse matrices which are loosely coupled. These systems are solved in succession using iterative techniques that take advantage of the matrix structure. In contrast, chemical reactions tightly couple field variables locally within the domain (e.g., within a single computational cell) resulting in low dimensional but dense, nonlinear systems that are better solved using direct integration techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl 2 from T m + 85 K to T m + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl 2 . We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. Furthermore, this result may contribute to new understandings of the structural relaxation of other network liquids.

36 MATERIALS SCIENCE↗

Preferential Solute Transport in Low Permeability Zones During Spontaneous Imbibition in Heterogeneous Porous Media

Multiphase flows in porous media|and the associated solute transport processes|are controlled by a combination of gravity, capillary, and viscous forces. Geologic heterogeneity influences flow and transport processes, yet gaps exist in our understanding of how heterogeneity impacts capillary-driven transport, such as during spontaneous imbibition. Here we use positron emission tomography, combined with a newly developed method for conducting spontaneous imbibition experiments, to observe solute transport into low permeability regions during imbibition. Using an experimentally parameterized 2D numerical model, we demonstrate how capillary-driven ow near the imbibition front can carry solutes into low permeability regions. This process displaces solutes from high permeability zones while the cumulative amount of solute in low permeability zones increases- opposite of what is observed under fully saturated solute transport conditions. These results highlight the complex flow and solute transport behavior that arise during multiphase displacements in heterogeneous geologic porous media.

58 GEOSCIENCES↗

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗

Evaluation of the Radioactive Material Release in the Harborview Research and Training Building and Implications for Emergency Response

On May 2, 2019, during the 137 Cs source recovery operation, a source capsule in a research irradiator containing approximately 77.1 TBq was breached. Based on a geometric reconstruction analysis of the damage to the capsule, approximately 46.3 GBq (0.04%) was impacted by the chop saw (grinder) inside a Mobile Hot Cell (MHC) on the loading dock at the University of Washington Harborview Research and Training (HRT) Building. A very small fraction of the material impacted, less than 1%) was released from the Mobile Hot Cell and then to the rest of the HRT Building. The objectives of this project were to assess the accidental release of 137 CsCl and its implications related to emergency response methods and the ramifications of 137 CsCl transport. The phenomenology of this event was also compared with past alkali halide dispersal events. The vast number of measurements and samples collected by the remediation contractors, the Department of Energy's Nuclear Emergency Support Team, and the small number of retrospective samples collected by the authors informed the analysis. The techniques included (1) autoradiography and electron microscopy of samples collected from the HRT Building and the irradiator, (2) 3D visualization of deposition on surfaces and within the ventilation system, and (3) a study of the damage to the source capsule to evaluate the Cs particle size and particle composition due to the grinding accident. Subsequently, the cesium contaminant transport through the numerous pathways in the building was reconstructed to assess the deposition on surfaces as a function of particle size. The implications for emergency response are relevant to data quality and management. A Data Quality Objective (DQO) guides data collection methods so that they have appropriate accuracy and precision for the intended application. Recommendations were made with respect to the sample collection protocols and sample archival.

61 RADIATION PROTECTION AND DOSIMETRY↗

Ab initio electron dynamics in high electric fields: Accurate prediction of velocity-field curves

Electron dynamics in external electric fields governs the behavior of solid-state electronic devices. First-principles calculations enable precise predictions of charge transport in low electric fields. However, studies of high-field electron dynamics remain elusive due to a lack of accurate and broadly applicable methods. In this work, we develop an efficient approach to solve the real-time Boltzmann transport equation with both the electric field term and ab initio electron-phonon collisions. These simulations provide field-dependent electronic distributions in the time domain, allowing us to investigate both transient and steady-state transport in electric fields ranging from low to high (> 10 kV/cm). The broad capabilities of our approach are shown by computing nonequilibrium electron occupations and velocity-field curves in Si, GaAs, and graphene, obtaining results in quantitative agreement with experiment. Our approach sheds light on microscopic details of transport in high electric fields, including the dominant scattering mechanisms and valley occupation dynamics. Our results demonstrate quantitatively accurate calculations of electron dynamics in low to high electric fields, with broad application to power and micro-electronics, optoelectronics, and sensing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Understanding fission gas bubble distribution, lanthanide transportation, and thermal conductivity degradation in neutron-irradiated α-U using machine learning

U—10Zr based metallic nuclear fuel is the leading candidate for next-generation sodium-cooled fast reactors in the United States. US research reactors have used and tested this fuel type since the 1960s and accumulated considerable experience and knowledge about the fuel performance. Most of the knowledge, however, remains empirical. The lack of mechanistic understanding of fuel performance puts a large burden on proof through experimental verification for the qualification of U—10Zr fuel for commercial use. Further, this paper proposes an image data-driven machine learning approach, coupled with domain knowledge provided by advanced post irradiation examination, to provide unprecedented quantified insights into the morphology, size, density and the connectivity of fission gas bubbles and their effects on the fission product transportation and thermal conductivity. Specifically, we developed a method to automatically detect, extract statistics, and classify ~19,000 fission gas bubbles into different categories, and quantitatively link the data to lanthanide transportation through connected bubbles and degradation of thermal conductivity along the radial temperature gradient in a neutron irradiated U—10Zr annular fuel. Results indicate the approach can be modified to study other irradiation effects, such as secondary phase redistribution and gaseous fuel swelling in other irradiated nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Charge Carrier Transport in Iron Pyrite Thin Films: Disorder-Induced Variable-Range Hopping

Here, the origin of p-type conductivity and the mechanism responsible for low carrier mobility were investigated in pyrite (FeS 2 ) thin films. Temperature-dependent resistivity measurements (10-400 K) were performed on polycrystalline and nanostructured thin films prepared by three different methods: (1) spray pyrolysis, (2) hot-injection synthesized and spin-coated nanocubes, and (3) pulsed laser deposition. The films have a high hole density (10 18 -10 19 ) cm -3 and low mobility (0.1-4 cm 2 V -1 s -1 ) regardless of the method used for their preparation. The charge transport mechanism is determined to be thermally activated conduction (TAC) at near room temperature, with Mott-type variable-range hopping (VRH) of holes via localized states occurring at lower temperatures. The density functional theory (DFT) predicts that sulfur vacancy induces localized defect states within the band gap and the charge remains localized around the defect. The data indicates that the electronic properties including hopping transport in pyrite thin films can be correlated to sulfur vacancy-related defects. The results provide insights into the electronic properties of pyrite thin films and their implications for charge transport.

36 MATERIALS SCIENCE↗

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the Radioactive Material Released in the Harborview Research and Training Building and Some Implications for Emergency Response

On 2 May 2019, during the 137 Cs source recovery operation, a source capsule in a research irradiator containing approximately 77.1 TBq was breached. Based on a geometric reconstruction analysis of the damage to the capsule, approximately 46.3 GBq (0.04%) was impacted by the chop saw (grinder) inside a mobile hot cell on the loading dock at the University of Washington Harborview Research and Training (HRT) Building. A very small fraction of the material impacted, less than 1%, was released from the mobile hot cell and then to the rest of the HRT Building. The objectives of this project were to assess the accidental release of 137 CsCl and its implications related to emergency response methods and the ramifications of 137 CsCl transport. The phenomenology of this event was also compared with past alkali halide dispersal events. The vast number of measurements and samples collected by the remediation contractors, the Department of Energy’s Nuclear Emergency Support Team, and the small number of retrospective samples collected by the authors informed the analysis. The techniques included (1) autoradiography and electron microscopy of samples collected from the HRT Building and the irradiator, (2) 3D visualization of deposition on surfaces and within the ventilation system, and (3) a study of the damage to the source capsule to evaluate the Cs particle size and particle composition due to the grinding accident. Subsequently, the cesium contaminant transport through the numerous pathways in the building was reconstructed to assess the deposition on surfaces as a function of particle size. Furthermore, the implications for emergency response are relevant to data quality and management. A Data Quality Objective guides data collection methods so that they have appropriate accuracy and precision for the intended application. Recommendations were made with respect to the sample collection protocols and archiving of samples.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Optimization of particle tracking methods for stochastic media

Random media emerge in several applications involving particle transport, encompassing e.g. photon propagation through Rayleigh-Taylor instabilities in fuel pellets for inertial confinement fusion, or neutron multiplication problems related to the assessment of re-criticality risk following severe accidents with fuel degradation. Reference calculations in such material configurations by means of Monte Carlo transport codes are particularly challenging, since high-density stochastic media might involve several hundreds of thousands of volumes and thus make particle tracking routines extremely cumbersome. In order to cope with these issues, two distinct strategies have been proposed so far: the use of neighbor maps, or the use of delta tracking. In this work we will compare these methods and illustrate their specific merits and drawbacks, as taken both alone and in combination with each other. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A simple, validated approach for design of two-dimensional periodic particle patterns via acoustophoresis

Two-dimensional patterning of microparticles enables a wide range of functional materials, including patterned energy storage electrodes, flexible electronics, and sensor arrays. Particle patterning via acoustics offers an attractive path to generate a wide variety of 2D periodic patterns that introduce tailorable hierarchical porosity, useful for controlling surface area, transport distances, and other properties. This method is most effective with micron scale particles and patterns of tens to hundreds of microns. To enable systematic exploration of the broad design space for such patterns, this work develops a model of 2D and 3D assembly of particles at high loadings and validates the obtained patterns against both experiments and more computationally intensive modeling techniques. Using this simple model, connections are mapped between input parameters (like actuation conditions, particle volume fraction, material properties) and output geometrical features (like void size and shape, pattern connectivity, and surface area) so that they can be tailored to given applications. The utility of this simple model is illustrated by predicting and then experimentally demonstrating new hierarchical patterns resulting from multiple waves of different frequencies interacting. These multiscale patterns offer the potential to lift the limits on surface area, diffusion distances, and other features.

42 ENGINEERING↗