Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion equation”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A computational thrombus formation model: application to an idealized two-dimensional aneurysm treated with bare metal coils

Cardiovascular implantable devices alter the biofluid dynamics and biochemistry of the blood in which they are placed. These perturbations can lead to thrombus formation which may or may not be desired, depending on the application. In this work, a computational model is developed that couples biofluid dynamics and biochemistry to predict the clotting response of blood to such devices. The model consists of 28 advection–diffusion–reaction partial differential equations to track proteins in the blood involved in clotting and utilizes boundary flux terms to model the initiation of the intrinsic clotting pathway at thrombogenic device surfaces. We use this model to simulate the transient clot growth within a 2D idealized bifurcation aneurysm filled with various distributions of bare metal coils with similar packing densities. The clot model predicts initial clot formation to occur in areas along coil surfaces where flow is minimal and where time-averaged shear rates are the smallest. Among the six coil-filled aneurysm cases simulated, maximum thrombus occlusion ranged between 80.8 and 92.2% of the post-treatment aneurysm volume and was achieved 325–450 s after treatment. Finally, with further refinement and validation, the computational clotting model will be a valuable engineering tool for evaluating and comparing the relative performance of cardiovascular implantable devices.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Continuous Fluorescence Depletion Anisotropy Measurement of Protein Rotation

Protein rotation in viscous environments can be measured by fluorescence depletion anisotropy (FDA) which combines long lifetimes of chromophore triplet states with the sensitivity of fluorescence excitation and detection. FDA achieves sensitivity well beyond that attainable by the more common technique of time-resolved phosphorescence anisotropy (TPA). We have now combined benefits of both time-domain and frequency-domain FDA into a single continuous technique (CFDA). Intensity and polarization of a single laser beam are modulated continuously according to a complex, repeating waveform. Fluorescence signals excited from triplet-forming fluorescent probes are digitized over recurring waveform periods by a high-speed signal averager. CFDA experiments typically involve substantial ground state depletion. Thus signals, unlike those of TPA, are not linear in the exciting light intensity and simple data analysis based on such linearity is not appropriate. An exact solution of the coupled diffusion and triplet production/decay equation describing CFDA within individual data points has been combined with simulated annealing optimization to extract triplet and anisotropy decay kinetics from experimental data. Related calculations compare possible excitation waveforms with respect to rotational information provided per fluorescence photon. We present CFDA results for the model system of eosin conjugates of carbonic anhydrase, BSA and immunoglobulin G in 90% glycerol at various temperatures and initial cellular results on eosin-IgE bound to 2H3 cell Type I Fcε receptors. We explore how CFDA reflects rotational parameters of heterogeneous systems and discuss challenges of extending this method to single cell microscopic measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Experimental and Modeling Studies for the Removal of Crystal Violet Dye from Aqueous Solutions using Eco-friendly Gracilaria corticata Seaweed Activated Carbon/Zn/Alginate Polymeric Composite Beads

This research article describes, an eco-friendly activated carbon prepared from the Gracilaria corticata seaweeds which was employed for the preparation of biodegradable polymeric beads for the efficient removal of crystal violet dye in an aqueous solution. The presence of chemical functional groups in the adsorbent material was detected using FTIR spectroscopy. The morphology and physical phases of the adsorbent materials were analyzed using SEM and XRD studies respectively. Batch mode dye adsorption behavior of the activated carbon/Zn/alginate polymeric beads was investigated as a function of dosage, solution pH, contact time, initial dye concentration and temperature. Maximum dye removal was observed at a pH of 5.0, 1 g of adsorbent dosage with 60 mg/L dye concentration, 50 min of contact time and at 30 °C. The equilibrium modeling studies were analyzed with Freundlich and Langmuir adsorption isotherms and the adsorption dynamics was predicted with Lagergren’s pseudo-first order, pseudo-second order equations and intra particle diffusion models. The process of dye removal followed a pseudo second-order kinetics rather than pseudo first order. The thermodynamic parameters like standard Gibbs free energy change (ΔG°), enthalpy change (ΔH°) and entropy change (ΔS°) were determined and the results imply that the adsorption process was spontaneous, endothermic and increases the randomness between the adsorbent and adsorbate. The results from the experimental and correlation data reveal that the Gracilaria corticata activated carbon/Zn/alginate polymeric beads have proved to be an excellent adsorbent material for the removal of CV dye.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Diffusion-driven GaInP/GaAs light-emitting diodes enhanced by modulation doping

