Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fickian diffusion”

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.

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Neutron Reflectometry Reveals Diffusion in Contrast-Matched Brush Particle Bilayers

A material system for performing layer-spread experiments on brush particle bilayers is presented and used to determine the diffusion constant of brush particles in the melt state. Selective deuteration of the core and shell of organo-silica nanoparticles grafted with poly(methyl methacrylate) was used to match the scattering length density of the core and the polymer canopy layer. This subdued the scattering of particle cores (i.e., formfactor scattering) and enabled the analysis of the interdiffusion kinetics using neutron reflectivity. For low molecular grafts, i.e., grafts with a molecular weight below the entanglement limit, the interdiffusion kinetics revealed both a sub- and Fickian diffusion regime. The former was attributed to the local dynamics that was constrained by the slow-moving cores of neighboring brush particles that acted as long-lived physical cross-links. No transition to Fickian diffusion was observed for entangled systems, even at prolonged annealing times. This suggested a higher level of kinetic restraint in entangled brush particle melts as compared to, for example, star polymers with a comparable chain length for which Fickian diffusion has been reported under similar conditions.

Diffusion, Polymer Grafted Nanoparticles, Neutron ↗
NASA NTRS2013

Conductance of Ion Channels - Theory vs. Experiment

Transmembrane ion channels mediate a number of essential physiological processes in a cell ranging from regulating osmotic pressure to transmission of neural signals. Kinetics and selectivity of ion transport is of critical importance to a cell and, not surprisingly, it is a subject of numerous experimental and theoretical studies. In this presentation we will analyze in detail computer simulations of two simple channels from fungi - antiamoebin and trichotoxin. Each of these channels is made of an alpha-helical bundle of small, nongenomically synthesized peptides containing a number of rare amino acids and exhibits strong antimicrobial activity. We will focus on calculating ionic conductance defined as the ratio of ionic current through the channel to applied voltage. From molecular dynamics simulations, conductance can be calculated in at least two ways, each involving different approximations. Specifically, the current, given as the number of charges transferred through the channel per unit of time, can be obtained from the number of events in which ions cross the channel during the simulation. This method works well for large currents (high conductance values and/or applied voltages). If the number of crossing events is small, reliable estimates of current are difficult to achieve. Alternatively, conductance can be estimated assuming that ion transport can be well approximated as diffusion in the external potential given by the free energy profile. Then, the current can be calculated by solving the one-dimensional diffusion equation in this external potential and applied voltage (the generalized Nernst-Planck equation). To do so three ingredients are needed: the free energy profile, the position-dependent diffusion coefficient and the diffusive flux of ions into the channel. All these quantities can be obtained from molecular dynamics simulations. An important advantage of this method is that it can be used equally well to estimating large and small currents. In addition, once the free energy profile becomes available the full current-voltage dependence can be readily obtained. For both channels we carried out calculations using both approaches. We also tested the main assumptions underlying the diffusive model, such as uncorrelated nature of individual crossing events and Fickian diffusion. The accuracy and consistency of different methods will be discussed. Finally we will discuss how comparisons between calculated and measured ionic conductance and selectivity of transport can be used for determining structural models of the channels.

Pohorille, Andrew↗
NASA NTRS1987

Nonlinear diffusion in the solid phase

A Fourier series expansion method has been used to develop simple and relatively accurate solutions for nonlinear Fickian diffusion in the solid phase. For the fundamental harmonic, the first term in the Fourier expansion is qualitatively similar but quantitatively quite different from the conventional exponential function. In contrast, the Fourier coefficients of the higher order harmonics are combinations of fast-decay and slow-decay components. The simple analytical results agree well with numerical solutions for one sample calculation.

Tsang, Tung↗
DOE OSTI2026

Henry’s Solubility and Diffusion Coefficients for 29 Volatile Organic Compounds in Polydimethylsiloxane Sylgard 184 at 293 K

