Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion models”

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,765 records · Page 98

An end-to-end deep learning method for solving nonlocal Allen–Cahn and Cahn–Hilliard phase-field models

Here, we propose an efficient end-to-end deep learning method for solving nonlocal Allen–Cahn (AC) and Cahn–Hilliard (CH) phase-field models. One motivation for this effort emanates from the fact that discretized partial differential equation-based AC or CH phase-field models result in diffuse interfaces between phases, with the only recourse for remediation is to severely refine the spatial grids in the vicinity of the true moving sharp interface whose width is determined by a grid-independent parameter that is substantially larger than the local grid size. In this work, we introduce non-mass conserving nonlocal AC or CH phase-field models with regular, logarithmic, or obstacle double-well potentials. Because of non-locality, some of these models feature totally sharp interfaces separating phases. The discretization of such models can lead to a transition between phases whose width is only a single grid cell wide. Another motivation is to use deep learning approaches to ameliorate the otherwise high cost of solving discretized nonlocal phase-field models. To this end, loss functions of the customized neural networks are defined using the residual of the fully discrete approximations of the AC or CH models, which results from applying a Fourier collocation method and a temporal semi-implicit approximation. To address the long-range interactions in the models, we tailor the architecture of the neural network by incorporating a nonlocal kernel as an input channel to the neural network model. We then provide the results of extensive computational experiments to illustrate the accuracy, predictive capabilities, and cost reductions of the proposed method.

42 ENGINEERING↗

The global distribution of thermospheric odd nitrogen for solstice conditions during solar cycle minimum

A two-dimensional model of odd nitrogen in the thermosphere and upper mesosphere is described. The global distributions of nitric oxide and atomic nitrogen are calculated for the solstice period for quiet and moderate magnetic activity during the solar minimum period. The effect of thermospheric transport by winds is investigated along with the importance of particle-induced ionization in the auroral zones. The results are compared with rocket and satellite measurements, and the sensitivity of the model to eddy diffusion and neutral winds is investigated. Downward fluxes of NO into the mesosphere are given, and their importance for stratospheric ozone is discussed. The results show that the summer-to-winter pole meridional circulation transports both NO and N(S-4) across the solar terminator into the polar night region where there is a downward vertical transport toward the mesosphere. The model shows that odd nitrogen densities at high winter latitudes are entirely controlled by particle precipitation and transport processes.

Gerard, J.-C.↗

Positron annihilation radiation from the Galactic center region

Observations show that there are two components of positron annihilation radiation from the region of the Galactic center: a variable component resulting from one or just a few compact sources at or near the Galactic center and a steady, diffuse component resulting from positron annihilation in the Galactic disk. The diffuse component is modeled using the observed longitude distributions of 70-150 MeV gamma rays, CO, and hot plasma revealed by Fe line emission. Recent results on positron annihilation in the interstellar medium are reviewed and the implications of the annihilation processes on the fraction of positrons annihilating via positronium and on the shape of the 511 keV annihilation line are discussed. The sources of diffuse Galactic positrons are also reviewed and the nature of the compact source of annihilation radiation near the Galactic center is discussed.

Ramaty, R.↗

Experimental and modeling study of calcium carbonate precipitation and its effects on the degradation of oil well cement during carbonated brine exposure

Decalcification of cement in solutions of carbonated brine is important to a host of engineering applications, especially in subsurface service environments where cementitious materials are frequently utilized as engineered barriers for wellbore seals, as well as shaft and drift seals and waste forms for nuclear waste disposal. Analysis of leaching simulations and experiments shows that, depending on solution compositions (dissolved CO{sub 2} concentration, pH, Ca ion concentration), calcite precipitation occurring during leaching of cement in contact with carbonated brine can have a significant impact on cement reactivity, in some instances leading to complete arrest of reactivity via calcium carbonate ''pore-clogging''. We present modeling and experimental results that examine the range of solution conditions that can lead to pore-clogging. Analysis of the results shows that distinct regimes of leaching behavior, based on pH and pCO2, can be used to form a framework to better understand the occurrence of pore-clogging.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An atomistically informed mesoscale model for growth and coarsening during discharge in lithium-oxygen batteries

