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 991 records · Page 55

Social Sensor Analytics: Making Sense of Network Models in Social Media

Social networks can be thought of as noisy sensor networks mapping real world information to the web. Owing to the extensive body of literature in sensor network analysis, this work sought to apply several novel and traditional methods in sensor network analysis for the purposes of efficiently interrogating social media data streams from raw data. We carefully revisit our definition of a social media signal from previous work both in terms of time-varying features within the data and the networked nature of the medium. Further, we detail our analysis of global patterns in Twitter over the months of November 2013 and June 2014, detect and categorize events, and illustrate how these analyses can be used to inform graph-based models of Twitter, namely using a recent network influence model called PhySense: similar to PageRank but tuned to behavioral analysis by leveraging a sociologically inspired probabilistic model. We ultimately identify forms of information dissemination via analysis of time series and dynamic graph spectra and corroborate these findings through manual investigation of the data as a requisite step in modeling the diffusion process with PhySense. We hope to sufficiently characterize global behavior in a medium such as Twitter as a means of learning global model parameters one may use to predict or simulate behavior on a large scale. We have made our time series and dynamic graph analytical code available via a GitHub repository https://github.com/cpatdowling/salsa and our data are available upon request.

event detection↗

Detailed high-resolution three-dimensional simulations of OMEGA separated reactants inertial confinement fusion experiments

Here, we present results from the comparison of high-resolution three-dimensional (3D) simulations with data from the implosions of inertial confinement fusion capsules with separated reactants performed on the OMEGA laser facility. Each capsule, referred to as a “CD Mixcap,” is filled with tritium and has a polystyrene (CH) shell with a deuterated polystyrene (CD) layer whose burial depth is varied. In these implosions, fusion reactions between deuterium and tritium ions can occur only in the presence of atomic mix between the gas fill and shell material. The simulations feature accurate models for all known experimental asymmetries and do not employ any adjustable parameters to improve agreement with experimental data. Simulations are performed with the RAGE radiation-hydrodynamics code using an Implicit Large Eddy Simulation (ILES) strategy for the hydrodynamics. We obtain good agreement with the experimental data, including the DT/TT neutron yield ratios used to diagnose mix, for all burial depths of the deuterated shell layer. Additionally, simulations demonstrate good agreement with converged simulations employing explicit models for plasma diffusion and viscosity, suggesting that the implicit sub-grid model used in ILES is sufficient to model these processes in these experiments. In our simulations, mixing is driven by short-wavelength asymmetries and longer-wavelength features are responsible for developing flows that transport mixed material towards the center of the hot spot. Mix material transported by this process is responsible for most of the mix (DT) yield even for the capsule with a CD layer adjacent to the tritium fuel. Consistent with our previous results, mix does not play a significant role in TT neutron yield degradation; instead, this is dominated by the displacement of fuel from the center of the implosion due to the development of turbulent instabilities seeded by long-wavelength asymmetries. Through these processes, the long-wavelength asymmetries degrade TT yield more than the DT yield and thus bring DT/TT neutron yield ratios into agreement with experiment. Finally, we present a detailed comparison of the flows in 2D and 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Current-induced spin torque resonance of magnetic insulators affected by field-like spin-orbit torques and out-of-plane magnetizations

The spin-torque ferromagnetic resonance (ST-FMR) in a bilayer system consisting of a magnetic insulator such as Y{sub 3}Fe{sub 5}O{sub 12} and a normal metal with spin-orbit interaction such as Pt is addressed theoretically. We model the ST-FMR for all magnetization directions and in the presence of field-like spin-orbit torques based on the drift-diffusion spin model and quantum mechanical boundary conditions. ST-FMR experiments may expose crucial information about the spin-orbit coupling between currents and magnetization in the bilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling the Formation of Hydrate-Filled Veins in Fine-Grained Sediments from in Situ Microbial Methane

