Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion model”

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

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

At least 541 records · Page 30

Performance characterization of nanofiltration, reverse osmosis, and ion exchange technologies for acetic acid separation

A major obstacle to widespread implementation of bio-based fuels and chemicals is the lack of efficient and costeffective separation methods. To purify acetic acid produced by biochemical conversion of biomass via anaerobic digestion, this work employs two commonly used separation technologies: (1) ion-exchange (IX) resin and (2) pressure-driven membranes. This study tested five commercially available strong- and weak-base anion-exchange resins and five commercially available nanofiltration (NF) and reverse osmosis (RO) membranes. The pH of the feed solution significantly affected adsorption capacity. At pH 6.3, a strong-base IX resin (IRN-78) performed best (95.1% acetate removal). With strong-base IX resins, the Langmuir isotherm model fit well, whereas for weak-base IX resins, the Freundlich isotherm provided a better fit. A pseudo-second-order kinetic model fit well for both IRN-78 and IRA-67. Regarding membrane separation, RO (BW30XFR membrane) achieved the highest rejection (98.6% acetate rejection), whereas an NF membrane (NF*) achieved the best combination of permeate flux (105 L/(h∙m 2 )) and rejection (83.1% acetate rejection). For membrane performance, the experimental data were fit using the solution diffusion model. Increased pH in the feed solution lowered permeate flux but increased acetic acid rejection. Finally, then the acetic acid concentration in the feed solution increased, both permeate flux and acetic acid rejection decreased for membrane NF*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accumulation of Antifreeze Proteins on Ice Is Determined by Adsorption

Antifreeze proteins (AFPs) enable the survival of organisms in freezing environments by binding to ice crystals and suppressing their growth. The rate at which AFPs accumulate on ice is determined by an interplay between their diffusion from the aqueous solution and their adsorption to the ice crystal. To interrogate the relative importance of these two processes, we combined fluorescence experiments with reaction–diffusion modeling. We found that as AFPs accumulate on ice, their concentration in the aqueous solution does not develop a gradient but remains equal to its bulk concentration throughout our experiments. Furthermore, these findings led us to conclude that AFP accumulation on ice crystals, which are smaller than 100 μm in radius, is not limited by the diffusion of AFPs, but by the kinetics of AFP adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of morphogen gradient: The role of dimension and discreteness

The fundamental processes of biological development are governed by multiple signaling molecules that create non-uniform concentration profiles known as morphogen gradients. It is widely believed that the establishment of morphogen gradients is a result of complex processes that involve diffusion and degradation of locally produced signaling molecules. We developed a multi-dimensional discrete-state stochastic approach for investigating the corresponding reaction-diffusion models. It provided a full analytical description for stationary profiles and for important dynamic properties such as local accumulation times, variances, and mean first-passage times. The role of discreteness in developing of morphogen gradients is analyzed by comparing with available continuum descriptions. It is found that the continuum models prediction about multiple time scales near the source region in two-dimensional and three-dimensional systems is not supported in our analysis. Using ideas that view the degradation process as an effective potential, the effect of dimensionality on establishment of morphogen gradients is also discussed. In addition, we investigated how these reaction-diffusion processes are modified with changing the size of the source region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗

Phase space effects on fast ion distribution function modeling in tokamaks

Integrated simulations of tokamak discharges typically rely on classical physics to model energetic particle (EP) dynamics. However, there are numerous cases in which energetic particles can suffer additional transport that is not classical in nature. Examples include transport by applied 3D magnetic perturbations and, more notably, by plasma instabilities. Focusing on the effects of instabilities, ad-hoc models can empirically reproduce increased transport, but the choice of transport coefficients is usually somehow arbitrary. New approaches based on physics-based reduced models are being developed to address those issues in a simplified way, while retaining a more correct treatment of resonant wave-particle interactions. The kick model implemented in the tokamak transport code TRANSP is an example of such reduced models. It includes modifications of the EP distribution by instabilities in real and velocity space, retaining correlations between transport in energy and space typical of resonant EP transport. The relevance of EP phase space modifications by instabilities is first discussed in terms of predicted fast ion distribution. Results are compared with those from a simple, ad-hoc diffusive model. It is then shown that the phase-space resolved model can also provide additional insight into important issues such as internal consistency of the simulations and mode stability through the analysis of the power exchanged between energetic particles and the instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling Tritium Production and Release at High-Energy Proton Accelerators