Diffusion-driven charge transport (DDCT) in III–V light-emitting diodes (LEDs) can enable unconventional optoelectronic devices and functionality by fundamentally changing device design and the current injection principle. In our recent study, an AlGaAs/GaAs DDCT–LED consisting of an array of lateral heterojunctions was studied for large-area applications at high powers. Here, we investigate the current spreading and recombination uniformity of a modulation doped GaInP/GaAs DDCT–LED. In particular, we analyze how the background doping of the lower GaInP cladding layer (CL) and the GaAs substrate changes the carrier distribution within the active region of the device. Our charge transport simulations based on the drift-diffusion current and continuity equations predict that modulation doping by a p-doped CL provides much higher recombination uniformity at high powers compared to an n-doped CL. Most importantly, improved current spreading is achieved while maintaining excellent device performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Exploring the impact of hydro-mechanical coupling on compaction localization

In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.

58 GEOSCIENCES↗
DOE OSTI2022

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.

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

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2015

Modeling of electron energy spectra and mobilities in semi-metallic Hg{sub 1−x}Cd{sub x}Te quantum wells

Electron mobility, energy spectra, and intrinsic carrier concentrations in the n-type Hg{sub 0.32}Cd{sub 0.68}Te/Hg{sub 1−x}Cd{sub x}Te/Hg{sub 0.32}Cd{sub 0.68}Te quantum well (QW) in semi-metallic state are numerically modeled. Energy spectra and wave functions were calculated in the framework of the 8-band k-p Hamiltonian. In our model, electron scattering on longitudinal optical phonons, charged impurities, and holes has been taken into account, and the mobility has been calculated by an iterative solution of the Boltzmann transport equation. Our results show that the increase of the electron concentration in the well enhances the screening of the 2D electron gas, decreases the hole concentration, and can ultimately lead to a high electron mobility at liquid nitrogen temperatures. The increase of the electron concentration in the QW could be achieved in situ by delta-doping of barriers or by applying the top-gate potential. Our modeling has shown that for low molar composition x the concentration of holes in the well is high in a wide range of electron concentrations; in this case, the purity of samples does not significantly influence the electron mobility. These results are important in the context of establishing optimal parameters for the fabrication of high-mobility Hg{sub 1−x}Cd{sub x}Te quantum wells able to operate at liquid nitrogen temperature and thus suitable for applications in terahertz detectors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Broadband light sources based on InAs/InGaAs metamorphic quantum dots

We propose a design for a semiconductor structure emitting broadband light in the infrared, based on InAs quantum dots (QDs) embedded into a metamorphic step-graded In{sub x}Ga{sub 1−x}As buffer. We developed a model to calculate the metamorphic QD energy levels based on the realistic QD parameters and on the strain-dependent material properties; we validated the results of simulations by comparison with the experimental values. On this basis, we designed a p-i-n heterostructure with a graded index profile toward the realization of an electrically pumped guided wave device. This has been done by adding layers where QDs are embedded in In{sub x}Al{sub y}Ga{sub 1−x−y}As layers, to obtain a symmetric structure from a band profile point of view. To assess the room temperature electro-luminescence emission spectrum under realistic electrical injection conditions, we performed device-level simulations based on a coupled drift-diffusion and QD rate equation model. On the basis of the device simulation results, we conclude that the present proposal is a viable option to realize broadband light-emitting devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Conforming and nonconforming virtual element methods for elliptic problems

Here we present, in a unified framework, new conforming and nonconforming virtual element methods for general second-order elliptic problems in two and three dimensions. The differential operator is split into its symmetric and nonsymmetric parts and conditions for stability and accuracy on their discrete counterparts are established. These conditions are shown to lead to optimal H 1 - and L 2 -error estimates, confirmed by numerical experiments on a set of polygonal meshes. The accuracy of the numerical approximation provided by the two methods is shown to be comparable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Anharmonic phonon damping enhances the T c of BCS-type superconductors