An atomistically informed mesoscale model is developed for the deposition of a discharge product in a Li-O{sub 2} battery. This mescocale model includes particle growth and coarsening as well as a simplified nucleation model. The model involves LiO{sub 2} formation through reaction of O{sub 2}{sup −} and Li{sup +} in the electrolyte, which deposits on the cathode surface when the LiO{sub 2} concentration reaches supersaturation in the electrolyte. A reaction-diffusion (rate-equation) model is used to describe the processes occurring in the electrolyte and a phase-field model is used to capture microstructural evolution. This model predicts that coarsening, in which large particles grow and small ones disappear, has a substantial effect on the size distribution of the LiO{sub 2} particles during the discharge process. The size evolution during discharge is the result of the interplay between this coarsening process and particle growth. The growth through continued deposition of LiO{sub 2} has the effect of causing large particles to grow ever faster while delaying the dissolution of small particles. The predicted size evolution is consistent with experimental results for a previously reported cathode material based on activated carbon during discharge and when it is at rest, although kinetic factors need to be included. The approach described in this paper synergistically combines models on different length scales with experimental observations and should have applications in studying other related discharge processes, such as Li{sub 2}O{sub 2} deposition, in Li-O{sub 2} batteries and nucleation and growth in Li-S batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cation–Ligand Interactions Dictate Salt Partitioning and Diffusivity in Ligand-Functionalized Polymer Membranes

Membranes are an attractive alternative to current thermal separations due to their scalability and energy efficiency in desalinating water. Unfortunately, many of the conventional membrane materials available today are unable to differentiate between ionic solutes, especially alkali cations, compromising their use in ion–ion separations. Inspired by the ion-specific interactions exhibited by biological ion channels, recent research efforts have focused on synthesizing and characterizing new polymeric materials that incorporate ligands into polymer networks to bias solubility and/or diffusivity of one cationic species over another. Despite these efforts, little is known about the influence of incorporating ligands into polymer membranes on solubility and diffusivity of the complexing species. In this study, we first build a qualitative model of salt partitioning, diffusivity, and permeability in generic cation-complexing ligand-functionalized polymer membranes. Next, to validate our model and hypotheses, we perform atomistic molecular dynamics simulations of a 12-crown-4-functionalized membrane in the presence of alkali halide salts at low concentration. Generally, cation complexation enhances cation solubility but decreases diffusivity. Interestingly, the reduction in diffusivity is predicted to be larger than the enhancement in solubility for materials which operate by the mechanisms proposed in our physical picture, ultimately resulting in a reduction in the permeability of the selectively complexing ion.

36 MATERIALS SCIENCE↗

Hydrogen uptake in graphite matrix at high temperature

Tritium management is a critical challenge for the next generation of nuclear reactors, such as Fluoride Salt Cooled High Temperature Reactors (FHRs) and High Temperature Gas-cooled Reactors (HTGRs), due to the higher production rate (up to 10,000 times) than conventional Light Water Reactors (LWRs). Graphitic materials employed as moderator, reflector, and fuel pebbles offer a potential pathway for tritium recovery by serving as a sink for tritium. Prediction of uptake capacity under reactor relevant conditions remains a challenge due to a lack of low partial pressure data and significant inter-grade variability of graphite. This study addresses these gaps by providing a comprehensive characterization of hydrogen (as a tritium surrogate) uptake and release behavior in the A3-3 graphite matrix (GM) used in fuel pebbles. Uptake measurements are performed at reactor relevant temperatures of 600- 800 °C, 1-200 Torr hydrogen pressure, and 15-120 min equilibration time, followed by thermal desorption spectroscopy up to 1100 °C. Uptake experiments at different equilibration times demonstrate the role of kinetics in hydrogen uptake, which can be modeled as a diffusion-with-trapping process. In the thermodynamics limit, the Sips adsorption model is shown to capture the uptake in A3-3 GM well. Our campaign provides a set of new results for hydrogen uptake in A3-3, including limiting uptake capacity at 600 °C, apparent diffusion coefficient at 600 °C, and the first estimates of the FHR/HTGR relevant (600 °C, 20 Pa partial pressure) equilibrium uptake capacity and time to saturation. Desorption data highlights a new site for hydrogen uptake, not observed in nuclear graphite, which we attribute to the non-graphitized binder. Using the Kissinger method, we estimate activation energy for release from the desorption peaks, confirming the activation energy for release from the basal planes and providing the first estimate for the activation energy of release from the binder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling the Thermostructural Stability of Melt-infiltrated Sic/sic Composites