Tritium is a well-known byproduct of particle accelerator operations. To keep levels of tritium below regulatory limits, tritium production is actively monitored and managed at Fermilab. We plan to study tritium production in the targets, beamline components, and shielding elements of the Fermilab facilities such as NuMI, BNB, and MI-65. To facilitate the analysis, we construct a simple model and use three Monte Carlo radiation codes, FLUKA, MARS, and PHITS, to estimate the amount of tritium produced in these facilities. The analysis could also serve as an intercomparison between these code results related to tritium production. To assess the actual amounts of tritium that would be released from various materials, we employ a semi-empirical diffusion model. The results of this analysis are compared to experimental data whenever possible. This approach also helps to optimize proposed target materials with respect to the tritium production and release.

Georgobiani, Dali [Fermilab]↗

Magneto-transport and microstructure of Co{sub 2}Fe(Ga{sub 0.5}Ge{sub 0.5})/Cu lateral spin valves prepared by top-down microfabrication process

Heusler alloy-based lateral spin valves with ohmic contacts are prepared for the Co{sub 2}Fe(Ga{sub 0.5}Ge{sub 0.5})/Cu system by means of the top-down microfabrication process. The magneto-transport and microstructure are characterized to investigate the influence of the microfabrication route on the spin dependent transport of lateral spin valves systematically. A large non-local spin signal (△R{sub S}) of 17.3 mΩ is observed at room temperature, which is attributed to the highly spin-polarized Co{sub 2}Fe(Ga{sub 0.5}Ge{sub 0.5}) ferromagnet and the clean Co{sub 2}Fe(Ga{sub 0.5}Ge{sub 0.5})/Cu interfaces confirmed by transmission electron microscopy. Based on the general expression of one-dimensional spin diffusion model, we discuss the importance of interfacial spin polarization in Heusler alloy-based lateral spin valves.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Current-based detection of nonlocal spin transport in graphene for spin-based logic applications

Graphene has been proposed for novel spintronic devices due to its robust and efficient spin transport properties at room temperature. Some of the most promising proposals require current-based readout for integration purposes, but the current-based detection of spin accumulation has not yet been developed. In this work, we demonstrate current-based detection of spin transport in graphene using a modified nonlocal geometry. By adding a variable shunt resistor in parallel to the nonlocal voltmeter, we are able to systematically cross over from the conventional voltage-based detection to current-based detection. As the shunt resistor is reduced, the output current from the spin accumulation increases as the shunt resistance drops below a characteristic value R*. We analyze this behavior using a one-dimensional drift-diffusion model, which accounts well for the observed behavior. These results provide the experimental and theoretical foundation for current-based detection of nonlocal spin transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Polarized imaging of visible synchrotron emission from runaway electron plateaus in DIII-D

A new approach to estimate the pitch angle of highly energetic postdisruption runaway electrons is presented by imaging the polarized runaway electron synchrotron emission. The vertical to horizontal polarization ratio is found to range between 3 and 14 in image-average studies looking at repeated shots as well as the polarization in a single discharge. Numerical studies with guiding center and full orbit simulations performed with the SOFT and KORC codes are given and compared to the experiment. The pitch angle evolution analyzed using a 1D impurity diffusion model in combination with kinetic simulations is also presented. Furthermore, these studies find agreement in the magnitude of the pitch angle, but not in the temporal rate of change of pitch angle. The measured pitch angle evolution timescale of order 50 ms falls between a modeled timescale of 20 ms in an equilibrium model and >>100 ms in a nonequilibrium model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimal Experimental Design for Inverse Problems in the Presence of Observation Correlations

