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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Structure and behavior of water-laden CH 4 /air counterflow diffusion flames

A counterflow configuration was used to measure thermal and species structure in water-vapor diluted nonpremixed methane–air flames. The motivation is to understand the chemical and thermal effects that water has when it is introduced as a diluent into the fuel side. Here, this work is relevant to combustion processes where water is incorporated naturally in the fuel; e.g., methane hydrates, and when water is added intentionally for emission reduction such as in flares and H 2 O/fuel emulsions combustion. Experimental data are compared to 1-D computations. The agreement is generally very good, but the one dimensional counterflow diffusion model overpredicts flame temperature and major radical, OH, concentration very near extinction in highly diluted H 2 O–methane/air diffusion flames. Changes in flame position, flame width, and peak temperature with the addition of water were measured. Flame temperatures were measured with thin filament pyrometry. OH-PLIF is used to characterize the flame reaction zone with water dilution; the OH distribution, flame position and thickness from the OH-PLIF images were measured. The results show that the OH intensity and reaction zone thickness decreases with the increase in water. Predictions and experiments demonstrate that water mainly acts thermally to lower the flame temperature until extinction. The OH maximum intensity shifts towards the air side of the counterflow burner with water addition. OH is also measured with CO 2 dilution of the fuel stream, and the results are compared with H 2 O addition, including comparisons with the OH molar peak predictions obtained using the GRI 3.0 mechanism and the CHEMKIN Pro one dimensional counterflow model. The study indicates that water’s chemical effects are to change the production and depletion of OH, H and O radicals, especially near extinction. Chemical kinetics simulation of the flame demonstrates good agreement in OH and flame temperatures over a wide range of dilution away from extinction, particularly for CO 2 . Lastly, an over prediction of the water carrying capacity near extinction is found for highly water-diluted flames.

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

Hybrid stochastic simulation of reaction-diffusion systems with slow and fast dynamics

In this paper, we present a novel hybrid method to simulate discrete stochastic reaction-diffusion models arising in biochemical signaling pathways. We study moderately stiff systems, for which we can partition each reaction or diffusion channel into either a slow or fast subset, based on its propensity. Numerical approaches missing this distinction are often limited with respect to computational run time or approximation quality. We design an approximate scheme that remedies these pitfalls by using a new blending strategy of the well-established inhomogeneous stochastic simulation algorithm and the tau-leaping simulation method. The advantages of our hybrid simulation algorithm are demonstrated on three benchmarking systems, with special focus on approximation accuracy and efficiency.

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

NEAMS-ATF M3 Milestone Report: Literature Review of Modeling of Radiation-Induced Swelling in Fe-Cr-Al Steels

Fe-Cr-Al steels are proposed as accident-tolerant-fuel (ATF) cladding materials in light water reactors due to their excellent oxidation resistance at high temperatures. Currently, the understanding of their performance in reactor environment is still limited. In this review, firstly we reviewed the experimental studies of Fe-Cr-Al based alloys with particular focus on the radiation effects in these alloys. Although limited data are available in literature, several previous and recent experimental studies have shown that Fe-Cr-Al based alloys have very good void swelling resistance at low and moderate irradiation doses but the growth of dislocation loops is very active. Overall, the behavior of radiation damage evolution is similar to that in Fe-Cr ferritic/martensitic alloys. Secondly, we reviewed the rate theory-based modeling methods for modeling the coevolution of voids and dislocation loops in materials under irradiation such as Frenkel pair three-dimensional diffusion model (FP3DM) and cluster dynamics. Finally, we summarized and discussed our review and proposed our future plans for modeling radiation damage in Fe-Cr-Al based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2013

The AN neutron transport by nodal diffusion

The two group diffusion model combined to a nodal approach in space is the preferred scheme for the industrial simulation of nuclear water reactors. The main selling point is the speed of computation, allowing a large number of parametric studies. Anyway, the drawbacks of the underlying diffusion equation may arise with highly heterogeneous interfaces, often encountered in modern UO{sub 2} and MO{sub x} fuel loading patterns, and boron less controlled systems. This paper aims at showing how the simplified AN transport model, equivalent to the well known SPN, can be implemented in standard diffusion codes with minor modifications. Some numerical results are illustrated. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Modeling MHD accretion-ejection—from the launching area to propagation scales