SiC/SiC composites developed by NASA with Sylramic-iBN fibers and melt-infiltrated (MI) SiC-Si matrices have demonstrated 1000-hour rupture life in air at 100 MPa and 1315OC. Recently it has been determined that a major factor controlling the long-term rupture life of these composites is not environment or stress, but an intrinsic microstructural and strength instability caused by a thermally-induced silicon attack of the Sic fibers. The objective of this paper is to present a simple diffusion-based analytical model which predicts well the observed effects of stress-free thermal exposure on the residual tensile strength of Sylramic-iBN/SiC-Si composites. The practical implications of the model for SiC/SiC composites with MI matrices are discussed.

DiCarlo, James A.↗

Internal protein motions in molecular-dynamics simulations of Bragg and diffuse X-ray scattering

Molecular-dynamics (MD) simulations of Bragg and diffuse X-ray scattering provide a means of obtaining experimentally validated models of protein conformational ensembles. This paper shows that compared with a single periodic unit-cell model, the accuracy of simulating diffuse scattering is increased when the crystal is modeled as a periodic supercell consisting of a 2 × 2 × 2 layout of eight unit cells. The MD simulations capture the general dependence of correlations on the separation of atoms. There is substantial agreement between the simulated Bragg reflections and the crystal structure; there are local deviations, however, indicating both the limitation of using a single structure to model disordered regions of the protein and local deviations of the average structure away from the crystal structure. Although it was anticipated that a simulation of longer duration might be required to achieve maximal agreement of the diffuse scattering calculation with the data using the supercell model, only a microsecond is required, the same as for the unit cell. Rigid protein motions only account for a minority fraction of the variation in atom positions from the simulation. The results indicate that protein crystal dynamics may be dominated by internal motions rather than packing interactions, and that MD simulations can be combined with Bragg and diffuse X-ray scattering to model the protein conformational ensemble.

59 BASIC BIOLOGICAL SCIENCES↗

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of a PBPK model of thiocyanate in rats with an extrapolation to humans: A computational study to quantify the mechanism of action of thiocyanate kinetics in thyroid

Thyroid homeostasis can be disturbed due to thiocyanate exposure from the diet or tobacco smoke. Thiocyanate inhibits both thyroidal uptake of iodide, via the sodium-iodide symporter (NIS), and thyroid hormone (TH) synthesis in the thyroid, via thyroid peroxidase (TPO), but the mode of action of thiocyanate is poorly quantified in the literature. The characterization of the link between intra-thyroidal thiocyanate concentrations and dose of exposure is crucial for assessing the risk of thyroid perturbations due to thiocyanate exposure. We developed a PBPK model for thiocyanate that describes its kinetics in the whole-body up to daily doses of 0.15 mmol/kg, with a mechanistic description of the thyroidal kinetics including NIS, passive diffusion, and TPO. The model was calibrated in a Bayesian framework using published studies in rats. Goodness-of-fit was satisfactory, especially for intra-thyroidal thiocyanate concentrations. Thiocyanate kinetic processes were quantified in vivo, including the metabolic clearance by TPO. The passive diffusion rate was found to be greater than NIS-mediated uptake rate. The model captured the dose-dependent kinetics of thiocyanate after acute and chronic exposures. Model behavior was evaluated using a Morris screening test. The distribution of thiocyanate into the thyroid was found to be determined primarily by the partition coefficient, followed by NIS and passive diffusion; the impact of the latter two mechanisms appears to increase at very low doses. Extrapolation to humans resulted in good predictions of thiocyanate kinetics during chronic exposure. The developed PBPK model can be used in risk assessment to quantify dose-response effects of thiocyanate on TH. - Highlights: • A PBPK model of thiocyanate (SCN{sup −}) was calibrated in rats in a Bayesian framework. • The intra-thyroidal kinetics of thiocyanate including NIS and TPO was modeled. • Passive diffusion rate for SCN{sup −} seemed to be greater than the NIS-mediated uptake. • The dose-dependent kinetics of SCN{sup −} was captured after an acute and chronic exposure. • The PBPK model of thiocyanate was successfully extrapolated to humans.