Optimal experimental design (OED) is the general formalism of sensor placement and decisions on the data collection strategy for engineered or natural experiments. This problem is prevalent in many critical fields such as battery design, numerical weather prediction, geosciences, environmental and urban studies. State-of-the-art computational methods for experimental design do not accommodate correlation structure in observational errors produced by many expensive-to-operate devices such as X-ray machines, radars, and satellites. Discarding evident data correlations leads to biased results, higher expenses, and waste of valuable resources. We present a general formulation of the OED formalism for model-constrained large-scale Bayesian linear inverse problems, where measurement errors are generally correlated. The proposed approach utilizes the Hadamard product of matrices to formulate the weighted-likelihood, and is valid for both finite as well as infinite-dimensional Bayesian inverse problems. Extensive numerical experiments are carried out for empirical verification of the proposed approach using an advection-diffusion model, where the objective is to optimally place a small set of sensors, under a limited budget, to predict the concentration of a contaminant in a closed and bounded domain.

Correlated observations↗

Spatial separation in a thermal mixture of ultracold {sup 174}Yb and {sup 87}Rb atoms

We report on the observation of unusually strong interactions in a thermal mixture of ultracold atoms which cause a significant modification of the spatial distribution. A mixture of {sup 87}Rb and {sup 174}Yb with a temperature of a few {mu}K is prepared in a hybrid trap consisting of a bichromatic optical potential superimposed on a magnetic trap. For suitable trap parameters and temperatures, a spatial separation of the two species is observed. We infer that the separation is driven by a large interaction strength between {sup 174}Yb and {sup 87}Rb accompanied by a large three-body recombination rate. Based on this assumption we have developed a diffusion model which reproduces our observations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

In-Situ Spatial Mapping of Hydrogen in Yttrium Hydrides at LANSCE (FY23 Version, Rev. 1)

This report summarizes the development of neutron imaging capabilities and experimental activities performed at the Los Alamos Neutron Science Center (LANSCE) with the main goal of measuring temperature-driven hydrogen diffusion within bulk-yttrium hydride (YH x ) materials. Yttrium hydride is the leading candidate to serve as a solid neutron moderator in microreactor cores, owing to its high density of hydrogen atoms as well as its superior thermal stability compared to all other metal hydrides. The experimental results and technique developments reported herein support the U.S. Department of Energy Office of Nuclear Energy’s (DOE-NE) Microreactor Program under Technology Maturation. In particular, it addresses the critical need to experimentally validate and verify hydrogen-diffusion models of metal hydrides used in high-temperature microreactor designs by means of high-spatial-resolution neutron imaging. These capabilities were designed to apply large temperature gradients across centimeter-sized YH x pellets to simulate conditions faced in the microreactor environment. In principle, neutron imaging, combined with in-situ sample heating, enables near real-time tracking of hydrogen diffusion in YH x on the sub-millimeter scale. In this report, an overview of neutron imaging methodology and technologies are given in the context of recent spatial measures of hydrogen concentrations in similar metal hydrides. Additionally, the commissioning and operation of a custom-built compact dual-zone furnace is given along with details on three in-situ heating measurements of YH x performed over the 2020 to 2022 LANSCE operation cycles. The aims of these experiments ranged from furnace commissioning, determining sample quality, i.e., hydrogen uniformity via neutron computed tomography, and studying the effects of applied temperature-gradients on YH x pellets. Analyses and results from these neutron imaging measurements are given along with outlooks and guidelines for optimal future hydrogen diffusion measurements. Our conclusions are as follows. Image analyses indicate that centimeter-sized yttrium hydride cylindrical pellets exhibit uniform, whole-body hydrogen desorption and absorption without clear temperature dependence as reflected in the image attenuation at the opposing ends of each sample. This suggests that despite the large magnitude in temperature gradients applied by the furnace heating elements, the sample equilibrates to an unknown intermediate temperature. The origin of this result is likely the combination of short sample length (∼1cm) and use of a TZM can for containment where the latter created a thermal short across the sample. Nevertheless, the results from the most recent measurements indicate that neither significant concentration gradients of hydrogen were formed in centimeter-sized samples through the entire temperature range (25 °C to 950 °C) nor any formed due to temperature gradients on the order of 50 °C/cm up to 700 °C/cm. Furthermore, images from the FY2021 and FY2022 measurements indicate that samples of YH x , fabricated from either the direct hydride or powder metallurgy methods, are highly uniform in their hydrogen concentration to within the measurements’ spatial resolutions. The following questions arise from these latest results: 1) What is the intermediate temperature of the pellets in the TZM cans? 2) How quickly does the temperature equilibrate within the sample? and, 3) Do the observed changes in image attenuation follow known pressure-composition-temperature relations of yttrium hydride?

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