We present the results of axisymmetric magnetohydrodynamic simulations investigating the launching of jets and outflows from a magnetically diffusive accretion disk. The time evolution of the disk structure is self-consistently taken into account. In contrast to previous works, we have applied spherical coordinates for the numerical grid, implying substantial benefits concerning the numerical resolution and the stability of the simulation. Due to the new setup, we were able to run simulations for more than 150,000 dynamical times on a domain extending 1500 inner disk radii with a resolution of up to 24 cells per disk height in the inner disk. Depending on the disk magnetization, jet launching occurs in two different but complementary regimes—jets driven predominantly by centrifugal or magnetic forces. These regimes differ in the ejection efficiency concerning mass, energy, and angular momentum. We show that it is the actual disk magnetization and not so much the initial magnetization that best describes the disk-jet evolution. Considering the actual disk magnetization, we also find that simulations starting with different initial magnetization evolve in a similar—typical—way, due to advection and diffusion, as the magnetic flux in the disk evolves in time. Exploring a new, modified diffusivity model, we confirm the self-similar structure of the global jet-launching disk, obtaining power laws for the radial profiles of the disk physical variables such as density, magnetic field strength, or accretion velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

High power density soft x-ray GaAs photodiodes with tailored spectral response

GaAs photodiode arrays have been designed for non-destructive monitoring of x-ray beam position in soft coherent beamline front ends in synchrotron light sources. A shallow p-on-n junction was employed to reduce the device photocurrent density to optimize the operation with beam power densities up to 20 W mm -2 , mainly coming from hard x-rays. With this approach, the photocurrent is primarily defined by the excess carriers generated by low-energy x-ray photons absorbed near the detector surface. The p–n junction structures were grown by molecular beam epitaxy and processed into 64-element linear photodiode arrays. The devices were characterized first in the visible range with a high-power Ar-ion laser and then tested in the soft and hard x-ray regions up to 10 keV at two beamlines of the National Synchrotron Light Source II. The responsivity was measured to be 0.16 A W -1 at 0.7 keV and 0.05 A W -1 at 6 keV in agreement with modeling. At higher x-ray energies the measured responsivity was lower than predicted in the framework of the carrier diffusion model; a possible explanation is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Interstitial and interlayer ion diffusion geometry extraction in graphitic nanosphere battery materials

Large-scale molecular dynamics (MD) simulations are commonly used for simulating the synthesis and ion diffusion of battery materials. A good battery anode material is determined by its capacity to store ion or other diffusers. However, modeling of ion diffusion dynamics and transport properties at large length and long time scales would be impossible with current MD codes. To analyze the fundamental properties of these materials, therefore, we turn to geometric and topological analysis of their structure. In this paper, we apply a novel technique inspired by discrete Morse theory to the Delaunay triangulation of the simulated geometry of a thermally annealed carbon nanosphere. We utilize our computed structures to drive further geometric analysis to extract the interstitial diffusion structure as a single mesh. Lastly, our results provide a new approach to analyze the geometry of the simulated carbon nanosphere, and new insights into the role of carbon defect size and distribution in determining the charge capacity and charge dynamics of these carbon based battery materials.

25 ENERGY STORAGE↗
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2021

The Role of Oxygen Concentration in Enabling High Gradients in Niobium SRF Cavities

We studied the role of O concentration with depth in the performance of Nb SRF cavities. An ensemble of electropolished 1.3 GHz cavities, which initially showed high field Q-slope (HFQS), was subjected to sequential testing and treatment with in-situ low temperature baking at various temperatures. We find that increasing the bake duration causes (i) an increase in the onset of HFQS until it is absent up to quench (ii) a non-monotonic relationship with the quench field (iii) an evolution of the RBCS toward a non-equilibrium behavior that drives anti-Q slope. Our data is qualitatively explained by assuming an O diffusion model and suggests that the mitigation of HFQS that arises from 120 C in-situ LTB is mediated by the diffusion of O from the native oxide which prevents the precipitation of proximity-coupled Nb nano-hydrides, in turn enabling higher quench fields. The decrease in quench field for cavities in which O has been diffused >90 nm from the RF surface may be due to a reduction of the field limit in the SS bilayer structure. We also suggest that the evolution of the RBCS occurs due to the absence of proximity coupled inclusions, bringing about non-equilibrium effects.

Bafia, D.↗
DOE OSTI2025

MoRE-Brain: Routed Mixture of Experts for Interpretable and Generalizable Cross-Subject fMRI Visual Decoding