54 ENVIRONMENTAL SCIENCES↗

Diagnosing Model Errors in Simulation of Solar Radiation on Inclined Surfaces

Transposition models have been widely used in the solar energy industry to simulate solar radiation on inclined PV panels. Following numerous studies comparing the performance of transposition models, this paper aims to understand the quantitative uncertainty in the state-of-the-art transposition models and the sources leading to the uncertainty. Our results show significant differences between two highly used isotropic transposition models with one substantially underestimating the diffuse plane-of-array (POA) irradiances when diffuse radiation is perfectly isotropic. In the empirical transposition models, the selection of empirical coefficients and land surface albedo can both result in uncertainty in the output. This study can be used as a guide for future development of physics-based transposition models.

14 SOLAR ENERGY↗

The Fermi-LAT gamma-ray excess at the Galactic Center in the singlet-doublet fermion dark matter model

The singlet-doublet fermion dark matter model (SDFDM) provides a good DM candidate as well as the possibility of generating neutrino masses radiatively. The search and identification of DM requires the combined effort of both indirect and direct DM detection experiments in addition to the LHC. Remarkably, an excess of GeV gamma rays from the Galactic Center (GCE) has been measured with the Fermi Large Area Telescope (LAT) which appears to be robust with respect to changes in the diffuse galactic background modeling. Although several astrophysical explanations have been proposed, DM remains a simple and well motivated alternative. In this work, we examine the sensitivities of dark matter searches in the SDFDM scenario using Fermi-LAT, CTA, IceCube/DeepCore, LUX, PICO and LHC with an emphasis on exploring the regions of the parameter space that can account for the GCE. We find that DM particles present in this model with masses close to ∼ 99 GeV and ∼ (173–190) GeV annihilating predominantly into the W{sup +}W{sup −} channel and t t-bar channel respectively, provide an acceptable fit to the GCE while being consistent with different current experimental bounds. We also find that much of the obtained parameter space can be ruled out by future direct search experiments like LZ and XENON-1T, in case of null results by these detectors. Interestingly, we show that the most recent data by LUX is starting to probe the best fit region in the SDFDM model.

79 ASTRONOMY AND ASTROPHYSICS↗

Integrability breaking in the Rule 54 cellular automaton

Cellular automata have recently attracted a lot of attention as testbeds to explore the emergence of many-body quantum chaos and hydrodynamics. We consider the Rule 54 model, one of the simplest interacting integrable models featuring two species of quasiparticles (solitons), in the presence of an integrability-breaking perturbation that allows solitons to backscatter. Here, we study the onset of thermalization and diffusive hydrodynamics in this model, compute perturbatively the diffusion constant of tracer particles, and comment on its relation to transport coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Large-cell Monte Carlo renormalization of irreversible growth processes

Monte Carlo sampling is applied to a recently formulated direct-cell renormalization method for irreversible, disorderly growth processes. Large-cell Monte Carlo renormalization is carried out for various nonequilibrium problems based on the formulation dealing with relative probabilities. Specifically, the method is demonstrated by application to the 'true' self-avoiding walk and the Eden model of growing animals for d = 2, 3, and 4 and to the invasion percolation problem for d = 2 and 3. The results are asymptotically in agreement with expectations; however, unexpected complications arise, suggesting the possibility of crossovers, and in any case, demonstrating the danger of using small cells alone, because of the very slow convergence as the cell size b is extrapolated to infinity. The difficulty of applying the present method to the diffusion-limited-aggregation model, is commented on.