Continental margin sediments are dominantly fine-grained silt and clay, and methane hydrates in these sediments are often found in semi-vertical veins and fractures. In several instances, these hydrate veins occupy discrete depth intervals that are a few tens of meters thick and are surrounded by hydrate-free sediments. As they are not connected with gas sources beneath the base of the gas hydrate stability zone (GHSZ), these isolated hydrate-bearing intervals have been interpreted as formed by in situ microbial methane. To investigate the formation of these hydrate deposits, we applied a time-dependent advection-diffusion-reaction model that includes the effects of sedimentation, compaction, solute diffusion, and microbial methane generation. Microbial methane generation depends on the amount of metabolizable organic carbon deposited at the seafloor, whose progressive degradation produces methane beneath the sulfate reduction zone. If the amount of organic carbon entering the methanogenic zone is kept constant in time, we found that the computed amounts of hydrate formed in discrete intervals within the GHSZ are well below those estimated from observations. On the other hand, if the deposition of organic carbon is higher in a given time interval, methane generation during burial is more intense in the corresponding sediment interval, resulting in enhanced hydrate formation. With variations in organic carbon deposition comparable to those generally observed in continental margins, our model was able to reproduce the methane hydrate contents that were estimated from drilling. These results support the suggestion that in situ microbial generation associated with transient organic carbon deposition is the source of methane that forms isolated intervals of hydrate-filled veins in fine-grained sediments.

58 GEOSCIENCES↗

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion↗

Strengthening model development and effects of low diffusing solutes to coarsening resistance in aluminum alloys

A modified Orowan strengthening model is proposed to account for finite rod-shaped precipitates with hemispherical caps in aluminum alloy systems. A combined computational and experimental approach is used to study the influences of anisotropic Orowan looping and solute-dislocation interaction on temperature-dependent yield strength. Here, the strengthening model is validated with a dataset containing 297 experimental precipitate geometries, chemistries, temperatures, and strength measurements, and achieves a strong predictive correlation of 0.8713 with experimentally measured yield strengths. Under conditions that are applicable to coarsening, constant particle volume fraction and aspect ratio, the model predicts that short rod precipitates provide far superior strengthening effects compared to plate precipitates. A cast Al-Si-Mg-Cu alloy with rod-shaped Q-phase (Al 3 Cu 2 Mg 9 Si 7 ) precipitates was developed with a novel chemistry exploring the use of low-diffusivity elements (Mn, Ni, V, Zr) to limit precipitate coarsening. The thermodynamic behavior of Mn, Ni, V, and Zr across the Q-phase interface is examined using transmission electron microscopy (TEM), first-principles density-functional theory (DFT) calculations. and atom-probe tomography (APT). DFT calculations utilizing TEM identified Q-phase/Al-matrix interfaces show that Mn, Ni, V, and Zr preferentially segregate to the Q-phase precipitate boundaries which suggests inhibition of precipitate coarsening and higher strengths after temperature exposure. Atom-probe tomography confirms solute atom partitioning/segregation at the Q-phase/Al-matrix interface, found in the modified commercial AS7GU alloy (A356 +0.5%Cu), which supports these observations.

36 MATERIALS SCIENCE↗

Data‐Efficient Generation of Synthetic Microstructures of Polymer‐Bonded Energetic Material With Fine‐Tuned Stable Diffusion

Among current deep learning approaches for synthetic image generation, diffusion-based models stand out in terms of algorithmic stability and ability to retain high-fidelity image features with detailed resolution. Here, in this work, we employ Dreambooth, a method for fine-tuning Stable Diffusion, on X-ray CT images of microstructure of the polymer-bonded form (PBX) of a commonly used high explosive, Pentaerythritol tetranitrate (PETN), which yields generative models for creating synthetic PBX images. The models developed here represent five classes (or ‘lots’) of microstructures and demonstrate successful generation of images of each class with high fidelity, as verified by computed classification accuracy of ∼ 94% or higher. Data augmentation afforded by such image synthesis can be used to more reliably decipher underlying statistics, build processing-structure correlations, recognize off-normal structural anomalies, and identify age-related changes. Ideas related to converting image data into appropriate density mapping and performing mesoscale simulation or surrogate modeling of detonation are also discussed.

Dreambooth↗

Titanophosphate glasses as lithium-free nonsilicate pH-responsive glasses—Compatibility between pH responsivity and self-cleaning properties