Decoding visual experiences from fMRI offers a powerful avenue to understand human perception and develop advanced brain-computer interfaces. However, current progress often prioritizes maximizing reconstruction fidelity while overlooking interpretability, an essential aspect for deriving neuroscientific insight. To address this gap, we propose MoRE-Brain, a neuro-inspired framework designed for high-fidelity, adaptable, and interpretable visual reconstruction. MoRE-Brain uniquely employs a hierarchical Mixture-of-Experts architecture where distinct experts process fMRI signals from functionally related voxel groups, mimicking specialized brain networks. The experts are first trained to encode fMRI into the frozen CLIP space. A finetuned diffusion model then synthesizes images, guided by expert outputs through a novel dual-stage routing mechanism that dynamically weighs expert contributions across the diffusion process. MoRE-Brain offers three main advancements: First, it introduces a novel Mixture-of-Experts architecture grounded in brain network principles for neuro-decoding. Second, it achieves efficient cross-subject generalization by sharing core expert networks while adapting only subject-specific routers. Third, it provides enhanced mechanistic insight, as the explicit routing reveals precisely how different modeled brain regions shape the semantic and spatial attributes of the reconstructed image. Extensive experiments validate MoRE-Brain’s high reconstruction fidelity, with bottleneck analyses further demonstrating its effective utilization of fMRI signals, distinguishing genuine neural decoding from over-reliance on generative priors. Consequently, MoRE-Brain marks a substantial advance towards more generalizable and interpretable fMRI-based visual decoding.

Wei, Yuxiang [Georgia Institute of Technology]↗
DOE OSTI2013

Estimation of Sorption Behavior of Europium(III) Using Biotite Flakes - 13272

The interaction of biotite and Eu(III) (europium (III)) was examined by using secondary ion-microprobe mass spectrometer (SIMS), fluorescence emission spectrum and decay behavior of fluorescence emission spectrum in addition to the time-changes of Eu(III) and potassium ions concentrations in a solution, using the flake form samples. The results of SIMS showed that the intensity of Eu was gradually decreasing with depth, while the intensity of Eu in the case shaken for 30 days exceeded that in the case for 1 day. Furthermore, the spatial distribution of Eu(III) and potassium ions in the flake of biotite suggested that Eu ions diffuse mainly from the edges of biotite flake, while Eu ions can slightly diffuse through some small cracks existing on the flake surface far from the edges. Besides, the elution amount of potassium from the biotite flakes into a solution was proportional to the sorption amount of Eu(III). The changes nearly revealed ion exchange between these ions, while muscovite flake sample did not show such ion exchange reaction. In addition, from the time-change of Eu(III) concentration, an apparent diffusion coefficient was estimated to be 8.0x10{sup -12} m{sup 2}/s, by using two-dimensional diffusion model coupled with a film between the solid phase and the liquid phase. Furthermore, the fluorescent intensity decreased with the shaking (contacting) time. This means that Eu(III) gradually diffuses into the inside of biotite edges of the biotite flakes, after the sorption of Eu(III) in the edges. This tendency was observed also in the powder samples. The observed fluorescence decay (at 592 nm in wave length) showed almost similar curve in any samples, indicating a certain sorption form of Eu(III) onto the edges of the biotite flakes. These results mentioned above suggest that the diffusion processes through internal layer in biotite mainly control the sorption behavior of multivalent ions. Such diffusion processes affect the retardation-effects on fracture surfaces in the rock matrix, depending on the fluid flow velocity of groundwater. That is, a more reliable model considering the mass transfer in the internal layer of biotite may be required to estimate the sorption behavior of RNs with biotite which controls the whole sorption behavior of granite. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Kinetics of calcium sulfoaluminate formation from tricalcium aluminate, calcium sulfate and calcium oxide

The formation kinetics of tricalcium aluminate (C{sub 3}A) and calcium sulfate yielding calcium sulfoaluminate (C{sub 4}A{sub 3}$) and the decomposition kinetics of calcium sulfoaluminate were investigated by sintering a mixture of synthetic C{sub 3}A and gypsum. The quantitative analysis of the phase composition was performed by X-ray powder diffraction analysis using the Rietveld method. The results showed that the formation reaction 3Ca{sub 3}Al{sub 2}O{sub 6} + CaSO{sub 4} → Ca{sub 4}Al{sub 6}O{sub 12}(SO{sub 4}) + 6CaO was the primary reaction < 1350 °C with and activation energy of 231 ± 42 kJ/mol; while the decomposition reaction 2Ca{sub 4}Al{sub 6}O{sub 12}(SO{sub 4}) + 10CaO → 6Ca{sub 3}Al{sub 2}O{sub 6} + 2SO{sub 2} ↑ + O{sub 2} ↑ primarily occurred beyond 1350 °C with an activation energy of 792 ± 64 kJ/mol. The optimal formation region for C{sub 4}A{sub 3}$ was from 1150 °C to 1350 °C and from 6 h to 1 h, which could provide useful information on the formation of C{sub 4}A{sub 3}$ containing clinkers. The Jander diffusion model was feasible for the formation and decomposition of calcium sulfoaluminate. Ca{sup 2+} and SO{sub 4}{sup 2−} were the diffusive species in both the formation and decomposition reactions. -- Highlights: •Formation and decomposition of calcium sulphoaluminate were studied. •Decomposition of calcium sulphoaluminate combined CaO and yielded C{sub 3}A. •Activation energy for formation was 231 ± 42 kJ/mol. •Activation energy for decomposition was 792 ± 64 kJ/mol. •Both the formation and decomposition were controlled by diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Protein–Protein and Protein–Ligand Interactions by the ClusPro Team in CASP16

