Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling

We have treated the modeling, analysis, numerical analysis, and algorithmic development for nonlocal models of diffusion and mechanics. Variational formulations were developed and finite element methods were developed based on those formulations for both steady state and time dependent problems. Obstacle problems and optimization problems for the nonlocal models were also treated and connections made with fractional derivative models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

A sub-pore model for multi-scale reaction–diffusion problems in porous media

Applications of reaction–diffusion systems in porous media pose a challenging problem for computational modeling approaches due to their multi-physics and multi-scale nature. The length scales usually span 3–4 orders of magnitude while physical phenomena involved include heat and mass transfer processes, and surface reactions. In this paper, a novel methodology that accounts for all the length scales and physical phenomena involved in a single framework is described. A length scale based dual approach is proposed – the larger pore channels (macro-pores) are resolved using conventional numerical techniques and a novel ‘sub-pore’ model is used to account for the unresolved pore channels (sub-pores) and the important physics therein. The porous network in the sub-pore system is composed of a fractal-like hierarchical system of straight cylindrical pores. Simplified governing equations for mass and energy transport are solved within the sub-pore system along with a reaction kinetics model to account for surface adsorption. An implicit coupling strategy is used to couple the macro-pore and the sub-pore systems so as to ensure conservation. The developed methodology is then applied to a few test cases and it is established that the proposed framework is necessary for problems where the adsorption time scale is much smaller than (diffusion-limited) or comparable to the diffusion time scale. Furthermore, it is also demonstrated that the framework can be potentially used to model the network of porous channels in its entirety thus significantly reducing computational costs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Efficient Reformulation of Linear and Nonlinear Solid-Phase Diffusion in Lithium-ion Battery Models using Symmetric Polynomials: Mass Conservation and Computational Efficiency

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms, along with suitable discretization or approximations for the solid-phase diffusion equation. The solid-phase diffusion equation represents the main computational burden for typical pseudo-2-dimensional (p2D) models since these equations in the pseudo r -dimension must be solved at each point in the computational grid. This substantially increases the complexity of the model as well as the computational time. Traditional approaches towards simplifying solid-phase diffusion possess certain significant limitations, especially in modeling emerging electrode materials which involve phase changes and variable diffusivities. A computationally efficient representation for solid-phase diffusion is discussed in this paper based on symmetric polynomials using Orthogonal Collocation and Galerkin formulation (weak form). A systematic approach is provided to increase the accuracy of the approximation (p form in finite element methods) to enable efficient simulation with a minimal number of semi-discretized equations, ensuring mass conservation even for non-linear diffusion problems involving variable diffusivities. These methods are then demonstrated by incorporation into the full p2D model, illustrating their advantages in simulating high C-rates and short-time dynamic operation of Lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

A Method for Obtaining Liquid–Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H 2 O

Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas–solid adsorption in both theory and experiment, such rates are still not well studied for liquid–solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid–solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid–solid adsorption systems.

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

Ab Initio Design of High-Entropy Thermal/ Environmental Barrier Coatings

Next generation thermal/environmental barrier coatings (TEBC) require carefully balancing various properties including phase stability, thermal conductivity, coefficient of thermal expansion (CTE), mechanical properties, and resistance against hot corrosion and water vapor recession. This work mainly focuses on rapid design of cost-effective high entropy rare-earth disilicates and aluminum garnets to protect SiC-based ceramic matrix composites and nickel-based superalloys in the hot section of gas turbine engines using density functional theory methods. Our calculations identify several low-cost high entropy TEBC exhibiting ultralow thermal conductivity at 1500 K and desirable CTE while maintaining good mechanical properties, including Er1/2Y3/4Yb3/4Si2O7, Gd1/4Er1/4Y3/4Yb3/4Si2O7, Eu1/4Er1/4Y3/4Yb3/4Si2O7, and (Y1/4Gd1/4Er1/4Yb1/4)3Al5O12. This work also aims to gain fundamental understanding of oxygen diffusion in model disilicates. Minimizing oxidizer (such as water vapor and oxygen) permeability through the EBC layer can significantly decrease the growth rate of thermally grown oxide and extend the service life of the coating system. Oxygen diffusion mechanisms including formation energy of defects under varying oxygen conditions and defect migration energy barriers will be presented.

coefficient of thermal expansion↗
DOE OSTI2019

Studying Nanosecond Discharge in Argon at Atmospheric Pressure with Preionization

The influence of initial conditions on specific features of the formation and development of a cathode-directed ionization wave between two flat electrodes in argon at atmospheric pressure was studied at nanosecond time resolution using a high-speed photoelectron monitor and numerical simulations based on the two-dimensional axisymmetric diffusion–drift model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2010