Highlights: ► Ti{sup 3+}-containing TP glasses are lithium-free nonsilicate pH-responsive ones. ► TP glasses with a large amount of Ti{sup 3+} ions show good pH responsivity. ► TP glasses with pH responsivity and self-cleaning properties are obtained. ► pH response of TP glasses is explained by phase boundary potential model. -- Abstract: Lithium silicate-based glasses have been widely used as commercially available pH glass electrodes. It was revealed that Ti{sup 3+}-containing titanophosphate (TiO{sub 2}–P{sub 2}O{sub 5}, TP) glasses are pH-responsive as lithium-free nonsilicate glasses for the first time. The absorption coefficient at 532 nm, α{sub 532} as a measure of Ti{sup 3+} content in TP glasses increased with increasing melting temperature. TP glasses with large α{sub 532} tended to give low electrical resistivity, high pH sensitivity and the short pH response time. The first post-annealing (oxidation of Ti{sup 3+}) of TP glasses at 600–620 °C for 60–240 h resulted in the occurrence of the photo-induced hydrophilicity along with the disappearance of pH responsivity and the increase of electrical resistivity. The second post-annealing (reduction of Ti{sup 4+}) of the first post-annealed TP glasses at 600–620 °C for 48 h under vacuum recovered both pH responsivity and electrical resistivity to the level of the as-prepared TP glasses with maintaining the photo-induced hydrophilicity. Moreover, the second post-annealed TP glasses had photocatalytic activity for methylene blue (MB) comparable to commercially available self-cleaning glass. Thus, TP glasses with the compatibility between pH responsivity and self-cleaning properties were obtained by the sequential post-annealing (oxidation and reduction) of as-prepared glasses. From some circumstantial evidences, pH response of TP glasses was explained in terms of phase boundary potential model related to hopping conduction of electron from Ti{sup 3+} to Ti{sup 4+} via O{sup 2−} ion in TP glasses rather than diffusion potential model.

36 MATERIALS SCIENCE↗

HIGH-ENERGY COSMIC-RAY DIFFUSION IN MOLECULAR CLOUDS: A NUMERICAL APPROACH

The propagation of high-energy cosmic rays (CRs) through giant molecular clouds constitutes a fundamental process in astronomy and astrophysics. The diffusion of CRs through these magnetically turbulent environments is often studied through the use of energy-dependent diffusion coefficients, although these are not always well motivated theoretically. Now, however, it is feasible to perform detailed numerical simulations of the diffusion process computationally. While the general problem depends upon both the field structure and particle energy, the analysis may be greatly simplified by dimensionless analysis. That is, for a specified purely turbulent field, the analysis depends almost exclusively on a single parameter-the ratio of the maximum wavelength of the turbulent field cells to the particle gyration radius. For turbulent magnetic fluctuations superimposed over an underlying uniform magnetic field, particle diffusion depends on a second dimensionless parameter that characterizes the ratio of the turbulent to uniform magnetic field energy densities. We consider both of these possibilities and parametrize our results to provide simple quantitative expressions that suitably characterize the diffusion process within molecular cloud environments. Doing so, we find that the simple scaling laws often invoked by the high-energy astrophysics community to model CR diffusion through such regions appear to be fairly robust for the case of a uniform magnetic field with a strong turbulent component, but are only valid up to {approx}50 TeV particle energies for a purely turbulent field. These results have important consequences for the analysis of CR processes based on TeV emission spectra associated with dense molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗

THE LONGITUDINAL TRANSPORT OF ENERGETIC IONS FROM IMPULSIVE SOLAR FLARES IN INTERPLANETARY SPACE

We present a study of the longitudinal spread of energetic charged particles from a localized instantaneous compact source on the Sun. Our study utilizes a diffusive-transport model for the propagation of energetic ions in interplanetary space. We show that even for very small values of the ratio of perpendicular to parallel diffusion coefficients-a few percent-the particles spread significantly in longitude. Spatial diffusion and adiabatic energy loss of ions in the interplanetary plasma cause impulsive particle events at Earth's orbit to last a few days. In this time, the combination of transport both along and across the local Parker-spiral magnetic field and the longitudinal motion of the magnetic lines of forces rooted at the Sun as it rotates leads to substantial longitudinal transport of the particles. We show that spacecraft separated by as much as 180 Degree-Sign or more may observe events associated with compact solar sources, such as those from impulsive solar flares. Our results are qualitatively consistent with recent multi-spacecraft observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Spin-scattering rates in metallic thin films measured by ferromagnetic resonance damping enhanced by spin-pumping

