Modeling fast diffusion processes in time integration of stiff stochastic differential equations.
Abstract not provided.
Han, Xiaoying↗
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.
Diffuse Interface Modeling Crystal Plasticity with Accommodation of Grain Boundary Sliding - Poster
Presented by NETL at GRC for Physical Metallurgy, 07/19
Cheng, Tianle↗
Diffuse Interface Modeling of Grain Boundary Segregation in Metallic Systems.
Abstract not provided.
Abdeljawad, Fadi F.↗
Nanoscopic dynamics of bicontinous microemulsions: effect of membrane associated protein
Bicontinous microemulsions (BμE) generally consist of nanodomains formed by surfactant in a mixture of water and oil at nearly equal proportions and are potential candidates for the solubilization and purification of membrane proteins. In this paper, we present the first time report of nanoscopic dynamics of surfactant monolayers within BμEs formed by the anionic surfactant sodium dodecyl sulfate (SDS) measured on the nanosecond to picosecond time scale using quasielastic neutron scattering (QENS). BμEs investigated herein consisted of middle phases isolated from Winsor-III microemulsion systems that were formed by mixing aqueous and oil solutions under optimal conditions. QENS data indicates that surfactants undergo two distinct motions, namely (i) lateral motion along the surface of the oil nanodomains and (ii) localized internal motion. Lateral motion can be described using a continuous diffusion model, from which the lateral diffusion coefficient is obtained. Internal motion of surfactant is described using a model which assumes that a fraction of the surfactants’ hydrogens undergoes localized translational diffusion that could be considered confined within a spherical volume. The effect of cytochrome c, an archetypal membrane-associated protein known to strongly partition near the surfactant head groups in BμEs (a trend supported by small-angle X-ray scattering [SAXS] analysis), on the dynamics of BμE has also been investigated. QENS results demonstrated that cytochrome c significantly hindered both the lateral and the internal motions of surfactant. The lateral motion was more strongly affected: a reduction of the lateral diffusion coefficient by 33% was measured. This change is mainly attributable to the strong association of cytochrome c with oppositely charged SDS. In contrast, analysis of SAXS data suggested that thermal fluctuations (for a longer length and slower time scale compared to QENS) were increased upon incorporation of cytochrome c. Finally, this study demonstrates the utility of QENS for evaluating dynamics of BμEs in nanoscopic region, and that proteins directly affect the microscopic dynamics, which is of relevance for evaluating release kinetics of encapsulated drugs from BμE delivery systems and the use of BμEs as biomembrane mimetic systems for investigating membrane protein–biomembrane interactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Modeling fission product diffusion in TRISO fuel particles with BISON
Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Comparison of experimental methods for estimating matrix diffusion coefficients for contaminant transport modeling
Here, diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating effective matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of effective matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than effective matrix diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields effective matrix diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.
54 ENVIRONMENTAL SCIENCES↗
Co-ion specific effect on sodium halides sorption and transport in a cross-linked poly(p-styrene sulfonate-co-divinylbenzene) for membrane applications
Here, sodium halide sorption and transport in a cation exchange membrane based on cross-linked poly(p-styrene sulfonate-co-divinylbenzene) (CR61) were investigated experimentally and theoretically. At fixed co-ion valence, sodium halide partitioning between aqueous electrolyte solutions and ion exchange polymers is only minimally sensitive to co-ion chemical identity. NaCl, NaF and NaI exhibit similar solubility coefficients in CR61, which agrees reasonably well with predictions of the Donnan-Manning model, especially at high external salt concentrations. Ion activity coefficients in the membrane were evaluated and ion sorption thermodynamics were investigated. Ion diffusion coefficients in CR61 were predicted using a recently developed framework based on the Manning diffusion model, which accounts for electrostatic effects on ion diffusion, and the Mackie-Meares model, which accounts for tortuosity/area effects on ion diffusion in swollen polymers. Model predictions agree reasonably well with the experimental data, especially at high external salt concentrations. The observed decrease in salt diffusion coefficients in CR61 relative to values in aqueous solutions is due primarily to tortuosity/area effects. Salt diffusion coefficients correlate with co-ion hydrated radius and membrane water content. Finally, salt permeability coefficients predicted using the solution-diffusion model exhibit good agreement with experimental data, especially at external salt concentration greater than 0.1 mol/L.
42 ENGINEERING↗
Sequential Infiltration Synthesis of Multinary Metal Oxides: The Role of Precursor Diffusion and Reactivity
Understanding volatile metal–organic precursor transport in polymer matrices during sequential infiltration synthesis (SIS) is critical to the design of hybrid materials of tailored composition and uniform structure. Here, we investigate the diffusion and reaction behavior of diethylzinc (DEZ) in poly(methyl methacrylate) (PMMA) matrices with variable indium and zinc incorporation. Unlike conventional SIS where precursors adduct or react directly with the polymer backbone, DEZ shows minimal interaction with PMMA. However, DEZ may interact strongly with In–OH or Zn–OH moieties that are previously incorporated into a PMMA thick film. Using spectroscopic ellipsometry, X-ray photoelectron spectroscopy depth profiling, and EDX line scans, we quantify Zn incorporation as a function of exposure time, polymer thickness, and previous infiltration steps. The Zn distribution evolves from surface-localized to deep within the ∼1 μm thick film with increasing DEZ exposures up to 30 min. A reaction-diffusion model incorporating a hindering factor to capture diffusivity decay due to product accumulation reproduces the Zn depth profile and yields physically meaningful parameters to guide future growth. Furthermore, this study provides a framework for interpreting metal–organic precursor infiltration behavior and highlights the utility of model-guided SIS process design.
diethylzinc↗
Deep Throttle Turbopump Technology Design Concepts
The objective of this project is to increase the throttling range of turbopumps from 30 to 120% of the design value, while maintaining high performance levels. Details are given on wide flow range issues, H-Q characteristics, stall characteristics, energy levels, pressure fluctuations at impeller exit, WFR impeller characteristics, commercial diffuser pumps, slotted or tandem vanes, leading edge characteristics, leading edge models, throat models, diffusion passage models, computational fluid dynamics (CFD) methodologies, and CFD flow cases.
Guinzburg, Adiel↗
Salt and ion transport in a series of crosslinked AMPS/PEGDA hydrogel membranes
Salt permeability and ionic conductivity are critical properties for membranes in water purification and energy applications. Both properties depend on individual ion sorption and diffusion coefficients, which are significantly influenced by polymer chemical and physical parameters such as fixed charge concentration and membrane water content. However, systematic studies connecting polymer structure to ion transport properties are still lacking. In this study, a series of uncharged and charged membranes were synthesized using poly(ethylene glycol) diacrylate (PEGDA) as a cross-linker and 2-acrylamido-2-methyl-1-propanesulfonic acid (AMPS) as a charged monomer. Membrane fixed charge concentration and water uptake were systematically varied by adjusting AMPS content in the pre-polymerization mixture. Salt sorption and permeability coefficients and ionic conductivity of these membranes were measured as a function of NaCl solution concentration (0.01–1 M). Combining the solution-diffusion model and Nernst-Planck equation, individual ion diffusion coefficients were calculated. Experimental Na + diffusion coefficients for all materials were well described by the Mackie and Meares tortuosity model, highlighting the strong influence of water content on ion diffusivity in both uncharged and charged polymers. Model predictions for Cl – diffusion coefficients agree reasonably well with experimental values, with some deviations occurring in more highly charged membranes. Finally, this discrepancy might result from interactions not captured by the Mackie and Meares model (e.g., fixed charge-ion interactions).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Flow matching meets biology and life science: a survey
Over the past decade, advances in generative modeling, such as generative adversarial networks, masked autoencoders, and diffusion models, have significantly transformed biological research and discovery, enabling breakthroughs in molecule design, protein generation, catalysis discovery, drug discovery, and beyond. At the same time, biological applications have served as valuable testbeds for evaluating the capabilities of generative models. Recently, flow matching has emerged as a powerful and efficient alternative to diffusion-based generative modeling, with growing interest in its application to problems in biology and life sciences. This paper presents the first comprehensive survey of recent developments in flow matching and its applications in biological domains. We begin by systematically reviewing the foundations and variants of flow matching, and then categorize its applications into three major areas: biological sequence modeling, molecule generation and design, and peptide and protein generation. For each, we provide an in-depth review of recent progress. We also summarize commonly used datasets and software tools, and conclude with a discussion of potential future directions.
59 BASIC BIOLOGICAL SCIENCES↗
Ge auto-doping and out-diffusion in InGaP grown on Ge substrate and their effects on the ordering of InGaP
We report on the Ge auto-doping and out-diffusion in InGaP epilayer with Cu-Pt ordering grown on 4-in. Ge substrate. Ge profiles determined from secondary ion mass spectrometry indicate that the Ge out-diffusion depth is within 100 nm. However, the edge of the wafer suffers from stronger Ge gas-phase auto-doping than the center, leading to ordering deterioration in the InGaP epilayer. In the edge, we observed a residual Cu-Pt ordering layer left beneath the surface, suggesting that the ordering deterioration takes place after the deposition rather than during the deposition and In/Ga inter-diffusion enhanced by Ge vapor-phase auto-doping is responsible for the deterioration. We thus propose a di-vacancy diffusion model, in which the amphoteric Ge increases the di-vacancy density, resulting in a Ge density dependent diffusion. In the model, the In/Ga inter-diffusion and Ge out-diffusion are realized by the random hopping of In/Ga host atoms and Ge atoms to di-vacancies, respectively. Simulation based on this model well fits the Ge out-diffusion profiles, suggesting its validity. By comparing the Ge diffusion coefficient obtained from the fitting and the characteristic time constant of ordering deterioration estimated from the residual ordering layer, we found that the hopping rates of Ge and the host atoms are in the same order of magnitude, indicating that di-vacancies are bound in the vicinity of Ge atoms.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Kinematic models of supergranular diffusion on the Sun
We develop kinematic models of diffusion generated by supergranulation at the solar surface. These models use current observations for the size, horizontal velocity, and lifetime of supergranules. Because there is no observational description of the appearance and disappearance of supergranules, we investigate models using several plausible evolution processes, including the effect of different lifetime distribution functions for the cells. The results are quite insensitive to the methods chosen to replace old supergranules, the distribution of cell lifetimes, and even the cell lifetime itself, for mean lifetimes between 15 to 30 hr. Calculated diffusion coefficients range between 500 and 700 sq km/s, in agreement with the best fit diffusion coefficients used by Sheeley and his collaborators to model the large-scale distribution of magnetic fields over the solar surface. However, our models do not explain the field distribution in plage, and they predict that virtually all the strong field in quiet Sun exists in relatively isolated clumps. We suggest possible mechanisms for the creation of plage and the bright network seen in quiet Sun.
Simon, G. W.↗
Room-Temperature Micron-Scale Exciton Migration in a Stabilized Emissive Molecular Aggregate
We report 1.6 ± 1 μm exciton transport in self-assembled supramolecular light-harvesting nanotubes (LHNs) assembled from amphiphillic cyanine dyes. We stabilize LHNs in a sucrose glass matrix, greatly reducing light and oxidative damage and allowing the observation of exciton–exciton annihilation signatures under weak excitation flux. Fitting to a one-dimensional diffusion model, we find an average exciton diffusion constant of 55 ± 20 cm 2 /s, among the highest measured for an organic system. We develop a simple model that uses cryogenic measurements of static and dynamic energetic disorder to estimate a diffusion constant of 32 cm 2 /s, in agreement with experiment. We ascribe large exciton diffusion lengths to low static and dynamic energetic disorder in LHNs. We argue that matrix-stabilized LHNS represent an excellent model system to study coherent excitonic transport.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Simulation of stochastic diffusion via first exit times
In molecular biology it is of interest to simulate diffusion stochastically. In the mesoscopic model we partition a biological cell into unstructured subvolumes. In each subvolume the number of molecules is recorded at each time step and molecules can jump between neighboring subvolumes to model diffusion. The jump rates can be computed by discretizing the diffusion equation on that unstructured mesh. If the mesh is of poor quality, due to a complicated cell geometry, standard discretization methods can generate negative jump coefficients, which no longer allows the interpretation as the probability to jump between the subvolumes. We propose a method based on the mean first exit time of a molecule from a subvolume, which guarantees positive jump coefficients. Two approaches to exit times, a global and a local one, are presented and tested in simulations on meshes of different quality in two and three dimensions.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Toroidal modeling of plasma flow damping and density pump-out by RMP during ELM mitigation in HL-2A
Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
LA-R43S6-L1 Postshot Report [Slides]
The Ranchero LA-R43S6-L1 test (formerly referred to as LA-43S-L1. The R specifically designates it is a Ranchero device and the 6 indicates a 6” high explosive charge diameter in anticipation of using larger diameters in the future.) was conducted on 3/8/17 with good success. The experiment was the first test of a “Swooped” Ranchero flux compression generator (FCG), and the load for the test was an aluminum imploding liner. In addition to the swoop, the stator was anodized to provide an insulating layer instead of one or more layers of polyethylene as had been the case for all earlier Ranchero tests. Using the Firing point TA-39-88 capacitor bank, 3.5 MA was delivered to the FCG for the initial magnetic field, and 36.4 MA were delivered to the load having an initial inductance of 2.5 nH. Implosion speeds over 1.1 cm/us were recorded, which exceeds previous LANL solid liner experimental results. The test was diagnosed with 18 PDV channels, 12 of which recorded liner implosion data, and the other 6 of which tracked the FCG armature during flux compression. These armature expansion data recorded during flux compression are the first of their kind, and it was previously unknown, due to unknown effects of the SF 6 in the flux compression volume, whether or not the expansion could be recorded using the PDV technique. Post shot 2D MHD calculations employed a complete external circuit model, and an improved flux diffusion model not available prior to the test. Results from these simulations show very good agreement with the experimental result. Four PDV channels recorded the liner implosion viewed radially outward from the center of the liner. All four of these probes recorded implosion velocities greater that 1 cm/µs, with the cylindrical center displaced by ~2 mm from actual center. Four PDV channels looked at +10° angles from the CMU wall, and four looked at -10° angles. A glide plane interaction is shown moving toward the center of the liner in MHD calculations, and PDV probes provide confirmation. The high explosive (HE) in the LA-R43S6-L1 test was PBX 9501 which had to be glued together in many pieces, and the test was preceded by a camera test to verify that the PBX 9501 could be assembled with acceptable tolerance in glue joints to prevent severing the armature during expansion. The camera test verified an acceptable armature performance, and results are given in complete detail in a post shot report LA-UR-19-20124, and summarized here. This report provides complete detail of the considerable body of data obtained on the test and in the post shot calculations. It is prepared in Power Point for ease of preparation and future review. Shot documentation available in the LANL on-line library are included as references, and non-referenceable documents will be cited and stored in an LA- R43S6-L1 post shot folder stored on an M-6 shared drive. Fabrication drawings are also maintained in an M-6 shared drive and paper files are available from M-6 personnel. This report will be maintained as a PowerPoint document on the shared drive, as well, since included movies will play and graphs from Excel files retain information in the PowerPoint versions. PDF versions are required for clearance and are not to large to send by e-mail.
42 ENGINEERING↗
The use of airborne and ground-based measurements of launch vehicle HCl puffs in a finite difference parameter estimation procedure
Model 4 of the multilayer diffusion models, version 5 (MDMV5) is discussed. MDMV5 is designed to help factor considerations of the potential environment impact of a given launch operation into the actual launch decision process. The model supports this decision process by providing quantitative estimates of particular air quality indices at a set of spatial positions in the launch area and surrounding communities, for specified meteorological conditions. Modeling concepts, the finite difference estimation algorithm, and data application and results are examined.
Englar, T. S., Jr.↗