In this study, a theory of superconductivity is presented where the effect of anharmonicity, as entailed in the acoustic, or optical, phonon damping, is explicitly considered in the pairing mechanism. The gap equation is solved including diffusive Akhiezer damping for longitudinal acoustic phonons or Klemens damping for optical phonons, with a damping coefficient which, in either case, can be directly related to the Grüneisen parameter and hence to the anharmonic coefficients in the interatomic potential. The results show that the increase of anharmonicity has a strikingly nonmonotonic effect on the critical temperature T c . The optimal damping coefficient yielding maximum T c is set by the velocity of the bosonic mediator. This theory may open up unprecedented opportunities for material design where T c may be tuned via the anharmonicity of the interatomic potential, and presents implications for the superconductivity in the recently discovered hydrides, where anharmonicity is very strong and for which the anharmonic damping is especially relevant.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Fast Multigrid Reduction-in-Time for Advection via Modified Semi-Lagrangian Coarse-Grid Operators

Many iterative parallel-in-time algorithms have been shown to be highly efficient for diffusion-dominated partial differential equations (PDEs) but are inefficient or even divergent when applied to advection-dominated PDEs. We consider the application of the multigrid reduction-in-time (MGRIT) algorithm to linear advection PDEs. Here, the key to efficient time integration with this method is using a coarse-grid operator that provides a sufficiently accurate approximation to the so-called ideal coarse-grid operator. For certain classes of semi-Lagrangian discretizations, we present a novel semi-Lagrangian-based coarse-grid operator that leads to fast and scalable multilevel time integration of linear advection PDEs. The coarse-grid operator is composed of a semi-Lagrangian discretization followed by a correction term, with the correction designed so that the leading-order truncation error of the composite operator is approximately equal to that of the ideal coarse-grid operator. Parallel results show substantial speed-ups over sequential time integration for variable-wave-speed advection problems in one and two spatial dimensions, and using high-order discretizations up to order five. The proposed approach establishes the first practical method that provides small and scalable MGRIT iteration counts for advection problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Electrochemical Reactivity of Hydrogen Storage Materials: Exploring Borohydride and Hydrazinium Salt Mixtures

Electrochemical characterization of hydrogen storage materials was conducted in a non-aqueous environment to investigate the direct electrochemical release and consumption of hydrogen and the potential for regeneration. We first address the challenge of minimal solubility of the synthetic precursors, sodium borohydride (NaBH 4 ) and hydrazinium bromide (N 2 H 5 Br), in both organic and inorganic solvents. We next determine and calibrate a reference electrode formulation compatible with our non-aqueous media and analytes that demonstrates a stable reference potential. We employ cyclic voltammetry (CV) to characterize the precursors and mixtures thereof. Each CV peak is assigned to a corresponding electrochemical reaction. Using the rate-dependent CV method and Randles–Ševčík equation, we calculate the diffusion coefficient of each chemical (NaBH 4 and N 2 H 5 Br). Analysis of the CVs, coupled with 11B NMR analysis, reveals a room temperature chemical transformation of NaBH 4 and N 2 H 5 Br mixtures into hydrazine borane (N 2 H 4 BH 3 ). These results are particularly significant, considering the limited information available on the electrochemical characterization of metal borohydride and hydrazinium salt in non-aqueous media. This work establishes a foundation for adapting a non-aqueous electrochemical system to further study the borohydride family of chemistries and to design and develop electrochemical devices for direct electrical and chemical energy interconversion with hydrogen storage materials.

08 HYDROGEN↗
DOE OSTI2021

TMAP7

TMAP7 was developed by the Fusion Safety Program at INL primarily as a safety analysis code, to analyze tritium retention and loss in fusion reactor structures and systems during normal operation and accident conditions. TMAP7 solves one-dimensional thermal- and mass-diffusive transport and trapping equations for transport through structures, across the interface between enclosures and structures, and zero dimensional fluid transport between enclosures.

Humrickhouse, PaulW.↗
DOE OSTI2021

TMAP4

TMAP4 was developed by the Fusion Safety Program at INL primarily as a safety analysis code, to analyze tritium retention and loss in fusion reactor structures and systems during normal operation and accident conditions. TMAP4 solves one-dimensional thermal- and mass-diffusive transport and trapping equations for transport through structures, across the interface between enclosures and structures, and zero dimensional fluid transport between enclosures.

Humrickhouse, PaulW.↗