MXE: A LAMMPS package for simulating long-term diffusive mass transport in nanoscale systems
Explore the source record for details and available documents.
Mendez Granado, Juan Pedro↗
SEARCH · Engineering Papers
Results for “mass diffusivity”
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.
Representing effects of surface heterogeneity in a multi-plume eddy diffusivity mass flux boundary layer parameterization
Explore the source record for details and available documents.
54 ENVIRONMENTAL SCIENCES↗
Effects of thermophoresis on high-pressure binary-species boundary layers with uniform and non-uniform compositions
Direct numerical simulations of high-pressure binary-species temporal boundary layers are performed to investigate the flow physics for three situations: (1) uniform and equal composition, (2) uniform but unequal compositions and (3) non-uniform composition. Both colder- and hotter-wall situations compared with the free stream are simulated. The working fluid is a nitrogen/methane mixture. The analysis is performed at a case-specific self-similar state. Even when the initial composition is uniform, the methane mean mass fraction decreases near the colder wall, whereas it increases near the hotter wall and the mass fraction fluctuates in the entire boundary layer. Analysis of the species-mass diffusion balance and flow structures reveal that both mass-fraction variation and fluctuations are induced by the Soret effect. When the initial composition is non-uniform and the wall is colder, the methane mean mass fraction monotonically increases from the wall akin to its initial profile. However, when the wall is hotter the mean mass fraction decreases near the wall in contrast to its initial profile, a fact traced through the species-mass diffusion balance to the Soret effect being large and enriching methane near the wall. In contrast, the direction of the Soret flux is opposite for the colder wall, thus keeping the methane concentration small. Although the initial magnitude of the difference between the wall and free-stream temperature is the same in all cases, the situation is not symmetric between colder- and hotter-wall cases; the flow structure exhibits much smaller scales when the wall is hotter than when the wall is colder.
Mechanics↗
Measurement and modeling of methane diffusion in hydrocarbon mixtures
Methane (CH 4 ) dissolution and diffusive mass transfer in liquid hydrocarbon mixtures is of key interest in the context of enhanced oil recovery from tight shales. In this paper, we have studied CH 4 dissolution and diffusion in normal alkane mixtures at a temperature of 50 ⁰C and at pressures of ~8 MPa. For the measurement of CH 4 dissolution/diffusion in bulk liquid hydrocarbon mixtures, we have utilized a high pressure and temperature, constant-volume diffusion (CV-D) setup. Here, we have studied CH 4 diffusion in mixtures with three long-chain normal alkanes: decane (C 10 ), dodecane (C 12 ) and hexadecane (C 16 ). We have measured CH 4 solubility and diffusion in its binary mixtures with each of these three normal alkanes, as well as in ternary and quaternary mixtures. During the experiments, the swelling of the liquid mixture due to CH 4 dissolution was measured in situ via a cathetometer and was subsequently integrated into the data analysis. A key conclusion from this study is that the solubility and transport properties of the multicomponent mixtures can be predicted accurately from binary mixture measurements using an appropriate Equation of State (EOS) and Wilke’s simplification of the classical Maxwell-Stefan (MS) diffusion theory. This observation can facilitate accurate prediction of diffusive mass transfer in more complex liquid hydrocarbon mixtures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A general synthesis of single atom catalysts with controllable atomic and mesoporous structures
The control of single metal atomic sites has been extensively studied in the field of single atom catalysts. By contrast, the precise control of the mesoporous structure in the matrix material, which directly correlates with mass diffusions and may play a dominant role in delivering industrially relevant reaction rates, has been overlooked. In this work, we report a general method for the synthesis of a single atom catalyst with control of the atomic structure of the single atomic site as well as the mesoporous structure of the carbon support for optimized catalytic performance. Various combinations of metal centres (Ni, Co, Mn, Zn, Cu, Sc and Fe) and mass diffusion channels in two dimensions and three dimensions were achieved. Using CO 2 reduction to CO as an example, our Ni single atom catalyst with three-dimensional diffusion channels delivered a practical current of 350 mA cm –2 while maintaining a 93% CO Faradaic efficiency, representing a sixfold improvement in turnover frequency compared to two-dimensional counterparts.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion
Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.
Science & Technology - Other Topics↗
Volumetric lattice Boltzmann method for pore-scale mass diffusion-advection process in geopolymer porous structures
Not Available
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Diffusion-based mass map reconstruction from weak lensing data
Not Available
Boruah, Supranta S. (ORCID:0000000309851488)↗
Verification of Bison fission product species conservation under TRISO reactor conditions
When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.
07 ISOTOPE AND RADIATION SOURCES↗
The effect of saturated thermal conduction on clouds in a hot plasma
We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.
79 ASTRONOMY AND ASTROPHYSICS↗
Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit
Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.
97 MATHEMATICS AND COMPUTING↗
Irradiation-enhanced Interactions at UO 2 /Al 2 O 3 /Al Interfaces
Combustion synthesis is used to prepare thin UO 2 films on aluminum alloy substrates. This simple preparation method involves electrospraying uranyl nitrate + acetylacetone + 2-methoxyethanol solution on the substrate, followed by a short annealing at 350 or 550 °C. The irradiation of films with a 40 Ar 2+ ion beam (energy of 1.7 MeV and fluences of 7.7 × 10 16 and 1.3 × 10 17 ions/cm 2 ) is conducted to investigate irradiation-induced restructuring processes. High-resolution transmission electron microscopy (TEM) and X-ray photoelectron spectroscopy (XPS) investigations show that the annealing temperature significantly influences the crystallinity and stability of materials during irradiation. A small amount of Mg in the alloy substrate diffuses into the amorphous Al 2 O 3 interfacial layer between the film and the substrate. Local thermal spikes from the incoming ions facilitate the irradiation-induced mixing of immiscible Al 2 O 3 and UO 2 for the materials prepared at 350 °C. This mass diffusion produces relatively large cavities at the interface. Selective diffusion of a more significant amount of Mg for the materials prepared at 550 °C suppresses the mixing of the Al 2 O 3 interlayer with the film but forms Mg y U 1–y O 2±x solid solutions during irradiation. Local thermal heating triggers the precipitation of a discontinuous crystalline MgO layer close to the film surface. As a result, the enhanced and selective diffusion of Mg into the film makes the materials prepared at 550 °C more robust and mechanically stable than those prepared at 350 °C.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Time-Varying Model for Predicting Formaldehyde Emission Rates in Homes
Recent studies have succeeded in relating emissions of various volatile organic compounds to material mass diffusion transfer using detailed empirical characteristics of each of the individual emitting materials. While significant, the resulting models are often scenario specific and/or require a host of individual component parameters to estimate emission rates. This study developed an approach to estimate aggregated emissions rates based on a wide number of field measurements. We used a multi-parameter regression model based on previous mass transfer models to predict formaldehyde emission rate for a whole dwelling using field-measured, time-resolved formaldehyde concentrations, air exchange rates, and indoor environmental parameters in 63 California single-family houses built between 2011 and 2017. The resulting model provides time-varying formaldehyde emission rates, normalized by floor area, for each study home, assuming a well-mixed mass balance transport model of the home, and a well-mixed layer transport model of indoor surfaces. The surface layer model asserts an equilibrium concentration within the surface layer of the emitted materials that is a function of temperature and RH; the dwelling ventilation rate serves as a surrogate for indoor concentration. We also developed a more generic emission model that is suitable for broad prediction of emission for a population of buildings. This model is also based on measurements aggregated from 27 homes from the same study. We showed that errors in predicting household formaldehyde concentrations using this approach were substantially less than those using a traditional constant emission rate model, despite requiring less unique building information.
54 ENVIRONMENTAL SCIENCES↗
Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field
Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.
09 BIOMASS FUELS↗
Matrimid substrates with bicontinuous surface and macrovoids in the bulk: A nearly ideal substrate for composite membranes in CO 2 capture
Scalable fabrication of thin-film composite (TFC) membrane for post-combustion carbon capture is often limited by the availability of high-performance nanoporous substrate. Ideally, the substrate should allow for fast gas transport at the selective layer/substrate interface as well as in the bulk of the substrate. In this study, highly permeable Matrimid substrates were prepared via vapor- and nonsolvent-induced phase separations. The phase separation mechanism of the Matrimid®/N-methyl-2-pyrrolidone casting solution was modulated by the addition of LiCl to control the solution thermodynamic stability, and nanoporous Matrimid substrates were formed with a bicontinuous surface and macrovoids in the bulk. This hierarchically optimized structure resulted in a CO 2 permeance of 2.60 × 10 5 GPU (1 GPU = 1 × 10 –6 cm 3 (STP) cm –2 s –1 cmHg –1 = 3.349 × 10 –10 mol m –2 s –1 Pa –1 ), which was ca. 11 times more permeable than a benchmark polyethersulfone substrate with cellular pores. The benefit of using this new substrate was demonstrated by coating a 170-nm amine-containing polymer to form a TFC facilitated transport membrane. The membrane exhibited a CO 2 permeance of 932 GPU at 57°C, which was 72 GPU higher than the counterpart coated on the benchmark substrate. Meanwhile, the CO 2 /N 2 selectivity was remained at 158. This permeance improvement could be well explained by the resistance-in-series model, where the improved permeance was attributed to the reductions in substrate and lateral diffusion mass transfer resistances. The upper bound (UB) analysis indicates that the Matrimid substrate is nearly ideal for the state-of-the-art polymers for CO 2 /N 2 separation. Here, the substrate improvement could also reduce the parasitic energy associated with the membrane process due to the reduced membrane cost and flue gas compression requirement.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Neutron capture of UO 2 targets prepared by spin-coating assisted combustion synthesis
Two uranium dioxide (UO 2 ) targets of (414 ± 23) nm and (1092 ± 93) nm thicknesses were prepared on 6061 aluminum alloy and puratronic grade aluminum backing materials. The targets were deposited with a novel method combining spin coating and solution combustion synthesis (SCS). The target layers consisted of small (3–7 nm) UO 2 grains and uniformly distributed ultra-small (1–3 nm) pores. The prepared targets were tested at the Los Alamos National Laboratory’s LANSCE facility for neutron irradiation damage and suitability for neutron capture experiments. The samples showed no signs of target material loss after the irradiation. However, irradiation caused a significant increase in the grain size (4–10 nm), as well as upward mass diffusion and coalescence of the pores due to the thermal spikes. The magnesium in the aluminum 6061 alloy backing also diffused into the UO 2 layer during neutron irradiation. The structural changes in the target after the irradiation do not affect the data from neutron capture. As a result, the new method can be used more broadly to prepare other actinide targets for nuclear physics experiments.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Versatile TRISO fuel particle modeling in Bison
Tri-structural isotropic (TRISO) fuel particles are a key component in several previous and current reactors as well as in a variety of novel nuclear reactor designs. Interest in TRISO fuel is on the rise, necessitating considerable computer modeling of TRISO fuel behavior in order to support related design and licensing activities. The Bison nuclear fuel performance code, which offers a full set of capabilities for modeling TRISO fuels, makes it easier to explore the various important aspects of TRISO fuel behavior. One key advantage of Bison is its ability to create meshes in 1D, 2D, and 3D. Users can customize these meshes for specific geometries, mesh densities, and use cases. This enables a wide variety of analyses, including thermal, structural, mass diffusion, homogenization, and statistical failure analyses. Furthermore, the meshing capability simplifies analysts’ workflows. The inherent mesh generation capability eliminates the need for separate mesh-generating software and mesh file management. Also, the fact that the meshes are customizable makes it straightforward to automate an investigation over a range of geometric parameters or mesh densities. Here, the present paper highlights the ease with which Bison may be used to create meshes for both simple and relatively complex TRISO fuel particles, and it explores the types of analyses enabled by these meshes.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage
Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.
08 HYDROGEN↗