Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equations”

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,261 records · Page 70

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 ↗

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↗

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↗

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↗

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↗

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↗

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↗

Effects of q-profile structure on turbulence spreading: A fluctuation intensity transport analysis

This paper studies effects of q-profile structure on turbulence spreading. It reports results of numerical experiments using global gyrokinetic simulations. We examine propagation of turbulence, triggered by an identical linear instability in a source region, into an adjacent, linearly stable region with variable q-profile. The numerical experiments are designed so as to separate the physics of turbulence spreading from that of linear stability. The strength of turbulence spreading is measured by the penetration depth of turbulence. Dynamics of spreading are elucidated by fluctuation intensity balance analysis, using a model intensity evolution equation which retains nonlinear diffusion and damping, and linear growth. It is found that turbulence spreading is strongly affected by magnetic shear s, but is hardly altered by the safety factor q itself. There is an optimal range of modest magnetic shear which maximizes turbulence spreading. For high to modest shear values, the spreading is enhanced by the increase of the mode correlation length with decreasing magnetic shear. However, the efficiency of spreading drops for sufficiently low magnetic shear even though the mode correlation length is comparable to that for the case of optimal magnetic shear. The reduction of spreading is attributed to the increase in time required for the requisite nonlinear mode-mode interactions. The effect of increased interaction time dominates that of increased mode correlation length. Our findings of the reduction of spreading and the increase in interaction time at weak magnetic shear are consistent with the well-known benefit of weak or reversed magnetic shear for core confinement enhancement. Weak shear is shown to promote locality, as well as stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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↗

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↗

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↗

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↗

SU-D-210-03: Limited-View Multi-Source Quantitative Photoacoustic Tomography

Purpose: This work is to investigate a novel limited-view multi-source acquisition scheme for the direct and simultaneous reconstruction of optical coefficients in quantitative photoacoustic tomography (QPAT), which has potentially improved signal-to-noise ratio and reduced data acquisition time. Methods: Conventional QPAT is often considered in two steps: first to reconstruct the initial acoustic pressure from the full-view ultrasonic data after each optical illumination, and then to quantitatively reconstruct optical coefficients (e.g., absorption and scattering coefficients) from the initial acoustic pressure, using multi-source or multi-wavelength scheme.Based on a novel limited-view multi-source scheme here, We have to consider the direct reconstruction of optical coefficients from the ultrasonic data, since the initial acoustic pressure can no longer be reconstructed as an intermediate variable due to the incomplete acoustic data in the proposed limited-view scheme. In this work, based on a coupled photo-acoustic forward model combining diffusion approximation and wave equation, we develop a limited-memory Quasi-Newton method (LBFGS) for image reconstruction that utilizes the adjoint forward problem for fast computation of gradients. Furthermore, the tensor framelet sparsity is utilized to improve the image reconstruction which is solved by Alternative Direction Method of Multipliers (ADMM). Results: The simulation was performed on a modified Shepp-Logan phantom to validate the feasibility of the proposed limited-view scheme and its corresponding image reconstruction algorithms. Conclusion: A limited-view multi-source QPAT scheme is proposed, i.e., the partial-view acoustic data acquisition accompanying each optical illumination, and then the simultaneous rotations of both optical sources and ultrasonic detectors for next optical illumination. Moreover, LBFGS and ADMM algorithms are developed for the direct reconstruction of optical coefficients from the acoustic data. Jing Feng and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000) and the Shanghai Pujiang Talent Program (#14PJ1404500)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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↗

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↗

Snow Micro-Structure Model

PIKA is a MOOSE-based application for modeling micro-structure evolution of seasonal snow. The model will be useful for environmental, atmospheric, and climate scientists. Possible applications include application to energy balance models, ice sheet modeling, and avalanche forecasting. The model implements physics from published, peer-reviewed articles. The main purpose is to foster university and laboratory collaboration to build a larger multi-scale snow model using MOOSE. The main feature of the code is that it is implemented using the MOOSE framework, thus making features such as multiphysics coupling, adaptive mesh refinement, and parallel scalability native to the application. PIKA implements three equations: the phase-field equation for tracking the evolution of the ice-air interface within seasonal snow at the grain-scale; the heat equation for computing the temperature of both the ice and air within the snow; and the mass transport equation for monitoring the diffusion of water vapor in the pore space of the snow.

Micah Johnson, AndrewSlaughter↗