LHC Beam Diffusion Dependence on RF Noise: Models And Measurements

Radio Frequency (RF) accelerating system noise and non-idealities can have detrimental impact on the LHC performance through longitudinal motion and longitudinal emittance growth. A theoretical formalism has been developed to relate the beam and RF loop dynamics with the bunch length growth [1]. Measurements were conducted at LHC to validate the formalism, determine the performance limiting RF components, and provide the foundation for beam diffusion estimates for higher energies and intensities. A brief summary of these results is presented in this work. During a long store, the relation between the energy lost to synchrotron radiation and the noise injected to the beam by the RF accelerating voltage determines the growth of the bunch energy spread and longitudinal emittance. Since the proton synchrotron radiation in the LHC is very low, the beam diffusion is extremely sensitive to RF perturbations. The theoretical formalism presented in [1], suggests that the noise experienced by the beam depends on the cavity phase noise power spectrum, filtered by the beam transfer function, and aliased due to the periodic sampling of the accelerating voltage signal V{sub c}. Additionally, the dependence of the RF accelerating cavity noise spectrum on the Low Level RF (LLRF) configurations has been predicted using time-domain simulations and models [2]. In this work, initial measurements at the LHC supporting the above theoretical formalism and simulation predictions are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Printable organic thin film transistors for glucose detection incorporating inkjet-printing of the enzyme recognition element

The effect of device architecture upon the response of printable enzymatic glucose sensors based on poly(3-hexythiophene) (P3HT) organic thin film transistors is presented. The change in drain current is used as the basis for glucose detection and we show that significant improvements in drain current response time can be achieved by modifying the design of the sensor structure. In particular, we show that eliminating the dielectric layer and reducing the thickness of the active layer reduce the device response time considerably. The results are in good agreement with a diffusion based model of device operation, where an initial rapid dedoping process is followed by a slower doping of the P3HT layer from protons that are enzymatically generated by glucose oxidase (GOX) at the Nafion gate electrode. The fitted diffusion data are consistent with a P3HT doping region that is close to the source-drain electrodes rather than located at the P3HT:[Nafion:GOX] interface. Finally, we demonstrate that further improvements in sensor structure and morphology can be achieved by inkjet-printing the GOX layer, offering a pathway to low-cost printed biosensors for the detection of glucose in saliva.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Iterative spectral method for solving electrostatic or magnetostatic problems in complex and evolving heterostructures

Electrostatic or magnetostatic problems involving complex heterogeneity are nontrivial for modeling and simulation. Most existing numerical methods focus on sharp interface models and the computational cost increases with increasing complexity of the geometry. Here we develop an iterative spectral method, the bound charge successive approximation algorithm, to solve electrostatic or magnetostatic heterogeneity problems in the context of diffuse-interface modeling. As tests and verifications, this algorithm is applied to calculation of the depolarization factor of an ellipsoid, and simulation of random dielectric mixtures and the dielectophoretic motion of multiple particles. The algorithm shows excellent efficiency and the computational cost mainly depends on the permittivity or permeability contrast in the whole system, regardless of the complexity of the geometry. Specifically, for evolving heterostructures the solution of bound charge in one time step can be used as input for the next, which could further significantly shorten the iteration (approximation) process, making it practical to simulate the long-range electrostatic or magnetostatic interaction in complex and evolving heterostructures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Reproducing curvature effects due to differential diffusion in tabulated chemistry for premixed flames

Tabulated chemistry reduces the cost of modeling reacting flows by transporting fewer scalars. It also has the benefit of introducing fewer terms that must be closed for Reynolds-averaged Navier-Stokes and large eddy simulation modeling. In this study of premixed lean hydrogen–air flames, an extension to an existing tabulated chemistry model was developed by including differential diffusion and thermal diffusion of multiple species, using a mixture-averaged diffusion model. The transport equations for the tabulation variables (a progress variable and a mixture fraction-like variable) have been formulated, and curvature effects are intrinsic to the transport equations derived from first principles. The mathematical derivation is general, such that other fuel/air mixtures could be examined. The new tabulated chemistry model was then validated using a range of configurations, including flames with curvature, extinction, and hot spots. The model validation indicated excellent agreement with detailed chemistry simulations for lean premixed hydrogen–air mixtures. Conditional means of the progress variable source term were shown to agree well with detailed chemistry simulations of two-dimensional freely propagating flames and three-dimensional turbulent flames. Finally, progress variable source terms three times higher than one-dimensional flat flame results were identified and validated with detailed chemistry, and super-adiabatic hot spots were also predicted in the turbulent flame configuration.

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