Two-dimensional (2D) inverse gas chromatography (IGC) enables simultaneous determination of Henry’s solubility and Fickian diffusion coefficients for volatile organic compounds (VOCs) in polymer films. This technique offers a significant advantage over traditional cylindrical column IGC by providing precise control and measurement of the film thickness (here, 0.064 ± 0.002 mm), which is the critical length scale for accurate diffusivity determination. We apply this methodology to characterize VOC transport in Sylgard 184, a widely used polydimethylsiloxane (PDMS)-based polymer containing substantial silica filler content. At room temperature (20 °C), we measured solubility and diffusion coefficients for 29 common VOCs spanning diverse chemical functionalities, including alkanes, aromatics, chlorinated solvents, ketones, esters, and alcohols. Comparison with literature data for pure PDMS reveals that VOC solubility in Sylgard 184 is generally higher; for most non-hydrogen-bonding compounds it remains within a factor of 2 of pure PDMS, whereas alcohols are enhanced by roughly 1.8 to 3.7 times, which we attribute to favorable interactions with residual silanol groups on the silanized silica filler. Diffusion coefficients range from 1.0 × 10 –6 cm 2 /s (n-undecane) to 8.9 × 10–5 cm 2 /s (acetonitrile) and align well with extrapolated literature values for PDMS. This study provides essential thermodynamic and transport data for predicting VOC permeation in Sylgard 184 while demonstrating the utility of 2D IGC as a robust technique for characterizing rubbery polymer membranes across diverse industrial applications.

organic↗
DOE OSTI2024

Hydrogen transport in yttrium hydride under asymmetric heat

Metal hydrides are a promising moderator material for high temperature fission reactors. Yttrium hydride can be loaded to a high hydrogen density with relatively high hydrogen stability at temperatures up to 800°C. This makes yttrium hydride a potential moderator material for microreactors as foreseen in the fission surface power program. However, during the operation of such advanced reactors temperature gradients are expected which can change the local hydrogen density in the moderator. Hydrogen diffusion in metals is driven by a concentration gradient (Fick’s law) and thermal diffusion (Soret diffusion). Thermal diffusion is the transport of hydrogen, or other species, due to a temperature gradient. For example, hydrogen might migrate from the hot side of a sample to the cold side of a sample. Measuring Fickian diffusion is achieved through various permeation or absorption experiments, however measuring thermal diffusion is challenging and has rarely been performed. The Hydrogen Experimental Apparatus for Thermal Diffusion (HEATD) experiment is designed to induce thermal diffusion in samples and quench those samples so that the hydrogen distribution can be analyzed using hot vacuum extraction (HVE). One side of the sample was heated to a high temperature e.g., 800°C, while the other side of the sample is at a lower temperature. The sample was held under the applied temperature gradient for a given time until the anticipated hydrogen diffusion has occurred. The actual time depends depend on the sample composition and hydrogen concentration. The results from HVE showed that thermal diffusion took place in the specimen and the Soret coefficient was calculated.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Morphological instability of a thermophoretically growing deposit

The stability of the planar interface of a structureless solid growing from a depositing component dilute in a carrier fluid is studied when the main solute transport mechanism is thermal (Soret) diffusion. A linear stability analysis, carried out in the limit of low growth Peclet number, leads to a dispersion relation which shows that the planar front is unstable either when the thermal diffusion factor of the condensing component is positive and the latent heat release is small or when the thermal diffusion factor is negative and the solid grows over a thermally-insulating substrate. Furthermore, the influence of interfacial energy effects and constitutional supersaturation in the vicinity of the moving interface is analyzed in the limit of very small Schmidt numbers (small solute Fickian diffusion). The analysis is relevant to physical vapor deposition of very massive species on cold surfaces, as in recent experiments of organic solid film growth under microgravity conditions.

Castillo, Jose L.↗
DOE OSTI2026

A New Electrode Paradigm Enabled by Co-Axial Electrospinning

A conventional electrode in a polymer electrolyte membrane fuel cell (PEMFC) often imposes competing requirements on the ionomer, which must enable efficient proton conduction while allowing sufficient oxygen transport. To address these material-level requirements, a co-axial electrospun electrode comprising of a pristine ionomer core and a multifunctional platinum-carbon-ionomer shell is introduced. By decoupling proton and oxygen transport pathways, this architecture enables independent tailoring of core and shell chemistries and their interfaces. The core-shell structure of the electrospun fibers is confirmed by microscopy, while electrochemical and transport properties are quantified using electrochemical impedance spectroscopy and limiting-current diagnostics. The co-axial fibers in their current version meet or exceed the performance of conventional ultrasonically sprayed electrode at relative humidities above 75%, enabled by improvements in non-Fickian diffusion processes. This work establishes co-axial electrospun electrodes as a versatile platform for advanced PEMFC electrodes and a wide range of electrochemical applications, including CO2 electrolysis, water electrolysis, and critical mineral separations.

