SEARCH · Engineering Papers
Results for “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.
Modeling Farmers’ Adoption Potential to New Bioenergy Crops: An Agent-Based Approach
The use of fossil fuels is the primary source of greenhouse gas emissions but there are alternatives to these especially in the form of biofuels, fuels derived from bioenergy crops. This paper aims to determine farmers’ potential adoption rates of newly introduced bioenergy crops with a specific example of carinata in the state of Georgia. The determination is done using an agent-based modeling technique with two principal assumptions—farmers are profit maximizer and they are influenced by neighboring farmers. Two diffusion parameters (traditional and expansion) are followed along with two willingness (high and low) scenarios to switch at varying production economics to carinata and other prominent traditional field crops (cotton, peanuts, corn) in the study region. We find that a contract prices around $9, $8 and $7 can be a viable option for encouraging farmers to adopt carinata in low, average, and high profit conditions, respectively. Expansion diffusion (that diffuses all over the geographical area), rather than centered to the few places like traditional diffusion at the early stage of adoption in conjunction with higher willingness conditions influences higher adoption rates in the short-term. As such, the model can be used to understand the behavioral economics of carinata in Georgia and beyond, as well as offering a potential tool to study similar bioenergy crops.
Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters
Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.
Grain-boundary structure and segregation in Nb 3 Sn coatings on Nb for high-performance superconducting radiofrequency cavity applications
Here, we report on atomic-scale analyses of grain boundary (GB) structures and segregation in Nb 3 Sn coatings on Nb, prepared by the vapor-diffusion process, for superconducting radiofrequency (SRF) cavity applications, utilizing atom-probe tomography, high-resolution scanning transmission electron-microscopy and first-principles calculations. We demonstrate that the chemical composition of Nb 3 Sn GBs is correlated strongly with the diffusion of Sn and Nb at GBs during the coating process. In a sample coated with a relatively large Sn flux, we observe an interfacial width of Sn segregation at a GB of ~3 nm, with a maximum concentration of ~35 at.%. After post-annealing at 1100 °C for 3 h, the Sn segregated at GBs disappears and Nb segregation is observed subsequently at GBs, indicating that Nb diffused into the Nb3Sn GBs from the Nb substrate. It is also demonstrated that the amount of Sn segregation in a Nb 3 Sn coating can be controlled by: (i) Sn flux; and (ii) the temperatures of the Nb substrates and Sn source, which may affect the overall kinetics including GB diffusion of Sn and Nb. An investigation of the correlation between the chemical compositions of GBs and Nb 3 Sn SRF cavity performance reveals that the Nb 3 Sn SRF cavities with the best performance (high-quality factors at high accelerating electric-field gradients) do not exhibit Sn segregation at GBs. Our results suggest that the chemical compositions of GBs in Nb 3 Sn coatings for SRF cavities can be controlled by GB engineering and can be utilized to optimize fabrication of high-quality Nb 3 Sn coatings for SRF cavities.
Permeation of niobium through grain boundaries in copper
Low mutual solubility is expected to avert interdiffusion in layered composites of phase-separating metals. However, we show that Nb diffuses through polycrystalline Cu—despite the minimal bulk solubility of these elements—due to short circuit transport of Nb along certain grain boundaries in Cu. Atomistic modeling demonstrates that Nb-permeable Cu grain boundaries exhibit negative enthalpy of mixing of Nb, resulting in an enthalpically-stabilized solution of highly mobile Nb atoms that easily diffuse through the boundary. By contrast, Mo, which also has minimal solubility with Cu, has positive enthalpy of mixing at Nb-permeable Cu grain boundaries and does not diffuse through them. Furthermore, our findings suggest material selection and grain boundary engineering as pathways for designing improved diffusion barriers and thermally stable laminate composites.
Permeation of niobium through grain boundaries in copper
Low mutual solubility is expected to avert interdiffusion in layered composites of phase-separating metals. However, we show that Nb diffuses through polycrystalline Cu—despite the minimal bulk solubility of these elements—due to short circuit transport of Nb along certain grain boundaries in Cu. Atomistic modeling demonstrates that Nb-permeable Cu grain boundaries exhibit negative enthalpy of mixing of Nb, resulting in an enthalpically-stabilized solution of highly mobile Nb atoms that easily diffuse through the boundary. By contrast, Mo, which also has minimal solubility with Cu, has positive enthalpy of mixing at Nb-permeable Cu grain boundaries and does not diffuse through them. Furthermore, our findings suggest material selection and grain boundary engineering as pathways for designing improved diffusion barriers and thermally stable laminate composites.
A chain stretch-based gradient-enhanced model for damage and fracture in elastomers
Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.
Effects of various parameters of different porous transport layers in proton exchange membrane water electrolysis
Porous transport layers (PTLs) play an important role in proton exchange membrane water electrolysis (PEMWE) cells. The PTL facilitates water and gas transport, as well as thermal and electrical conduction, and is required to sustain good contact with adjacent components. It is expected that using PTLs with variations in material properties such as structure, composition, thickness and wettability results in performance changes of the PEMWE. Here, a general mathematical PEMWE model is developed that separates and analyzes the contributions of ohmic, activation, diffusion and Nernst potentials. For model validation, three inherently different anode PTL structures (carbon paper, sintered titanium particles, and titanium felt) are operated over a range of conditions. Additionally, the effects of PTL wettability were used to verify the model using Polytetrafluoroethylene (PTFE) treated Toray papers with PTFE loading ranging from 0% to 20%. The modeling results of both PTFE treated and untreated materials show good agreement with the experimental data. Mass transport or diffusion loss is the primary reason for performance differences between PTFE treated and untreated PTLs. Sintered titanium PTLs with thicknesses above 1 mm suffer from up to 33% increased ohmic losses without indicating any obvious changes in activation and diffusion losses when compared to untreated PTLs. The losses of the cell increase when using PTFE treated Toray paper. Individual contributions are quantified and assigned to increased ohmic, activation, and diffusion losses. In conclusion, the proposed model offers insights into the overpotential contributions of a PEMWE. It is a useful tool for predicting performance of various PTL materials and can be applied for PTL development and optimization efforts.
Methodologies for Design, Characterization and Testing of Electrolytes that Enable Extreme Fast Charging of Lithium-ion Cells
Selection, testing and validation of electrolyte candidates for Li-ion cells are discussed, based on a 10-minute target for extreme fast charge (XFC). A combination of modeling and laboratory measurements create a timely and synergistic approach to identifying candidate electrolyte formulations. Multi-solvent systems provide a balanced set of properties, wherein lower molecular-weight solvents offer reduced viscosity, increased species diffusivity, and mitigation of concentration polarization at high charge rates. Carefully selected formulations can exhibit peak conductivity and usable conductivity range of two to three times that of the baseline EC-EMC (3:7, wt.) + LiPF 6 . Candidates are also chosen based on stability and longevity within the cell environment. Lab testing coincides with property predictions from the Advanced Electrolyte Model (AEM) and a macro-scale cell model. Furthermore, cell testing utilized coin and pouch cells having NMC532 or NMC811 cathodes with graphite electrodes. Results indicate combinations of low-molecular weight solvents are key for fast-charge electrolytes as they extend the useful conductivity range to both low and higher salt concentrations, and possess higher self-diffusivities compared to conventional solvents. This reduces impacts from concentration polarization. The choice of electrolyte also influences the tendency for lithium metal deposition at the anode, as showcased by experimental and modeling results herein.
The effects of turbulence on the flame structure and NO formation of ammonia turbulent premixed combustion at various equivalence ratios
Ammonia is carbon-free and is regarded as a potential fuel to address global warming issues. Here, in this work, three-dimensional direct numerical simulations (DNS) of ammonia/air turbulent premixed flames were performed to explore the influence of turbulence and equivalence ratio on the flame structure and NO formation characteristics. Two equivalence ratios were considered, i.e. Φ=0.9 and Φ=1.1. The general flame structures were presented and species distributions were examined. The NO mass fraction was found to be the highest in the product for the lean case and in the reaction zone for the rich case. The conditional mean values of species mass fractions and reaction rates were compared with those of the unstrained and strained laminar premixed flames to explore how well the laminar flame structures can approach those of the turbulent flame. The budget analysis of the species transport equations showed that turbulent diffusion plays an important role in species transport. The turbulent diffusivity D T was estimated using the gradient hypothesis based on the DNS data. Various laminar flame simulations with different diffusivities were carried out. It was shown that the conditional means of the DNS agree well with those of the laminar flames including D T in the transport property calculation. The global and local NO formation characteristics were investigated. It was found that the mean NO production conditioned on the progress variable is lower compared with the corresponding laminar flame in the rich case. However, the relative contributions from various NO pathways are rarely affected by turbulence. The NO mass fraction is higher in negative curvature regions compared with positive curvature regions of the flame surface for the rich case, which is due to the preferential diffusion of H 2 and other radicals and the enhanced NO pathways in negatively curved regions.
Achieving clean production with nanostructured coated milling tools dry machining low carbon steel
The advancement of clean production using nanostructured materials in subtractive manufacturing processes has focused on replacing the use of liquid lubrication with solid lubricants coated to conventional cutting tools machining low carbon steel (~ 0.2 wt. % carbon). However, little is known about the wear mechanisms that dominate such tools when coated with functional graded nanostructured coatings. The present study follows an international standard for characterizing such wear (ISO 8688) but uses advanced measurement techniques such as x-ray fluorescence technology to measure the diffusion of chemical species from tool to workpiece and vice versa, that contributes to tool wear with the view to enhance the principles of cleaner production. The results show that when commercially-available tools coated with the appropriate functionally graded industrial coating, diffusion does not take place and creates an effective chemical barrier during machining. The experimental methods used in the study include machining workpieces using a CNC milling machine, tachometer for measuring spindle speed, dynamometer for measuring cutting forces, infra-red camera for measuring cutting temperatures, chemical species diffusion wear using x-ray florescence detector, and an optical microscope for measuring tool wear. Coupled with the physical measurement of wear, it is concluded that mechanical wear dominates the milling of a low carbon steel and that thermal properties of the coatings do not correlate with such wear. It is also concluded that flank wear as a function of volume of workpiece removed and machining duration (machining index) is a sound method for assessing the transition of the stages of wear during clean machining operations. The present study not only has implication for the design of better cutting tools, but also advances dry machining processes that eliminate the use of liquid lubricants making subtractive production processes cleaner. The main conclusions drawn from this study show that coated tools have a longer life under dry conditions compared to uncoated tools and that there is no correlation between coating properties and tool wear. In conclusion, diffusion of chemical species into the cut chips did not occur because of the lack of thermally induced chemical wear of the cutting tool and that wear is caused by the gradual erosion of the flank face due to the chip abrading the cutting tool.
Physics-informed machine learning with conditional Karhunen-Loève expansions
Not Available
Behavior of B in passivating gels formed on International Simple Glass in acid and basic pH
International Simple Glass altered for years in silica saturated conditions develops a passivating gel, which retains some mobile elements such as B and Ca in its inner region. It has been suggested that the retention of these elements impacts the dissolution rate of the glass. However, the elements' retention mechanism, speciation, and diffusion properties are still unknown. This study finds that simultaneous presence of B and Ca in the solution can dramatically slow down glass alteration. However, splitting the B and Ca into two different solutions with other glass elements did not give similar protection against glass alteration. These experiments along with other experiments conducted in basic or acid pH, in which isotopically tagged B can diffuse in gels revealed that (i) water diffusion is not dramatically affected by the retention of B and Ca in passivating gels, and (ii) slow diffusion and high retention of B is observed only in gels developed in basic pH during the tracing experiments. Finally, this last mechanism therefore appears as one of the many mechanisms controlling the glass dissolution rate in the residual regime.
Hydrogen and deuterium permeation in Hastelloy N
Hastelloy N was chosen as the fluoride salt-contacting structural material for the Molten Salt Reactor Experiment due to its excellent compatibility with the fuel salt FLiBe. FLiBe is currently investigated for several advanced fusion and fission reactor concepts where tritium generation in the FLiBe is anticipated. Knowledge of hydrogen transport properties through Hastelloy N is important to understand how tritium would permeate through this material and result in an unintentional release. In this study, the hydrogen and deuterium permeability, diffusivity, and solubility were measured from 500 to 700 °C at a primary-side pressure of 10 kPa in a well-characterized sample of Hastelloy N. The prepared polycrystalline Hastelloy N had C and O impurities present on the surface. These impurities were investigated using Auger Emission Spectroscopy and Ar depth profiling. Here, the adventitious C was removed upon the first Ar sputter cycle whereas O persisted deeper into the sample. For permeation experiments, applied deuterium pressures ranged from 13 Pa to 100 kPa and deuterium transport remained in the diffusion-limited regime (J ∝ P 0.5 ) throughout the pressure range examined. Two methods are employed to measure the effective hydrogen and deuterium diffusivity: rise and decline. The decline method produced improved statistical model fits for calculating the effective diffusion coefficient compared to the rise method. The resultant transport properties compared well to published values for other nickel alloys.
A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON
TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.
The influence of preexisting host rock damage on fault network localization
The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.
Atomic scale understanding of organic anion separations using ion-exchange resins
A combination of ab initio and classical molecular dynamic simulations was used to explore the adsorption/desorption and diffusion characteristics of ion-exchange resins for extraction of organic anions in electrically-driven separation processes. We considered two classes of carboxylate mixtures that are commonly encountered in bioprocessing separations: a short chain fatty acid mixture (acetate/butyrate) and an aromatic mixture (ferulate/coumarate). The suitability of several resin materials including PFC100E, IRC86, PFA444 and IRA67 was interrogated. The decomposition of the interaction energies by the symmetry-adapted perturbation theory and the classical molecular dynamic simulations of organic diffusion together reveal that the geometries of the organic anions and the functional groups of the resins, as well as their Columbic interactions, are the controlling factors in the diffusion process of the organic compounds in these resins. Additionally, classical simulations also show that modifying the functionality of resin beads, the magnitude of electric fields, and the ratio of organic mixtures may provide an effective way to control the diffusion rate and obtain selective separation of these organic mixtures. Finally, a general suggestion for favorable separation conditions is summarized.
Framework for predicting the fractionation of complex liquid feeds via polymer membranes
The separation of complex liquid hydrocarbon mixtures was recently demonstrated using the glassy polymer SBAD-1, showing that small molecule fractionation is possible by such organic membrane materials. Here, in this work, we develop a framework that will enable workable predictions of permeate flux and composition in complex hydrocarbon liquids through intrinsically porous glassy polymers. The predictions are made by incorporating experimentally-derived unary sorption and diffusion parameters in a Maxwell-Stefan framework coupled with multicomponent sorption models and various distinct diffusion phenomena. Across the range of sorption and diffusion phenomena considered, both the conventional Flory-Huggins model and the proposed Langmuir + Flory-Huggins sorption model combined with a simple average guest diffusivity or a more complex free-volume theory-based transport resulted in the lowest prediction error for three chosen multicomponent separations. The proposed Maxwell-Stefan framework simply requires pure component transport parameters to allow a fast approximation of the separation of multicomponent liquid hydrocarbon feeds that can potentially be extended to more complex feeds such as crude oil fractions.