A Route to Design Novel Functional Peptides by Applying a Denoising Diffusional Model to mRNA Display Libraries

In vitro directed evolution techniques, such as mRNA display, enable peptide ligand discovery and optimization. However, physical libraries that rely on a genetic code can only search a small fraction of sequence space due to inherent biases in the genetic code and experimental limitations. To address this challenge, denoising diffusion implicit models (DDIMs) are applied to generate novel peptide ligands against B‐cell lymphoma extra‐large (Bcl‐x L ), a key cancer target. Starting with high‐throughput sequencing data from previous selections, a DDIM is trained to produce novel sequences with high affinity binding. Experimental validation confirms that most generated sequences are functionally equivalent to the original library members for Bcl‐x L binding and demonstrated comparable binding kinetics and affinity relative to the wildtype and nearest original neighbors. Importantly, this approach generated rare sequences not easily accessible via mutation and directed evolution. These results indicate that DDIMs can complement and expand directed evolution data, efficiently exploring underrepresented regions of sequence space. This approach provides a broadly applicable framework for accelerating ligand discovery and optimizing molecular properties across diverse targets.

Qi, Pearl [Mork Family Department of Chemical Engi↗
DOE OSTI2010

FERMI LARGE AREA TELESCOPE FIRST SOURCE CATALOG

We present a catalog of high-energy gamma-ray sources detected by the Large Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray Space Telescope (Fermi), during the first 11 months of the science phase of the mission, which began on 2008 August 4. The First Fermi-LAT catalog (1FGL) contains 1451 sources detected and characterized in the 100 MeV to 100 GeV range. Source detection was based on the average flux over the 11 month period, and the threshold likelihood Test Statistic is 25, corresponding to a significance of just over 4{sigma}. The 1FGL catalog includes source location regions, defined in terms of elliptical fits to the 95% confidence regions and power-law spectral fits as well as flux measurements in five energy bands for each source. In addition, monthly light curves are provided. Using a protocol defined before launch we have tested for several populations of gamma-ray sources among the sources in the catalog. For individual LAT-detected sources we provide firm identifications or plausible associations with sources in other astronomical catalogs. Identifications are based on correlated variability with counterparts at other wavelengths, or on spin or orbital periodicity. For the catalogs and association criteria that we have selected, 630 of the sources are unassociated. Care was taken to characterize the sensitivity of the results to the model of interstellar diffuse gamma-ray emission used to model the bright foreground, with the result that 161 sources at low Galactic latitudes and toward bright local interstellar clouds are flagged as having properties that are strongly dependent on the model or as potentially being due to incorrectly modeled structure in the Galactic diffuse emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Criticality Calculations with MCNP6 - Practical Lectures

These slides are used to teach MCNP (Monte Carlo N-Particle) usage to nuclear criticality safety analysts. The following are the lecture topics: course information, introduction, MCNP basics, criticality calculations, advanced geometry, tallies, adjoint-weighted tallies and sensitivities, physics and nuclear data, parameter studies, NCS validation I, NCS validation II, NCS validation III, case study 1 - solution tanks, case study 2 - fuel vault, case study 3 - B&W core, case study 4 - simple TRIGA, case study 5 - fissile mat. vault, criticality accident alarm systems. After completion of this course, you should be able to: Develop an input model for MCNP; Describe how cross section data impact Monte Carlo and deterministic codes; Describe the importance of validation of computer codes and how it is accomplished; Describe the methodology supporting Monte Carlo codes and deterministic codes; Describe pitfalls of Monte Carlo calculations; Discuss the strengths and weaknesses of Monte Carlo and Discrete Ordinants codes; The diffusion theory model is not strictly valid for treating fissile systems in which neutron absorption, voids, and/or material boundaries are present. In the context of these limitations, identify a fissile system for which a diffusion theory solution would be adequate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Determining acceptable limits of fast-electron preheat in direct-drive-ignition–scale target designs

In direct-drive–ignition designs, preheat by fast electrons created by the two-plasmon–decay or stimulated Raman instabilities can increase the adiabat in the fuel layer and reduce compression and neutron yield. Since eliminating the preheat entirely is a major challenge, it is necessary to understand the levels of preheat that preclude ignition in a direct-drive target. Two 1-D ignition-scale target designs serve as the basis for examining the effects of synthetically increasing the levels of fast electrons using the 1-D radiation–hydrodynamic code LILAC , which includes two models of fast-electron transport. The first is an ignition design adapted from a 2-D polar-direct-drive (PDD) design for the National Ignition Facility (NIF). The second is derived from the first but uses a laser pulse that generates stronger shocks and a higher fuel adiabat. This more stable “alpha-burner” design approaches ignition and achieves yield multiplication as a result of alpha heating. The designs are then re optimized to recover performance. The igniting design, when fast-electron transport was modeled with diffusion, was found to tolerate 50% more fast-electron preheat of the cold (sub-50 eV) DT ice layer when the laser pulse was optimized using the optimizer Telios . When a straight-line fast-electron transport model was used, the effects of optimization were negligible. For the alpha-burner design, an increase of over a factor of at least 3 in the tolerable level of fast-electron preheat was obtained for both transport models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Methane production and emissions in trees and forests

Forest ecosystem methane (CH 4 ) research has focused on soils, but trees are also important sources and sinks in forest CH 4 budgets. Living and dead trees transport and emit CH 4 produced in soils; living trees and dead wood emit CH 4 produced inside trees by microorganisms; and trees produce CH 4 through an abiotic photochemical process. Here, we review the state of the science on the production, consumption, transport, and emission of CH 4 by living and dead trees, and the spatial and temporal dynamics of these processes across hydrologic gradients inclusive of wetland and upland ecosystems. Emerging research demonstrates that tree CH 4 emissions can significantly increase the source strength of wetland forests, and modestly decrease the sink strength of upland forests. Scaling from stem or leaf measurements to trees or forests is limited by knowledge of the mechanisms by which trees transport soil-produced CH 4 , microbial processes produce and oxidize CH 4 inside trees, a lack of mechanistic models, the diffuse nature of forest CH 4 fluxes, complex overlap between sources and sinks, and extreme variation across individuals. Understanding the complex processes that regulate CH 4 source–sink dynamics in trees and forests requires cross-disciplinary research and new conceptual models that transcend the traditional binary classification of wetland vs upland forest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Effects of Moisture on Diffusion in Unmodified Wood Cell Walls: A Phenomenological Polymer Science Approach

Despite the importance of cell wall diffusion to nearly all aspects of wood utilization, diffusion mechanisms and the detailed effects of moisture remain poorly understood. In this perspective, we introduce and employ approaches established in polymer science to develop a phenomenological framework for understanding the effects of moisture on diffusion in unmodified wood cell walls. The premise for applying this polymer-science-based approach to wood is that wood polymers (cellulose, hemicelluloses, and lignin) behave like typical solid polymers. Therefore, the movement of chemicals through wood cell walls is a diffusion process through a solid polymer, which is in contrast to previous assertions that transport of some chemicals occurs via aqueous pathways in the cell wall layers. Diffusion in polymers depends on the interrelations between free volume in the polymer matrix, molecular motions of the polymer, diffusant dimensions, and solubility of the diffusant in the polymer matrix. Because diffusion strongly depends on whether a polymer is in a rigid glassy state or soft rubbery state, it is important to understand glass transitions in the amorphous wood polymers. Through a review and analysis of available literature, we conclude that in wood both lignin and the amorphous polysaccharides very likely have glass transitions. After developing and presenting this polymer-science-based perspective of diffusion through unmodified wood cell walls, suggested directions for future research are discussed. A key consideration is that a large difference between diffusion through wood polymers and typical polymers is the high swelling pressures that can develop in unmodified wood cell walls. This pressure likely arises from the hierarchical structure of wood and should be taken into consideration in the development of predictive models for diffusion in unmodified wood cell walls.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Model-Independent Characterization of Charge Diffusion in Thick Fully Depleted CCDs

We present a new method to measure charge diffusion in charge-coupled devices (CCDs). The method is based on a statistical characterization of the shapes of charge clouds produced by low-energy X-rays using known properties of the two-dimensional Gaussian point-spread function (PSF). The algorithm produces reliable upper and lower bounds on the size of the PSF for photons converting near the entrance window of a device. It is optimally suited to the case of thick back-illuminated CCDs where the X-ray absorption length is smaller than the silicon thickness and the diffusion scale is of the same order as the pixel size. The only assumptions are that the charge cloud width is a monotonically increasing function of distance from the conversion point to the buried channel, and that the conversion probability is peaked at the surface. Otherwise, no physical models of carrier transport or knowledge of the electric field profile in the CCD are needed. In suboptimal conditions, the upper bound increases and the lower bound is unaffected, so confidence in the correctness of results is retained. The new method has been benchmarked against Monte Carlo simulations and tested on X-ray images measured on thick high-resistivity prototype CCDs for the Large Synoptic Survey Telescope. In Monte Carlo simulations of noiseless images having the optimal diffusion scale, the upper bound approximated the true PSF within 5%, increasing to 10% in simulations with added noise. Even with severely undersampled or truncated PSFs, the method brackets the true value to within 25%. Our method is accurate and computationally efficient and offers a fast and simple experimental setup.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