We determined the spin-transport properties of Pd and Pt thin films by measuring the increase in ferromagnetic resonance damping due to spin-pumping in ferromagnetic (FM)-nonferromagnetic metal (NM) multilayers with varying NM thicknesses. The increase in damping with NM thickness depends strongly on both the spin- and charge-transport properties of the NM, as modeled by diffusion equations that include both momentum- and spin-scattering parameters. We use the analytical solution to the spin-diffusion equations to obtain spin-diffusion lengths for Pt and Pd. By measuring the dependence of conductivity on NM thickness, we correlate the charge- and spin-transport parameters, and validate the applicability of various models for momentum-scattering and spin-scattering rates in these systems: constant, inverse-proportional (Dyakanov-Perel), and linear-proportional (Elliot-Yafet). We confirm previous reports that the spin-scattering time appears to be shorter than the momentum scattering time in Pt, and the Dyakanov-Perel-like model is the best fit to the data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Surface enhanced Raman spectroscopy for microfluidic pillar arrayed separation chips

Numerous studies have addressed the challenges of implementing miniaturized microfluidic platforms for chemical and biological separation applications. However, the integration of real time detection schemes capable of providing valuable sample information under continuous, ultra low volume flow regimes has not fully been addressed. In this report we present a chip based chromatography system comprising of a pillar array separation column followed by a reagent channel for passive mixing of a silver colloidal solution into the eluent stream to enable surface enhanced Raman spectroscopy (SERS) detection. Our design is the first integrated chip based microfluidic device to combine pressure driven separation capability with real time SERS detection. With this approach we demonstrate the ability to collect distinctive SERS spectra with or without complete resolution of chromatographic bands. Computational fluidic dynamic (CFD) simulations are used to model the diffusive mixing behavior and velocity profiles of the two confluent streams in the microfluidic channels. We evaluate the SERS spectral band intensity and chromatographic efficiency of model analytes with respect to kinetic factors as well as signal acquisition rates. Additionally, we discuss the use of a pluronic modified silver colloidal solution as a means of eliminating contamination generally caused by nanoparticle adhesion to channel surfaces.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A POSSIBLE APPROACH TO THREE-DIMENSIONAL COSMIC-RAY PROPAGATION IN THE GALAXY. IV. ELECTRONS AND ELECTRON-INDUCED {gamma}-RAYS