PHYSICAL STRUCTURE OF THE PLANETARY NEBULA NGC 3242 FROM THE HOT BUBBLE TO THE NEBULAR ENVELOPE

One key feature of the interacting stellar winds model of the formation of planetary nebulae (PNe) is the presence of shock-heated stellar wind confined in the central cavities of PNe. This so-called hot bubble should be detectable in X-rays. Here we present XMM-Newton observations of NGC 3242, a multiple-shell PN whose shell morphology is consistent with the interacting stellar winds model. Diffuse X-ray emission is detected within its inner shell with a plasma temperature of {approx}2.35 x 10{sup 6} K and an intrinsic X-ray luminosity of {approx}2 x 10{sup 30} erg s{sup -1} at the adopted distance of 0.55 kpc. The observed X-ray temperature and luminosity are in agreement with 'ad hoc' predictions of models including heat conduction. However, the chemical abundances of the X-ray-emitting plasma seem to imply little evaporation of cold material into the hot bubble, whereas the thermal pressure of the hot gas is unlikely to drive the nebular expansion as it is lower than that of the inner shell rim. These inconsistencies are compounded by the apparent large filling factor of the hot gas within the central cavity of NGC 3242.

79 ASTRONOMY AND ASTROPHYSICS↗

Status of Modifications to the AGR-3/4 Fission Product Transport Model

A one-dimensional (1D) finite-element model of the AGR 3/4 experiment fission product distributions based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework has been revised and implemented within BISON, with improvements which allow for detailed temperature histories to be used for analyses. Empirically determined concentration profiles and ring inventories from the AGR-3/4 experiment have been compiled for comparison purposes. Ongoing efforts to build fission product diffusion models which can describe the applicable physics for isotopic concentration profiles in each AGR-3/4 capsule are discussed, as well as the suitability of each isotopic concentration profile in each capsule for obtaining transport parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parallel computation safety analysis irradiation targets fission product molybdenum in neutronic aspect using the successive over-relaxation algorithm