08 HYDROGEN↗
DOE OSTI2023

Fundamental study of gas species transport in the oxygen electrode of solid oxide fuel and electrolysis cells

A fundamental analysis of multicomponent gas transport models was performed in application to the oxygen electrodes of solid oxide cells. It is common practice to neglect the effect of pressure gradients within oxygen electrodes, even though a net molar flux at the electrolyte surface implies that a pressure gradient must exist. The influence of both Darcy velocity and Knudsen flux are considered in the context of ordinary (Fickian) diffusion, the dusty gas model, and the binary friction model. Comparisons between the models and different sets of assumptions are made via parametric studies on operating load, oxygen partial pressure, microstructural properties, and electrode thickness. Results show that the pressure gradient will have a significant impact on the oxygen concentration distribution and therefore the concentration overpotential. In electrolysis mode, pressure increases up to 1 atm are predicted, indicating that pressure at the electrode/electrolyte interface could contribute to electrode delamination. Additionally, it is found that Darcy's law is insufficient for calculating the pressure distribution without accounting for the flux due to Knudsen diffusion. Additionally, it is found that for the range of properties typical of oxygen electrodes, there is a negligibly small difference between the dusty gas model and binary friction model from a practical standpoint.

08 HYDROGEN↗
NASA NTRS1985

The closure problem for turbulence in meteorology and oceanography

The dependent variables used for computer based meteorological predictions and in plans for oceanographic predictions are wave number and frequency filtered values that retain only scales resolvable by the model. Scales unresolvable by the grid in use become 'turbulence'. Whether or not properly processed data are used for initial values is important, especially for sparce data. Fickian diffusion with a constant eddy diffusion is used as a closure for many of the present models. A physically realistic closure based on more modern turbulence concepts, especially one with a reverse cascade at the right times and places, could help improve predictions.

Pierson, W. J., Jr.↗
DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
DOE OSTI2023

Development and application of sorption mass transfer models for fission product transport in TRISO fuel systems using BISON

Fission product mass transfer within and between TRISO particle and compact layers is an important phenomenon. It directly impacts fission product release predictions, which are used as source terms for safety and licensing calculations. Modeling fission product transport within layers and across bonded interfaces is relatively straightforward under the assumptions of isotropic Fickian diffusion, concentration continuity, and flux continuity. Modeling fission product transport across gaps and debonded layers is more difficult. Gaps are known to form between the buffer and inner PyC (IPyC), and debonding may occur between the IPyC and silicon carbide (SiC). A new mass transfer model was developed to provide a more accurate and robust fission product release calculation by accounting for interlayer sorptivity. One importance of the new model is to account for temperature and material changes across the gap based on sorption isotherm. This work presents the development of the sorption behavioral model and a fission product trapping model, and analyses the subsequent fission product diffusion behavior, specifically cesium, through TRISO particles and compact materials using the BISON fuel performance code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2001

Diffusion Time-Scale of Porous Pressure-Sensitive Paint

Pressure-sensitive paint (PSP) is an optical pressure sensor that utilizes the oxygen quenching of luminescence. PSP measurements in unsteady aerodynamic flows require fast time response of the paint. There are two characteristic time-scales that are related to the time response of PSP. One is the luminescent lifetime representing an intrinsic physical limit for the achievable temporal resolution of PSP. Another is the time-scale of oxygen diffusion across the PSP layer. When the time-scale of oxygen diffusion is much larger than the luminescent lifetime, the time response of PSP is controlled by oxygen diffusion. In a thin homogenous polymer layer where diffusion is Fickian, the oxygen concentration 1021 can be described by the diffusion equation in one-dimension.

Liu, Tianshu↗
NASA NTRS2000