ABSTRACT In the CASP16 experiment, our team employed hybrid computational strategies to predict both protein–protein and protein–ligand complex structures. For protein–protein docking, we combined physics‐based sampling—using ClusPro FFT docking and molecular dynamics—with AlphaFold (AF)‐based sampling, followed by AF‐based refinement. Our method produced numerous high‐accuracy complex models, including cases where AF alone failed, underscoring the critical role of physics‐based sampling alongside deep learning‐based refinement. For protein–ligand docking, we integrated the ClusPro LigTBM template‐based approach with a machine learning‐based confidence model for rescoring. The method preserves conserved interaction fragments derived from homologous complexes, followed by local resampling using physics‐based sampling and a diffusion model. Our template‐based strategy achieved a mean lDDT‐PLI of 0.69 across 233 targets, which was highly competitive. These results demonstrate that combining physics‐based modeling with AI‐driven refinement can significantly enhance the accuracy of both protein–protein and protein–ligand structure predictions.

Ashizawa, Ryota [Department of Applied Mathematics↗
DOE OSTI2019

Equilibrium, Kinetic and Thermodynamic Studies on Methylene Blue Adsorption by Konjac Glucomannan/Activated Carbon Aerogel

Activated carbon (AC) is a common adsorbent used for the purification of aqueous solutions, yet its poor aqueous dispersion may negatively affects its practicability significantly. In the present study, the fabrication of an AC aerogel immobilized with konjac glucomannan (KGM/AC) for environmental protection was reported, and the parameters (e.g. AC content in aerogel, solution pH, adsorption duration, adsorption temperature and adsorbent dosage) were analyzed systematically using methylene blue (MB) as a model adsorbate. Further, to explore the adsorption equilibrium isotherm at equilibrium, the adsorption data were respectively fitted to Langmuir and Freundlich models, while kinetics study was conducted to explain the adsorption mechanism and predict the adsorption efficiency. Besides, the effect of temperature on the adsorption behavior was assessed using Van’t Hoff equation. According to the experimental and calculative results, the adsorption data at equilibrium could be well geared into the Langmuir model, and the adsorption kinetics behavior was better consistent with the pseudo-second-order expression as well as the intra-particle diffusion model. The Brunauer, Emmett and Teller surface area of KGM/AC aerogel was 352.44 m{sup 2}/g. KGM/AC aerogel exhibited the maximal capacity of adsorption of 416.67 mg/g at 323 K, and its adsorption behavior toward MB was endothermic and spontaneous in nature.

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

Dynamical properties of the hydration shell of fully deuterated myoglobin

Freeze-dried perdeuterated sperm whale myoglobin was kept in a water-saturated atmosphere in order to obtain a hydration degree of 335 {sup 1}H{sub 2}O molecules per one myoglobin molecule. Incoherent neutron scattering was performed at the neutron spectrometer TOFTOF at the FRM II in an angular range of q from 0.6 to 1.8 A{sup -1} and a temperature range from 4 to 297 K. We used neutrons with a wavelength of {lambda}{alpha}E 6 A and an energy resolution of about 65 {mu}eV corresponding to motions faster than 10 ps. At temperatures above 225 K, broad lines appear in the spectra caused by quasielastic scattering. For an explanation of these lines, we assumed that there are only two types of protons, those that are part of the hydration water (72%) and those that belong to the protein (28%). The protons of the hydration water were analyzed with the diffusion model of Singwi and Sjoelander [Phys. Rev. 119, 863 (1960)]. In this model, a water molecule stays for a time {tau}{sub 0} in a bound state performing oscillatory motions. Thereafter, the molecule performs free diffusion for the time {tau}{sub 1} in a nonbound state followed again by the oscillatory motions for {tau}{sub 0} and so forth. We used the general formulation with no simplifications as {tau}{sub 0}>>{tau}{sub 1} or {tau}{sub 1}>>{tau}{sub 0}. At room temperature, we obtained {tau}{sub 0} {alpha}E 104 ps and {tau}{sub 1} {alpha}E 37 ps. For the protein bound hydrogen, the dynamics is described by a Brownian oscillator where the protons perform overdamped motions in limited space.

36 MATERIALS SCIENCE↗