One of the research activities in support of commercial radioisotope production program is a safety research on target FPM (Fission Product Molybdenum) irradiation. FPM targets form a tube made of stainless steel which contains nuclear-grade high-enrichment uranium. The FPM irradiation tube is intended to obtain fission products. Fission materials such as Mo{sup 99} used widely the form of kits in the medical world. The neutronics problem is solved using first-order perturbation theory derived from the diffusion equation for four groups. In contrast, Mo isotopes have longer half-lives, about 3 days (66 hours), so the delivery of radioisotopes to consumer centers and storage is possible though still limited. The production of this isotope potentially gives significant economic value. The criticality and flux in multigroup diffusion model was calculated for various irradiation positions and uranium contents. This model involves complex computation, with large and sparse matrix system. Several parallel algorithms have been developed for the sparse and large matrix solution. In this paper, a successive over-relaxation (SOR) algorithm was implemented for the calculation of reactivity coefficients which can be done in parallel. Previous works performed reactivity calculations serially with Gauss-Seidel iteratives. The parallel method can be used to solve multigroup diffusion equation system and calculate the criticality and reactivity coefficients. In this research a computer code was developed to exploit parallel processing to perform reactivity calculations which were to be used in safety analysis. The parallel processing in the multicore computer system allows the calculation to be performed more quickly. This code was applied for the safety limits calculation of irradiated FPM targets containing highly enriched uranium. The results of calculations neutron show that for uranium contents of 1.7676 g and 6.1866 g (× 10{sup 6} cm{sup −1}) in a tube, their delta reactivities are the still within safety limits; however, for 7.9542 g and 8.838 g (× 10{sup 6} cm{sup −1}) the limits were exceeded.

07 ISOTOPE AND RADIATION SOURCES↗

Localization landscape theory of disorder in semiconductors. I. Theory and modeling

Here, we present here a model of carrier distribution and transport in semiconductor alloys accounting for quantum localization effects in disordered materials. This model is based on the recent development of a mathematical theory of quantum localization which introduces for each type of carrier a spatial function called localization landscape. These landscapes allow us to predict the localization regions of electron and hole quantum states, their corresponding energies, and the local densities of states. We show how the various outputs of these landscapes can be directly implemented into a drift-diffusion model of carrier transport and into the calculation of absorption/emission transitions. This creates a new computational model which accounts for disorder localization effects while also capturing two major effects of quantum mechanics, namely, the reduction of barrier height (tunneling effect) and the raising of energy ground states (quantum confinement effect), without having to solve the Schrödinger equation. Finally, this model is applied to several one-dimensional structures such as single quantum wells, ordered and disordered superlattices, or multiquantum wells, where comparisons with exact Schrödinger calculations demonstrate the excellent accuracy of the approximation provided by the landscape theory.

36 MATERIALS SCIENCE↗

Modeling the Effects of Drift Orbit Bifurcation on the Magnetopause Shadowing Loss of Radiation Belt Electrons

Abstract Last closed drift shell (LCDS) has been identified as a crucial parameter for investigating the magnetopause shadowing loss of radiation belt electrons. However, drift orbit bifurcation (DOB) effects have not been physically incorporated into the LCDS calculation. Here we calculate event‐specific LCDS using different approaches to dealing with the DOB effects, that is, tracing field lines ignoring DOB, tracing test particles rejecting field lines with DOB, and tracing particles including field lines with DOB, and then incorporate them into a radial diffusion model to simulate the fast electron dropout observed by Van Allen Probes in May 2017. The model effectively captures the fast dropout at high L* and exhibits the best agreement with data when LCDS is calculated by tracing test particles with DOB more physically included. This study represents the first quantitative modeling of the DOB effects on radiation belt magnetopause shadowing loss via a more physical specification of LCDS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nonisothermal theory of the positive column of an electric discharge in the axial magnetic field

A nonisothermal model of the positive column allowing for electron energy balance is analyzed. The influence of the axial magnetic field on the characteristics of the cylindrical positive column of a low-pressure discharge is investigated in the hydrodynamic approximation. It is shown that the magnetic field affects the plasma density distribution, plasma velocity, and electron energies. The radial dependences of the plasma density, electron energy, and plasma velocity, as well as the azimuthal velocities of electrons and ions, are calculated for helium at different values of the magnetic field strength. It is established that inertia should be taken into account in the equations for the azimuthal motion of electrons and ions. The results obtained in the hydrodynamic approximation differ significantly from those obtained in the framework of the common diffusion model of the positive column in the axial magnetic field. It is shown that the distributions of the plasma density and radial plasma velocity in the greater part of the positive column tend to those obtained in the diffusion approximation at higher values of the axial magnetic field and gas density, although substantial differences remain in the near-wall region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