Insert Concepts for the Material Science Research Rack (MSRR-1) of the Material Science Research Facility (MSRF) on the International Space Station (ISS)

The Material Science Research Rack I (MSRR-1) of the Material Science Research Facility (MSRF) contains an Experiment Module (EM) being developed collaboratively by NASA and the European Space Agency (ESA). This NASA/ESA EM will accommodate several different removable and replaceable Module Inserts (MIs) which are installed on orbit NASA's planned inserts include the Quench Module Insert (QMI) and the Diffusion Module Insert (DMI). The QMI is a high-gradient Bridgman-type vacuum furnace with quench capabilities used for experiments on directional solidification of metal alloys. The DMI is a vacuum Bridgman-Stockbarger-type furnace for experiments on Fickian and Soret diffusion in liquids. This paper discusses specific design features and performance capabilities of each insert. The paper also presents current prototype QMI hardware analysis and testing activities and selected results.

Crouch, Myscha↗
DOE OSTI2023

Interpreting inorganic compositional depth profiles to understand the rate-limiting step in vapor phase infiltration processes

Vapor phase infiltration (VPI) is a post-polymerization modification technique that infuses inorganics into polymers to create organic–inorganic hybrid materials with new properties. Much is yet to be understood about the chemical kinetics underlying the VPI process. The aim of this study is to create a greater understanding of the process kinetics that govern the infiltration of trimethyl aluminum (TMA) and TiCl 4 into PMMA to form inorganic-PMMA hybrid materials. To gain insight, this paper initially examines the predicted results for the spatiotemporal concentrations of inorganics computed from a recently posited reaction–diffusion model for VPI. This model provides insight on how the Damköhler number (reaction versus diffusion rates) and non-Fickian diffusional processes (hindering) that result from the material transforming from a polymer to a hybrid can affect the evolution of inorganic concentration depth profiles with time. Subsequently, experimental XPS depth profiles are collected for TMA and TiCl 4 infiltrated PMMA films at 90 °C and 135 °C. The functional behavior of these depth profiles at varying infiltration times are qualitatively compared to various computed predictions and conclusions are drawn about the mechanisms of each of these processes. TMA infiltration into PMMA appears to transition from a diffusion-limited process at low temperatures (90 °C) to a reaction-limited process at high temperatures (135 °C) for the film thicknesses investigated here (200 nm). While TMA appears to fully infiltrate these 200 nm PMMA films within a few hours, TiCl 4 infiltration into PMMA is considerably slower, with full saturation not occurring even after 2 days of precursor exposure. Infiltration at 90 °C is so slow that no clear conclusions about mechanism can be drawn; however, at 135 °C, the TiCl 4 infiltration into PMMA is clearly a reaction-limited process, with TiCl 4 permeating the entire thickness (at low concentrations) within only a few minutes, but inorganic loading continuously increasing in a uniform manner over a course of 2 days. Near-surface deviations from the uniform-loading expected for a reaction-limited process also suggest that diffusional hindering is high for TiCl 4 infiltration into PMMA. In conclusion, these results demonstrate a new, ex situ analysis approach for investigating the rate-limiting process mechanisms for vapor phase infiltration.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Flux-difference split parabolized Navier-Stokes algorithm for non-equilibrium chemically reacting flows

A flux-difference split explicit finite-difference algorithm is presented for solving the parabolized form of the equations governing three-dimensional nonequilibrium chemically reacting flows. The algorithm is based on an explicit noniterative, upwind space-marching scheme developed by Korte, but differs in that the unsteady Riemann problem, rather than the steady Riemann problem, is solved. The algorithm allows either a second or an approximately third-order accurate upwind treatment of the convection terms by employing the unsteady approximate Riemann solver of Roe. The source terms of the species transport equations are treated in either an explicit or implicit manner, and the species diffusion terms are modeled with either a Fickian or a multicomponent model. A validation of the algorithm is performed by comparing computational results with the 2-D Mach 14, 15 degree compression-corner data of Holden. The three-dimensional capability of the algorithm is demonstrated by computing Mach 2.7 flow over a swept wedge scramjet fuel injector, and three-dimensional reacting flow capability is demonstrated by a computing a shock-jet interaction concept for mixing and combustion enhancement.

White, J. A.↗