Based on the diffusion-halo model for cosmic-ray (CR) propagation, including stochastic reacceleration due to collisions with hydromagnetic turbulence, we study the behavior of the electron component and the diffuse {gamma}-rays (D{gamma}'s) induced by them. The galactic parameters appearing in these studies are essentially the same as those appearing in the hadronic CR components, while we additionally need information on the interstellar radiation field, taking into account dependences on both the photon energy, E{sub ph}, and the position, r. We compare our numerical results with the data on hadrons, electrons, and D{gamma}'s, including the most recent results from Fermi, which gives two remarkable results: (1) the electron spectrum falls with energy as E{sup -3}{sub e} up to 1 TeV, and does not exhibit prominent spectral features around 500 GeV, in contrast to the dramatic excess appearing in both ATIC and PPB-BETS spectra and (2) the EGRET GeV-excess in the D{gamma} spectrum is neither due to an astronomical origin (much harder CR spectrum in the galactic center) nor a cosmological one (dark matter annihilation or decay), but due to an instrumental problem. In the present paper, however, we focus our interest rather conservatively upon the internal relation between these three components, using common galactic parameters. We find that they are in reasonable harmony with each other within both the theoretical and experimental uncertainties, apart from the electron-anomaly problem, while some enhancement of D{gamma}'s appears in the high galactic latitude with |b|>60{sup 0} in the GeV region.

79 ASTRONOMY AND ASTROPHYSICS↗

Recovering doping profiles in semiconductor devices with the Boltzmann-Poisson model

We investigate numerically an inverse problem related to the Boltzmann-Poisson system of equations for transport of electrons in semiconductor devices. The objective of the (ill-posed) inverse problem is to recover the doping profile of a device, presented as a source function in the mathematical model, from its current-voltage characteristics. To reduce the degree of ill-posedness of the inverse problem, we proposed to parameterize the unknown doping profile function to limit the number of unknowns in the inverse problem. We showed by numerical examples that the reconstruction of a few low moments of the doping profile is possible when relatively accurate time-dependent or time-independent measurements are available, even though the later reconstruction is less accurate than the former. We also compare reconstructions from the Boltzmann-Poisson (BP) model to those from the classical drift-diffusion-Poisson (DDP) model, assuming that measurements are generated with the BP model. We show that the two type of reconstructions can be significantly different in regimes where drift-diffusion-Poisson equation fails to model the physics accurately. However, when noise presented in measured data is high, no difference in the reconstructions can be observed.

97 MATHEMATICS AND COMPUTING↗

Model for Deriving the Direct Normal and Diffuse Horizontal Irradiance from the Global Tilted Irradiance

We developed a separation model for deriving the direct normal irradiance (DNI) and the diffuse horizontal irradiance (DHI) from the global tilted irradiance (GTI). The model is based on the DIRINT model, but substitutes the use of the GTI in place of the global horizontal irradiance (GHI) for determining the global clearness index (Kt). Additionally, an iterative method is used to adjust Kt to improve the derived values of the DNI and the DHI. The model is referred to as the GTI-DIRINT model. When using for model input the GTI measured with a pyranometer tilted at a small angle from the horizontal (10°), the deviations between the measured DNI and DHI and the GTI-DIRINT-modeled DNI and DHI were essentially the same as those for the DIRINT model when using the GHI for model input. Increasing the tilt angle from horizontal (as great as 44° for Eugene, Oregon) increased the deviations between the modeled and measured DNI and DHI, but our results were still reasonable. The GTI-DIRINT model provides reasonable values of DNI and DHI from the GTI, and if they are subsequently used with the Perez tilted surface model to model the GTI for orientations within 30° of the azimuth and tilt of the GTI used for model input, the result is as good or better than may be obtained when using a measured DNI and DHI, and without their significant instrumentation cost.

14 SOLAR ENERGY↗

Observation of high-energy neutrinos from the Galactic plane

The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth’s atmosphere, is unknown. Because of deflection by interstellar magnetic fields, cosmic rays produced within the Milky Way arrive at Earth from random directions. However, cosmic rays interact with matter near their sources and during propagation, which produces high-energy neutrinos. We searched for neutrino emission using machine learning techniques applied to 10 years of data from the IceCube Neutrino Observatory. By comparing diffuse emission models to a background-only hypothesis, we identified neutrino emission from the Galactic plane at the 4.5σ level of significance. The signal is consistent with diffuse emission of neutrinos from the Milky Way but could also arise from a population of unresolved point sources.

Science & Technology - Other Topics↗

Neutronics Modeling of the High Flux Isotope Reactor using COMSOL

The High Flux Isotope Reactor located at the Oak Ridge National Laboratory is a versatile 85 MWth research reactor with cold and thermal neutron scattering, materials irradiation, isotope production, and neutron activation analysis capabilities. HFIR staff members are currently in the process of updating the thermal hydraulic and reactor transient modeling methodologies. COMSOL Multiphysics has been adopted for the thermal hydraulic analyses and has proven to be a powerful finite-element-based simulation tool for solving multiple physics-based systems of partial and ordinary differential equations. Modeling reactor transients is a challenging task because of the coupling of neutronics, heat transfer, and hydrodynamics. This paper presents a preliminary COMSOL-based neutronics study performed by creating a two-dimensional, two-group, diffusion neutronics model of HFIR to study the spatially-dependent, beginning-of-cycle fast and thermal neutron fluxes. The 238-group ENDF/B-VII neutron cross section library and NEWT, a two-dimensional, discrete-ordinates neutron transport code within the SCALE 6 code package, were used to calculate the two-group neutron cross sections required to solve the diffusion equations. The two-group diffusion equations were implemented in the COMSOL coefficient form PDE application mode and were solved via eigenvalue analysis using a direct (PARDISO) linear system solver. A COMSOL-provided adaptive mesh refinement algorithm was used to increase the number of elements in areas of largest numerical error to increase the accuracy of the solution. The flux distributions calculated by means of COMSOL/SCALE compare well with those calculated with benchmarked three-dimensional MCNP and KENO models, a necessary first step along the path to implementing two- and three-dimensional models of HFIR in COMSOL for the purpose of studying the spatial dependence of transient-induced behavior in the reactor core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