Nakanishi, H.↗

Machine Learning Predictions of Simulated Self-Diffusion Coefficients for Bulk and Confined Pure Liquids

Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. Here in this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optical Fier Based System for Multiple Thermophysical Properties for Glove Box, Hot Cell and In-Pile Application

Thermal diffusivity of materials is of interest in nuclear applications at temperatures in excess of 2000°C. Commercial laser flash apparatus (LFA) that heats samples with a furnace typically do not reach these elevated temperatures nor are they easily adapted to a glove-box or hot cell environment. In this research, we performed work on an experimental technique using single laser surface heating, i.e. heating the disk sample only at its front surface with the continuous wave (CW) laser, to allow measurement of thermal diffusivity at very high temperatures within a small chamber. Thermal diffusivity is measured using a separate pulsed laser on the front side and IR detector on the rear side. The new way of heating provides easy operation in comparison to other heating methods. The measurement of sample reference temperature is needed for the measured thermal diffusivity. A theoretical model was developed to describe transient heat transfer across the sample due to the laser pulse, starting from the steady state temperature of the sample heated by the CW laser. The experimental setup was established with a 500W CW laser and maximum 50 Joule pulse laser irradiated at the front surface of the sample. The induced temperature rise at the rear surface, along with the steady-state temperature at the front surface, was recorded for the determination of thermal diffusivity and the sample temperature. Three samples were tested in vacuum over a wide temperature range of 500°C to 2100°C, including graphite, Inconel 600 and tungsten. The latter two samples were coated with sprayed graphite on their front surfaces in order to achieve surface absorption/emission needs, i.e. high absorptivity of the front surface against relatively low emissivity of the rear surface. Thermal diffusivity of graphite determined by our system are within a 5% difference of the commercial LFA data at temperatures below 1300°C and agree well with its trend at higher temperatures. Good agreement would also exist for Inconel 600 and tungsten. Despite large uncertainty of measuringthe sample temperature, the uncertainties of thermal diffusivity are less than 6% for all samples at elevated temperatures. The results indicate that single laser surface heating could be convenient and practical for the application of the LFA measurements without extra uncertainty, as temperature dependence of thermal diffusivity is usually negligible in the sample. Moreover, it is concluded that unequal surface treatment, i.e., high absorption on the front side and low emission on the rear side, greatly improves the measurement in serval aspects: less power requirement of the CW laser, less uncertainty of measured thermal diffusivity, and more uniform temperature distribution in the sample. The result of this research can be used as a general guideline for the design of this type of measurement system for nuclear applications. It can also be used directly to design and build a system similar to the one implemented in this project.

36 MATERIALS SCIENCE↗

On the Radio-emitting Particles of the Crab Nebula: Stochastic Acceleration Model

The broadband emission of pulsar wind nebulae (PWNe) is well described by non-thermal emissions from accelerated electrons and positrons. However, the standard shock acceleration model of PWNe does not account for the hard spectrum in radio wavelengths. The origin of the radio-emitting particles is also important to determine the pair production efficiency in the pulsar magnetosphere. Here, we propose a possible resolution for the particle energy distribution in PWNe; the radio-emitting particles are not accelerated at the pulsar wind termination shock but are stochastically accelerated by turbulence inside PWNe. We upgrade our past one-zone spectral evolution model to include the energy diffusion, i.e., the stochastic acceleration, and apply the model to the Crab Nebula. A fairly simple form of the energy diffusion coefficient is assumed for this demonstrative study. For a particle injection to the stochastic acceleration process, we consider the continuous injection from the supernova ejecta or the impulsive injection associated with supernova explosion. The observed broadband spectrum and the decay of the radio flux are reproduced by tuning the amount of the particle injected to the stochastic acceleration process. The acceleration timescale and the duration of the acceleration are required to be a few decades and a few hundred years, respectively. Our results imply that some unveiled mechanisms, such as back reaction to the turbulence, are required to make the energies of stochastically and shock-accelerated particles comparable.

79 ASTRONOMY AND ASTROPHYSICS↗